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Abstract 


Because  of  their  high  strength  and  stiffness  together  with  their  light  weight,  fiber 
reinforced  composite  materials  offer  great  potential  for  applications  particularly  in  the 
aerospace  industry.  The  weight  savings  translate  into  increased  performance  and 
decreased  fuel  costs.  In  addition,  the  use  of  these  materials  avoids  the  dependence  on 
foreign  sources  for  the  critical  elements  needed  in  the  new  exotic  metallic  alloys  that 
otherwise  might  be  used.  Early  research  has  proven  the  usefulness  of  these  materials, 
but  the  need  to  quantitatively  characterize  important  material  properties  and  develop 
applicable  nondestructive  evaluation  techniques  remains.  Important  physical  properties 
(i.e.  mechanical,  thermal,  electrical)  need  to  be  measured.  Moreover,  relationships 
between  these  physical  properties  and  important  engineering  properties  such  as 
strength,  residual  strength  after  impact  and  fatigue  loading,  and  fiber-matrix  interfacial 
strength  need  to  be  examined  to  provide  a  basis  for  quantitative  nondestructive  evalua¬ 
tion  of  these  materials. 

Toward  this  goal,  linear  and  nonlinear  elastic  properties  have  been  demonstrated 
to  be  important  physical  properties  in  conventional  materials.  In  particular,  nonlinear 
properties  are  important  in  the  nondestructive  determination  of  applied  and  residual 
stress  (strain)  as  well  as  measuring  the  interatomic  bonding  forces  in  crystalline  solids. 
Also,  several  investigations  have  established  a  possible  relationship  between  nonlinear 
elastic  properties  and  ultimate  strength  in  aluminum  and  carbon  steel. 

This  work  presents  the  theoretical  treatment  of  linear  and  nonlinear  elasticity  in 
a  unidirectionally  fiber  reinforced  composite  as  well  as  measurements  for  a  unidirec¬ 
tional  graphite/epoxy  composite  (T300/5208).  Linear  elastic  properties  were  measured 
by  both  ultrasonic  and  strain  gauge  measurements.  The  nonlinear  properties  were 
determined  by  measuring  changes  in  ultrasonic  "natural"  phase  velocity  with  a  pulsed 
phase  locked  loop  interferometer  as  a  function  of  stress  and  temperature.  These  meas¬ 
urements  provide  the  basis  for  further  investigations  into  the  relationship  between 
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nonlinear  elastic  properties  and  other  important  properties  such  as  strength  and  fiber- 
matrix  interfacial  strength  in  graphite/epoxy  composites. 
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L  INTRODUCTION 


The  theory  of  nonlinear  elasticity  and  the  measurement  of  nonlinear  elastic  pro¬ 
perties  of  materials  have  been  subjects  of  great  interest  for  many  years.  Nonlinear  elas¬ 
ticity  is  important  in  many  areas  of  physical  measurements  including  the  characteriza¬ 
tion  of  interatomic  bonding  forces  and  anharmonic  behavior  of  a  crystal  lattice,  nondes¬ 
tructive  measurement  of  applied  and  residual  stress  (strain)  using  ultrasonic  techniques, 
and  stress-strain  relations  in  finite  strain  of  materials.  This  will  be  discussed  in  more 
detail  later. 

One  of  the  first  successful  theoretical  treatments  of  the  subject  of  nonlinear  elas¬ 
ticity  was  presented  by  Murnaghan  [l]  in  1937.  He  extended  the  theory  of  infinitesimal 
elasticity  to  the  general  case  including  strains  of  any  magnitude.  Birch  [2],  in  1938, 
pointed  out  the  importance  of  this  theory  for  the  case  of  seismic  waves  superimposed  on 
the  large  hydrostatic  pressures  within  the  earth.  Since  that  time  a  number  of  research¬ 
ers  have  presented  volumes  of  material  on  the  subject  of  nonlinear  elasticity.  Although 
a  detailed  review  of  this  information  is  beyond  the  scope  of  this  essay,  important  work 
in  several  areas  will  be  discussed.  First,  literature  discussing  the  basic  property  of  non¬ 
linear  elasticity  and  its  effects  on  ultrasonic  wave  propagation  will  be  reviewed.  This 
will  be  followed  by  a  discussion  of  some  of  the  literature  involving  the  development  of 
theoretical  and  measurement  techniques  as  well  as  some  of  the  applications  for  non¬ 
linear  elastic  research.  Also,  previous  measurements  of  the  nonlinear  properties  of  poly¬ 
mers  and  composites  will  be  discussed. 

The  study  of  the  mechanical  properties  of  materials  is  concerned  with  the  defor¬ 
mation  of  a  material  upon  application  of  an  external  force  as  well  as  its  behavior  upon 
removal  of  that  force.  There  are  three  types  of  mechanical  behavior  exhibited  by 
materials  under  loading.  These  are  elastic,  anelastic  and  plastic.  In  elastic  deforma¬ 
tion,  the  material  deforms  instantaneously  upon  application  of  load  and  returns  to  its 
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oripnal  shape  instantaneously  upon  removal  of  load.  This  behavior  can  be  either  linear 
or  nonlinear.  This  is  demonstrated  in  Figs,  1.1  (a)  and  1.1  (b)  for  the  simple  case  of 
uniamal  loading  v/here  the  elongation  seir/es  as  a  measure  of  the  deformation  of  the 
material.  In  anelastic  deformation^  the  material  deforms  over  time  upon  the  application 
of  load  and  returns  to  its  original  shape  over  time  upon  the  removal  of  load.  This  is 
demonstrated  in  Fig.  1.2.  ¥/lien  a  material  undergoes  plastic  deformatioBj  it  deforms 
instantaneously  upon  application  of  load  but  does  not  return  to  its  original  shape  upon 
the  removal  of  the  load.  Fig.  1.3  shows  the  load^elongation  characteristics  of  plastic 
deformation. 

Hov/ever,  the  behavior  of  real  materials  can  be  a  combination  of  these  three 
types  of  mechanical  behavior  and  is  generally  dependent  on  the  amount  of  load  applied. 
The  material  may  be  elastic  for  small  loads  but  the  deformation  may  be  anelastic  or 
plastic  at  higher  loads.  Also  changes  in  other  parameters  such  as  temperature  may 
affect  the  response  of  the  material  upon  loading.  The  present  work  is  interested  in  the 
case  of  nonlinear  elastic  behavior  where  the  deformation  is  recoverable  but  not  linearly 
related  to  the  applied  load. 

/m  understanding  of  linear  and  nonlinear  elasticity  can  be  gained  by  considering 
theories  of  interatomic  bonding  in  molecules  such  as  the  Born  [3]  and  Madelung  [4] 
theories  for  ionic  crystals.  In  these  theories,  the  total  energy  of  bonding  is  the  energy 
due  to  the  sum  of  the  attractive  and  repulsive  forces  between  the  atoms  multiplied  by 
the  distance  over  v/hich  they  act.  The  general  shape  of  such  a  curve  is  shown  in  Fig.  1.4 
v/ith  the  potential  energy  (E)  versus  the  separation  distance  (r)  plotted.  The  equili¬ 
brium  separation  distance  is  Tq  where  the  potential  energy  is  at  a  minimum  Eq  .  To 
move  the  atoms  closer  together  or  further  apart  raises  the  energy  of  the  system  and 
thus  requires  the  application  of  an  external  force.  This  force  can  be  determined  by  tak¬ 
ing  the  derivative  of  E  v/ith  respect  to  separation  distance  for  the  curve  in  Fig.  1.4.  The 
resulting  force  versus  separation  distance  curve  is  then  a  simple  analogy  on  the  atomic 


ELONGATION 

a) 

Figure  1.1  -  (a)  Load-Elongation  curve  for  linear  elastic 
deformation. 


b) 


Figure  1.1  -  (b)  Load-Elongation  curve  for  nonlinear  elast 
deformation. 
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scale  to  the  load-displacement  curves  discussed  earlier. 

For  small  displacements  about  ro  ,  the  potential  energy  curve  appears  to  be  and 
can  be  approximated  as  a  quadratic  function.  Therefore  its  derivative  is  a  linear  curve. 
That  is,  the  relationship  between  force  and  separation  distance  (load  and  elongation)  is 
linear  for  small  displacements  about  equilibrium.  However,  as  the  amplitude  of  the  dis¬ 
placement  becomes  larger,  the  potential  energy  curve  is  no  longer  quadratic.  This 
implies  that  the  force-separation  curve  is  no  longer  linear  which  gives  rise  to  nonlinear 
elasticity. 

Another  interesting  point  about  Fig.  1.4  is  the  asymmetry  of  the  potential  energy 
curve  about  ro  .  This  provides  a  simple  explanation  for  the  understanding  of  thermal 
expansion  in  materials.  Since  the  curve  is  asymmetric,  as  the  material  is  heated  and 
the  atomic  lattice  vibrations  increase,  the  mean  separation  distance  between  atoms 
changes.  For  the  curve  in  Fig.  1.4,  the  mean  separation  distance  increases  with  increas¬ 
ing  amplitude  vibrations  about  Tq  .  Thus,  in  agreement  with  empirical  observations  in 
most  materials,  expansion  occurs  as  the  material  is  heated. 

Because  of  the  small  strains  (on  the  order  of  10“®  )  imposed  on  the  sample  during 
ultrasonic  wave  propagation,  nonlinear  elastic  behavior  is  generally  neglected  in  the 
theories  of  the  propagation  of  ultrasonic  waves.  However,  as  pointed  out  by  Green  [5], 
there  are  three  cases  where  nonlinear  elasticity  may  become  important  in  elastic  wave 
propagation.  First,  the  amplitude  of  the  wave  may  become  large  enough  to  cause  finite 
strains  in  the  material.  In  the  second  case,  nonlinear  behavior  may  occur  when  a  small 
amplitude  wave  is  superimposed  on  a  large  external  static  stress.  Also,  nonlinear  effects 
may  be  caused  by  defects  in  the  material  which  cause  localized  regions  of  finite  strain. 

The  effects  of  nonlinear  elasticity  on  ultrasonic  wave  propagation  were  also 
pointed  out  by  Green  [5].  These  include  the  distortion  of  finite  amplitude  sinusoidal 
longitudinal  waves.  As  they  propagate,  energy  from  the  fundamental  frequency  is 
transferred  into  the  harmonics  that  are  generated.  The  effect  has  been  used  by 
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numerous  investigators  to  measure  the  nonlinear  properties  of  solids.  Measurements  on 
aluminum  single  crystals  as  well  as  NaCl,  KCl,  LiF  and  a  magnesium  aluminum  alloy 
were  presented  by  Gedroils  and  Krasil’Nilcov  [6]  in  1963.  Breazeale  and  Thompson  [7] 
reported  measurements  of  harmonic  distortion  in  poly  crystalline  aluminum  also  in  1963 
which  was  followed  by  numerous  other  investigations. 

Another  effect  of  nonlinear  elasticity  on  ultrasonic  wave  propagation  is  that  a 
pure  mode  nonlinear  longitudinal  wave  may  propagate  alone,  but  a  nonlinear  transverse 
wave  must  have  a  longitudinal  component.  Also,  nonlinear  transverse  waves  do  not  dis¬ 
tort  when  propagating  in  a  defect  free  solid. 

Nonlinear  elastic  waves  can  interact  with  each  other  in  a  solid  to  produce  other 
waves.  This  interaction  can  also  take  place  with  thermal  phonons  to  cause  energy  loss 
from  the  elastic  wave.  The  amount  of  interaction  depends  on  the  amplitude  of  the  elas¬ 
tic  waves.  Interactions  of  two  nonlinear  ultrasonic  waves  to  produce  a  third  ultrasonic 
wave  were  demonstrated  by  Rollins  [8]  in  1963  in  fused  silica,  polycrystalline  aluminum 
and  polycrystalline  magnesium.  Rollins,  Taylor  and  Todd  [9]  demonstrated  this 
phenomenon  again  in  1964  and  reported  the  measurements  to  be  correct  to  within  an 
order  of  magnitude  with  predictions  based  on  the  measured  third  order  elastic  constants 
of  polycrystalline  magnesium.  Since  this  time,  further  examinations  of  nonlinear  wave 
interactions  have  followed. 

Another  effect  of  nonlinearity  in  elastic  wave  propagation  is  that  the  velocity  of 
a  small  amplitude  ultrasonic  v/ave  superimposed  on  a  large  static  stress  is  dependent  on 
the  amount  of  static  stress  applied.  It  is  this  effect  which  provides  the  basis  for  the 
measui  ement  technique  used  in  this  research  and  will  thus  be  discussed  in  more  detail. 
The  effect  can  be  accounted  for  by  two  simple  theoretical  explanations.  In  the  first,  the 
equation  for  the  velocity  of  an  elastic  wave  propagating  along  a  long  thin  rod  is  recalled 
to  be 


(1.1) 


-  9- 


where  v  is  the  ultrasonic  phase  velocity, 

E  is  Young’s  modulus 
and  p  is  the  density. 

This  is  based  on  the  assumption  that  the  material  is  homogeneous,  isotropic  and 
behaves  in  a  linear  elastic  fashion.  These  assumptions  imply  that  the  constitutive  equa¬ 
tion  of  the  material  has  the  form 

(7  =  E€  (1.2) 

where  c  is  the  stress 
and  e  is  the  strain. 

If  a  nonlinear  constitutive  equation  is  used  such  as 

O’  =  Ki6  +  £2^^,  (1-3) 

it  follows  that  the  velocity  will  now  depend  on  strain  and  therefore  stress. 

An  alternative  approach  to  understanding  this  is  based  on  the  potential  energy 
theory  for  interatomic  bonding  discussed  earlier.  Since  the  velocity  of  ultrasonic  waves, 
together  with  the  density  of  the  medium  give  a  measure  of  the  modulus  (second  deriva¬ 
tive  of  the  potential  energy),  this  modulus  can  be  evaluated  as  a  function  of  strain  (or 
stress)  by  measuring  the  velocity  as  a  function  of  strain  (or  stress).  Also,  since  for  large 
displacements  (i.e.  finite  strain),  the  potential  energy  function  is  not  a  simple  quadratic, 
the  modulus  and  thus  the  velocity  would  be  expected  to  change  with  applied  stress.  As 
will  be  seen  later,  the  relationship  between  velocity  and  stress  should  be  linear  if  terms 
up  to  third  order  in  strain  are  included  in  the  elastic  theory. 

The  theoretical  development  of  the  equations  relating  ultrasonic  velocity  to 
stress  began  with  the  original  work  on  finite  deformation  of  materials  by  Murnaghan  [l]. 
This  work  provided  the  basis  for  the  theory  presented  by  Hughes  and  Kelly  [lO]  in  1953 
in  which  they  derived  wave  speeds  as  a  function  of  stress  in  homogeneous,  isotropic 
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materials.  They  considered  both  uniaxial  and  hydrostatic  compression  for  different  com¬ 
binations  of  longitudinal  and  transverse  (shear)  mode  waves.  The  following  seven  equa¬ 
tions  v;ere  presented: 

=  X  -}-  2^  —  -“{61  +  4m  7\  -!-  10/i)  (1.3) 

3Kq 

P  1 

=  P-  —  (3m  —  — n  -h  SX  4-  6/.()  (1.4) 

OIVq  2 

WVx  ~  ^  (2!  -!-  X  -|-  ^  (4m  -f-  4X  4-  10/i))  (1.5) 

oAq  p 

P  Xn 

o  P  9\ 

PoWy  =  X  2/i  —  (21 - (m  -1-  X  4-  2p))  (1.7) 

oi\..Q  /i 

9  P  Xp 

Pop-  =  p-  ^  Ti"  ^ 

^  P  r  /X-h/i> 


v/here  X,/i  are  the  second  order  linear  elastic  moduli  of  Lame’  for  a  homogeneous, 
isotropic  m-aterial, 

l,m,n  are  the  third  order  elastic  moduli  for  a  homogeneous,  isotropic 
material  as  defined  by  Murnaghan, 
pQ  is  the  density  of  the  unstressed  material, 

Kg  is  the  bulk  modulus, 

P  is  the  uniaxial  or  hydrostatic  compression, 
and  V  is  the  velocity  of  an  ultrasonic  wave  propagating  along  the  x  axis.  The 
nrst  subscript  on  v  refers  to  the  direction  of  polarization,  while  the 
second  gives  the  direction  of  loading.  0  for  the  second  subscript  implies 
hydrostatic  loading. 
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They  made  measurements  of  the  nonlinear  constants  of  polystyrene,  Axmco  Iron  and 
Pyrex.  Toupin  and  Bernstein  [ll]  followed  this  work  with  an  extension  of  the  theory  to 
include  a  material  of  arbitrary  crystal  symmetry.  This  work  was  published  in  1961. 

The  next  major  theoretical  developments  were  presented  in  a  series  of  papers  by 
Brugger  [12,13]  and  Thurston  and  Brugger  [14|  in  1964  and  1965.  These  papers  provided 
the  formulation  and  notation  for  most  of  the  research  presented  in  following  years  on 
the  topic.  In  the  first  paper  by  Brugger  [12],  a  formal  thermodynamic  definition  of  the 
higher  order  elastic  constants  of  a  solid  was  given.  The  relationships  between  his 
definition  and  those  used  by  Nlurnaghan  [l]  and  other  authors  were  also  presented.  In 
the  paper  of  Thurston  and  Brugger  [14],  the  concept  of  an  ultrasonic  "natural"  velocity 
was  introduced.  The  "natural"  velocity  was  defined  as  the  initial  (unstressed)  path 
length  divided  by  the  time  of  flight  of  the  elastic  wave.  The  concept  of  "natural"  velo¬ 
city  is  of  great  importance  in  the  measurement  of  nonlinear  properties  because  it 
requires  that  the  change  in  only  one  parameter,  the  time  of  flight,  be  monitored  as  a 
function  of  stress  in  the  experiments.  In  determining  normal  ultrasonic  phase  velocity 
changes,  the  time  of  flight  as  well  as  the  propagation  distance  must  be  measured  mak¬ 
ing  the  experimental  measurements  more  difficult  and  uncertain.  Also  presented  in  this 
paper  were  general  equations  for  the  relationship  between  the  natural  velocity 
changes  and  applied  stress  for  different  crystal  symmetries,  wave  modes,  and  stress  con¬ 
ditions.  In  the  final  paper  by  Brugger  [13],  explicit  equations  were  presented  for  these 
relations  for  a  number  of  crystal  symmetries,  wave  modes  and  stress  conditions.  These 
equations  allow  the  calculation  of  all  of  the  third  order  elastic  coefficients  for  all  of  the 
crystal  point  groups. 

Thus,  based  on  the  theory  of  Thurston  and  Brugger  [14],  measurements  need  only 
be  made  of  the  change  in  "natural"  velocity  (i.e.  changes  in  time  of  flight)  of  ultrasonic 
waves  as  a  function  of  stress  to  determine  the  third  order  elastic  coefficients.  However, 
the  changes  involved  are  generally  quite  small  and  require  very  accurate  ultrasonic 
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techniques.  Several  such  techniques  have  been  developed  over  the  years  which  make 
these  measurements  possible.  Early  measurements  by  Hughes  et  al.  [15]  in  1950  used  a 
simple  pulse  transmission  technique.  A  voltage  spike  was  used  to  excite  a  piezoelectric 
ultrasonic  transducer  which  launched  an  ultrasonic  wave  into  the  sample.  A  receiving 
transducer  on  the  opposite  side  of  the  sample  converted  the  elastic  wave  into  an  electri¬ 
cal  signal  which  was  displayed  on  an  oscilloscope.  The  time  of  flight  was  measured 
directly  on  the  scope  with  an  uncertainty  reported  to  be  -h/—  0.03  fis.  This  led  to  frac¬ 
tional  uncertainties  of  larger  than  one  part  in  10^  for  samples  with  a  20  ps  travel  time. 


Improvements  in  the  ability  to  measure  more  accurately  the  ultrasonic  wave 
time  of  flight  follov/ed  rapidly.  In  1953,  McSkimin  [16,17]  reported  on  a  phase  com¬ 
parison  ultrasonic  technique.  This  technique  was  a  pulse-echo  technique  which  used  the 
same  transducer  as  a  generator  and  receiver  of  the  elastic  waves.  A  tone-burst  {burst  of 
several  cycles  of  high  frequency  sinusoidal  v/aves)  was  used  to  excite  the  transducer  and 
thus  generate  the  ultrasonic  pulse.  The  elastic  pulse  traveled  through  the  sample  and 
reflected  off  of  the  other  side.  This  occurred  a  number  of  times  producing  an  oscillo¬ 
scope  pattern  of  a  number  of  decaying  amplitude  pulses  which  were  ail  similar  except 
for  the  exponential  decrease  in  their  amplitude.  In  McSkimin’s  technique  he  amplified 
and  rectified  the  received  echos  and  displayed  them  on  the  oscilloscope.  The  time  of 
flight  was  determined  by  varying  the  frequency  of  the  tone-burst  to  find  frequencies 
where  the  received  echos  were  all  in  phase.  For  this  condition  the  time  of  flight  was 


given  by 


360 

fn 


(1.10) 


where  t  was  the  time  of  flight, 
n  was  an  integer, 


f  was  the  frequency  of  the  tone-burst 
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and  7  was  a  phase  angle,  in  degrees,  correcting  for  phase  shifts  at  the 
specimen-transducer  interface. 

The  integer  n  was  determined  by  measuring  two  sequential  values  of  the  frequency  for 
the  "in  phase"  condition.  In  this  case,  n  was  given  by 


n 


Af 


(1.11) 


where  f^  and  were  the  two  measured  values  of  frequency 
and  Af  was  the  difference  between  the  two  measured  values. 

An  expression  for  the  phase  angle  7  based  on  standard  equivalent  electrical  transmission 
line  theory  was  presented.  It  was  calculated  from  known  material  properties  (mechani¬ 
cal  impedances  of  specimen  and  transducer  and  resonant  frequency  of  transducer).  The 
overall  uncertainty  in  measuring  velocity  was  reported  to  be  0.14  %. 

In  1957  Williams  and  Lamb  [18]  reported  a  technique  using  a  pair  of  ultrasonic 
bursts  and  a  through  transmission  technique.  The  first  burst  was  launched  through  the 
sample  and  reflected  multiple  times  through  the  specimen.  The  second  burst  was 
launched  through  the  same  transducer  a  short  time  later  so  that  the  first  arrival  of  the 
second  pulse  at  the  receiving  transducer  would  coincide  with  the  arrival  of  the  first  echo 
of  the  first  pulse.  The  excitation  frequency  was  then  adjusted  so  that  there  was  cancel¬ 
lation  between  the  two  signals.  Ajialysis  was  then  presented  which  allowed  calculation 
of  the  travel  time  in  the  specimen  from  the  frequency  which  included  corrections  for 
phase  shifts  at  the  specimen-transducer  interface.  Analysis  of  the  phase  angle  also 
included  the  effect  of  the  coupling  medium  between  the  transducer  and  specimen.  The 
measured  velocity  was  considered  to  be  accurate  to  one  part  in  10^  . 

Other  techniques  developed  later  such  as  the  pulse  superposition  method 
(McSkimin  [19])  in  1961  and  the  pulse  echo  overlap  technique  (May  [20]  and  Papadakis 
[21])  increased  absolute  accuracy  of  velocity  measurements  under  certain  conditions  to 
several  parts  per  million  according  to  Papadakis  [22] .  He  also  reported  that  changes  in 


time  of  flight  could  be  measured  accurate  to  parts  in  10^  using  these  techniques.  Details 
of  the  pulse  echo  overlap  technique  will  be  discussed  later. 


The  final  technique  discussed  was  reported  by  Hsyman  [23,24]  in  1980.  It 
allov/ed  changes  in  "natural”  velocity  accurate  to  one  part  in  10^  to  be  measured.  The 
technique  uses  an  ultrasonic  instrument  known  as  a  pulsed  phase  locked  loop  (P2L2) 
which  measures  changes  in  resonant  frequency  of  the  specimen,  transducer,  bond  compo» 
site  resonator  system  as  a  function  of  some  other  parameter  such  as  stress  or  tempera¬ 
ture.  This  can  be  shov/n  to  be  equivalent  to  measuring  changes  in  ''natural”  velocity 
(W)  as  follov;s.  Since 


VV 


t 


(1.12) 


v/here  Lq  is  the  initial  length  of  propagation  v/hich  is  a  constant 
and  t  is  the  time  of  flight, 


then 


AW 


(iALq  —  LoAt) 


(1.13) 


which  reduces  to 


AW  = 


LnAt 


(1.14) 


since 


ALo  =  0. 


Therefore, 


AW  _  _  At 
W  ~  t 


(1.15) 


(1.16) 


The  norma!  phase  velocity  (v)  is  given  by 
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where  L  is  the  length  of  propagation  which  is  not  a  constant. 
Therefore, 

AW  ^ 

W  “  t  “  V  L  ■ 


(1.18) 


However  the  resonant  frequency  of  the  test  system  is  given  by 


F 


n 


nv 

2L 


where  n  is  an  integer  denoting  the  number  of  the  harmonic. 
Therefore, 

ZlF  _  Av  zXL 

F  V  L  ^ 


(1.19) 


(1.20) 


and  thus, 


AF  _  AW  ^  _  At 
F  W  t 


(1.21) 


The  way  in  which  the  P2L2  measures  changes  in  resonant  frequency  will  be  discussed  in 
a  later  section  as  it  is  the  instrument  used  for  the  nonlinear  measurements  in  this  work. 

Although  most  of  the  work  discussed  thus  far  deals  with  the  development  of 
theory  and  measurement  techniques  for  nonlinear  behavior  of  materials,  much  work  in 
the  literature  points  out  the  importance  of  nonlinear  measurements  and  their  applica¬ 
tions.  One  of  the  first  applications  of  nonlinear  elasticity  was  in  the  area  of  materials 
characterization.  Nonlinear  elastic  constants,  through  their  representation  in  terms  of 
interatomic  potentials  in  crystalline  solids,  provide  much  information  about  the  nature 
of  the  bonding  of  atoms  in  a  crystal.  Einspruch  and  Manning  [25]  pointed  out  the 
importance  of  nonlinear  and  anharmonic  phenomena  in  solid  materials  in  relation  to 
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other  important  materia!  properties  such  as  thermal  coefficients  of  expansion,  Griineisen 
constants  and  ultrasonic  attenuation.  They  also  pointed  out  the  need  for  a  nonlinear 
description  of  elasticity  in  theories  dealing  v/ith  regions  of  high  stress  or  finite  strain  in 
solids  such  as  around  dislocations. 

One  of  the  most  useful  and  important  applications  of  nonlinear  elasticity,  how¬ 
ever,  has  been  in  the  field  of  nondestructive  evaluation  of  applied  and  residual  stress 
(strain)  in  materials.  This  application  is  based  on  the  fact  that  ultrasonic  wave  speeds 
are  a  function  of  the  stress  state  of  a  material  due  to  nonlinear  effects.  Thus,  if  the 
relationship  between  stress  and  velocity  is  previously  measured  for  a  given  material, 
theoretically,  the  state  of  stress  in  the  material  can  be  monitored  by  measuring  the 
ultrasonic  velocity.  A  variation  of  this  is  called  the  acoustoelastic  effect  or  acoustical 
birefringence.  The  acoustoelastic  effect  arises  from  the  fact  that  tv/o  shear  waves,  one 
polarized  parallel  to  stress  and  the  other  perpendicular  to  stress,  have  a  difference  in 
velocity  proportional  to  the  applied  stress.  Again  this  is  due  to  nonlinear  effects. 
Acoustical  birefringence  is  similar  to  optical  birefringence  or  the  photoelastic  effect 
which  has  been  used  to  determine  strain  in  transparent  materials  for  many  years.  As 
early  as  1959  the  acoustoelastic  effect  was  being  touted  for  the  measurement  of  residual 
stress.  Benson  and  Raeison  [26]  stated  that  acoustoelasticity  could  be  as  effective  as 
photoelastic  methods  with  the  additional  benefit  of  applications  to  opaque  materials. 

However,  in  the  years  since,  problems  have  plagued  the  application  of  acoustical 
birefringence  for  the  measurement  of  residual  stress.  The  effect  of  velocity  difference 
betv/een  the  shear  v/aves  due  to  material  anisotropy  is  often  much  larger  than  changes 
observed  due  to  stress  making  residual  stress  measurements  difficult,  .\nother  problem 
involves  energy  flux  deviation.  In  an  anisotropic  body,  the  energy  flux  vector  (the  direc¬ 
tion  of  the  fiov/  of  energy  per  unit  time  per  unit  area)  of  an  elastic  wave  does  not  in 
general  coincide  with  the  wave  normal.  This  refraction  of  the  wave  is  different  for  the 
two  shear  waves  which  means  the  tv/o  waves  do  not  follow  the  same  path  through  the 
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material.  It  also  depends  on  the  degree  of  anisotropy  and  may  also  change  as  a  func¬ 
tion  of  stress.  These  factors  make  velocity  measurements  for  acoustoelastic  determina¬ 
tion  of  stress  extremely  difficult.  Also,  as  pointed  out  by  Hsu  [27],  although  there  are 
analogies  between  the  photoelastic  effect  and  the  acoustoelastic  effect,  a  comparison  of 
the  important  parameters  (i.e.  wavelengths,  frequencies,  velocity  changes,  etc.)  for  both 
acoustical  and  optical  wave  propagation  shows  that  the  acoustic  technique  is  not  as 
efficient  as  the  optical  technique. 

Nevertheless,  numerous  investigators  are  still  attempting  to  overcome  these  prob¬ 
lems.  Theory  has  been  developed  to  attempt  to  account  for  the  effect  of  special  textures 
which  has  met  with  limited  success  in  applications.  AJso,  as  noted  earlier,  measurement 
technology  in  ultrasonics  has  progressed  substantially.  The  progress  in  acoustoelasticity 
was  reviewed  in  detail  by  Pao  et  al.  [28].  The  general  problem  of  the  experimental 
determination  of  residual  stress  using  ultrasonic  methods  has  yet  to  be  solved.  However, 
the  detection  of  applied  stress  has  been  successfully  reported  by  numerous  authors 
[24,29,30]. 

Nonlinear  elastic  material  properties  may  also  be  important  in  the  nondestruc¬ 
tive  evaluation  of  important  engineering  properties  of  materials.  For  engineering  appli¬ 
cations,  it  is  important  to  determine  material  properties  such  as  strength  and  residual 
strength  after  impact  or  fatigue  loading' and  thermal  cycling.  However,  most  nondes¬ 
tructive  measurements  yield  only  physical  properties  such  as  modulus,  density  or 
coefficient  of  thermal  expansion.  To  determine  the  important  engineering  properties 
relationships  must  be  developed  between  the  measured  physical  quantities  and  the 
desired  engineering  properties.  Nonlinear  properties  may  be  useful  in  this  respect  as 
they  are  parameters  which  are  significant  in  large  deformations  (i.e.  near  failure  strains) 
of  materials.  To  examine  this,  Heyman  et  al.  [31]  studied  the  relationship  between  the 
Stress  Acoustic  Constant  (SAC)  which  is  a  measure  of  a  mixture  of  second  and  third 
order  elastic  moduli  and  carbon  content  in  carbon  steels.  They  found  the  SAC  not  only 
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seemed  to  be  related  to  the  carbon  content  but  to  the  strength  of  the  material  as  well. 
In  aluminum,  Heyman  and  Chern  [32]  examined  the  relationship  between  the  SAC  and 
heat  treatment.  They  found  that  while  the  second  order  moduli  were  insensitive  to  the 
different  heat  treatments  used,  the  SAC  was  able  to  differentiate  them.  Thus,  a  simple 
test  method  for  detecting  improperly  heat  treated  and  thus  inferior  strength  aluminum 
was  demonstrated  based  on  measurements  of  nonlinear  elastic  properties. 

Although  previous  discussion  has  considered  primarily  only  metallic  and  non- 
metallic  crystalline  solids,  numerous  studies  have  presented  measurements  of  nonlinear 
phenom_ena  in  polymeric  and  composite  materials.  In  1950,  Hughes  et  al.  [15]  reported 
changes  in  wave  speed  as  a  function  of  pressure  and  temperature  for  polyethylene, 
polystyrene  and  Lucite.  In  1953,  Hughes  and  Kelly  [10]  presented  the  three  third  order 
elastic  constants  for  polystyrene.  In  1959,  Singh  and  Nolls  [33]  measured  the  velocity  as 
a  function  of  hydrostatic  pressure  and  temperature  for  polyisobutylene.  Assay  et  al. 
[34,35]  presented  the  change  in  velocity  for  polymethylm-ethacrylate  as  a  function  of 
stress  and  temperature.  The  relationship  between  stress  and  velocity  became  nonlinear 
at  high  pressures  indicating  higher  order  nonlinear  effects.  Lamberson  [35],  in  1969, 
reported  on  the  temperature  and  hydrostatic  pressure  dependence  of  velocity  in  polys¬ 
tyrene  and  two  composite  materials.  One  of  the  composites  was  carbon  phenolic  v/hile 
the  other  was  a  tape  wmund  silica  phenolic.  He  also  reported  nonlinear  stress-velocity/ 
curves  for  these  materials  at  high  pressures.  Zarembo  and  Shklovskaya  [37]  reported 
harmonic  generation  data  for  polystyrene,  plexiglass  and  rubber  in  1971.  They  also 
examined  the  effect  of  hydrostatic  and  uniaxial  stress  on  harmonic  generation  in  these 
materials.  Other  authors  [38,39]  have  reported  nonlinear  measurements  on  different 
polymeric  materials  since  this  time. 

The  present  work  deals  with  the  m.easurement  of  the  linear  and  nonlinear  elastic 
properties  of  a  graphite/epoxy  composite  material.  More  specifically,  measurements 
were  made  on  unidirectional  laminates  of  T300/5208  (Thornel  300  graphite  fibers  in  a 
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Narmco  5208  epoxy  resin).  The  choice  of  a  unidirectional  lay  up  was  made  because  it 
offers  the  highest  elastic  symmetry  and  therefore  is  the  least  complex  orientation  to 
analyze.  The  elastic  symmetry  of  unidirectional  fiber  reinforced  composites  is  usually 
taken  to  be  that  of  transverse  isotropy  if  the  fibers  are  randomly  distributed  in  the 
plane  perpendicular  to  the  fiber  axis.  A  transversely  isotropic  material  is  one  which  has 
a  preferred  axis,  perpendicular  to  which  is  a  plane  in  which  the  material  behaves  in  an 
elastically  isotropic  fashion.  In  the  case  of  a  unidirectional  composite,  the  fiber  axis  is 
the  preferred  axis  having  a  much  larger  elastic  stiffness,  and  the  isotropic  plane  is  the 
plane  perpendicular  to  the  fiber  axis  with  a  much  lower  stiffness. 

For  a  transversely  isotropic  material  there  are  five  linear  elastic  (second  order) 
moduli.  This  is  the  same  as  for  the  case  of  a  hexagonal  single  crystal.  However,  there 
are  only  nine  nonlinear  (third  order)  moduli  for  transverse  isotropy  as  opposed  to  either 
ten  or  twelve  for  the  different  point  group  symmetries  of  hexagonal  single  crystals.  In 
the  present  work,  the  five  linear  (second  order)  elastic  stiffness  moduli  for  the  unidirec¬ 
tional  composites  were  computed  from  ultrasonic  velocity  measurements.  These  values 
were  then  compared  with  available  data  from  the  literature  as  well  as  checked  for  con¬ 
sistency  with  the  assumption  of  transverse  isotropy.  From  the  linear  elastic  stiffness 
moduli,  the  compliance  moduli  were  calculated  and  compared  with  those  obtained  by 
strain  gauge  measurements.  These  were  also  compared  with  reported  values  for  a  simi¬ 
lar  material.  Measurements  were  then  made  to  determine  the  change  in  ’’natural"  wave 
velocity  as  a  function  of  applied  stress  and  temperature.  A  variety  of  combinations  of 
propagation  direction  and  loading  configurations  were  used  with  longitudinal  and 
transverse  waves.  This  data  was  also  checked  for  consistency  with  the  predictions  of 
transverse  isotropy.  Then  the  velocity-stress  data  was  used  to  calculate  some  of  the 
nonlinear  coefficients  of  this  material. 
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n.  LINEAR  ELASTICITY 


n^Ao  intirodiuiction 

In  order  to  determine  tiie  iionlinesr  eLstic  modulij  the  Imesr  elnstic  stiffness  End 
compliance  moduli  mmst  first  be  measured.  In  this  chapter,  measurements  of  these  pro¬ 
perties  in  a  unidirectional  graphite/epo:ry  composite  will  be  discussed.  Definitions  of  the 
terms  used  in  linear  elastic  theory  v/ill  be  presented  first,  followed  by  a  theoretical 
explanation  of  the  methods  used.  Then,  the  measurements  will  be  reported  and  com¬ 
pared  v/ith  theoretically  predicted  and  previously  measured  values  for  this  material. 


IIoBo  ThLeoiry 

As  discussed  earlier,  the  theoretical  basis  for  linear  elasticity  shows  that  it  is 
only  applicable  for  small  or  infinitesimal  deformations.  This  permits  the  use  of  simplify¬ 
ing  assumptions  in  the  definition  of  strain  for  linear  elasticity.  The  coordinates  of  a 
point  of  material  in  an  undeformed  body  are  defined  to  be  aj  with  respect  to  the  origin 
in  an  orthogonal  coordinate  system  (Note:  Einstein  convention  of  summation  over 
repeated  indices  is  assumed  throughout  this  paper).  The  same  point  of  material  defined 
by  aj  will  be  displaced  to  a  new  position  Xj  in  the  material  after  deformation.  The  dis¬ 
placement  is  a  vector  defined  by 

Tf=Y~a'’  (2.1) 


having  components  u^.  The  strain  can  now  be  defined  in  two  ways.  The  Lagrangian 


strain  tensor  which  provides  a  measure  of  the  deformation  with  respect  to  the  unde- 
formed  material  coordinates  is  defined  by 


du^  (9uj 

5aj  ^  5a,  ^  Saj^aj 


1 

2 


(2.2) 


In  the  Eulerian  formulation  of  strain,  the  strain  is  described  in  reference  to  the 


deformed  state  and  is  given  by 


2 


(9ui 

5xj  SajSaj 


(2.3) 


The  Lagrangian  formulation  is  often  referred  to  as  a  material  description  of  strain  while 
the  Eulerian  is  called  the  spatial  description.  The  Lagrangian  is  more  often  used  in 
solid  mechanics  and  the  Eulerian  in  fluid  mechanics.  The  assumption  of  infinitesimal 
deformation  can  now  be  applied.  This  assumption  means  that  since 


5ui 
(9a  j 


and 


«1 


(2.4) 


for  infinitesimal  deformations,  then  the  quadratic  term 


5ui5uk 


is  much  smaller 


(  ^aj^aj  J 

than  the  linear  terms  and  can  be  dropped  from  the  expression.  Thus  the  small  or 
infinitesimal  strain  tensor  is  defined  by 


1 

“  2 


du;  du, 
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(2.5) 


The  small  strain  assumption  also  means  that 
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<9ui  5u, 
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da,  da, 


2. 
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dx. 
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dx 


J ) 


(2.6) 


or  that  there  is  no  difference  between  the  Lagrangian  and  Eulerian  formulations  for 
strain.  It  can  also  be  seen  from  the  definition  of  the  small  strain  tensor  that  it  is  a  sym¬ 
metric  tensor  (i.e.  =  Cjj  ). 

The  physical  description  of  the  components  £22,  and  €33  is  that  they  are  the 
normal  strains  which  provide  the  measure  of  the  normalized  change  in  length  of  the 
body  along  the  respective  coordinate  axis.  For  example,  if  Iq  is  the  unstrained  length 
along  axis  Xj  and  is  the  strained  length  along  Xj  ,  then 
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(h  ~  ^o) 


(2.7) 


The  remaining  strain  components  are  the  shear  strains  v/hich  are  a  measure  of  the 
change  in  angle  between  tv/o  originally  orthogonal  axes  in  the  undeformed  medium. 
This  is  demonstrated  in  Fig.  2.1  where  the  strain  is  given  by 

^12  “ -  ~  tan(0)  ^  9  (for  small  strain).  (2.8) 

X2 


The  stress  tensor  ((7y)  provides  a  m_easure  of  the  applied  forces  on  the  body.  It  is 
defined  by 

_  force  on  the  i^th  face  in  the  j^th  direction 
cross  sectional  area  of  the  i'th  face  * 

A  diagram  of  the  nine  components  of  the  stress  tensor  is  shov/n  in  Fig.  2.2.  The 
engineering  or  nominal  stress  is  used  where  the  cross  sectional  area  is  that  measured  in 
the  unstressed  state.  The  true  stress  uses  the  instantaneous  cross  sectional  area  which 
is  much  m.ore  difficult  to  measure.  The  components  (7^,  (722  and  are  the  normal 
stress  components  while  the  remainder  are  shear  stress  components.  Laws  of  statics  can 
be  used  to  show  that  if  a  body  is  to  be  in  rotational  equilibrium,  the  stress  tensor  must 
also  be  symmetric. 

Having  defined  both  stress  and  strain,  the  relationship  between  the  two  can  now 
be  discussed.  For  linear  elasticity,  the  generalized  Hooke’s  law  is  used.  This  expresses 
the  relationship  between  stress  and  strain  as 

^\]  ^ijkl^kl  (2.10) 

where  Cyi^j  is  the  fourth  rank  linear  elastic  stiffness  tensor.  In  the  most  general  form  this 
tensor  has  eighty-one  components.  Flowever,  because  the  stress  and  strain  tensors  are 
both  symmetric,  the  following  relations  are  valid 
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^ijkl  ^jikl 

and 

^ijkl  ”  ^ijlk*  (2.12) 

This  reduces  the  number  of  independent  components  to  thirty-six.  The  assumed 
existence  of  a  strain  energy  density  function  (<;6)  which  for  the  case  of  linear  elasticity 
has  the  form 

(j)  =  (2-13) 

provides  further  reduction  of  the  number  of  independent  elastic  coeflBcients.  It  must  be 
true  that 

Equating  this  expression  term  by  term  with  respect  to  strain  yields  the  relation 

^ijkl  ^  ^klij-  (2-15) 

This  reduces  the  number  of  independent  coefficients  to  twenty-one.  At  this  point,  the 
Voigt  notation  can  be  used  to  simplify  things  somewhat.  For  each  pair  of  indices  on 
stress,  strain  and  elastic  moduli,  the  following  substitutions  are  made 

11  -►  1  23-^4 

22  — ►  2  31  5  (2.16) 

33  — ►  3  12  — ►  6 

Using  this  notation,  the  generalized  Hooke’s  law  can  be  rewritten  as 


where  capital  subscripts  are  summed  from  one  to  six.  The  elastic  stiffness  moduli  can 
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now  be  written  in  matrix  form  for  the  most  general  anisotropic  solid  as 


^12 

Ci4 

Cl5 

^16 

^12 

^^22 

^23 

^24 

^25 

*^26 

^13 

<^23 

C33 

C34 

C35 

^36 

^14 

^24 

C34 

C44 

•^45 

^46 

^15 

^25 

C35 

^45 

^55 

^56 

^16 

^26 

C36 

*^46 

^56 

^66 

The  reduction  in  the  number  of  independent  elastic  constants  can  be  carried 
further  if  material  symmetries  are  taken  into  consideration.  This  is  accomplished  by 
examining  the  effect  on  the  strain  energy  density  function  of  a  rotation  of  the  material 
which  places  the  material  in  an  elastically  equivalent  condition  due  to  symmetries  of  the 
material.  Due  to  the  scalar  nature  of  the  strain  energy  density  function,  it  should  have 
the  same  value  in  both  configurations.  An  example  of  this  procedure  is  now  given.  The 
strain  energy  density  function  is  written  out  in  terms  of  the  twenty-one  independent 
coefficients  as 


<P 


■(^11^1  V  ^22^1  -r  Cggef  -i-  C44e|  4-  Cggef  4-  cggcl)  4-  c 


'12^1^2 


V  CiseiCa  4-  4-  Cjgejeg  Cigeife  V  ^23^2^3  +  ^24^2^4  ^25^2^5 

V  ^26^2^6  V  ^34^3^4  V  -l-  csgegCg  4-  +  C4ge4eg  4-  cggegCg  (2.19) 


In  an  isotropic  material  any  rotation  is  elastically  equivalent.  Assuming  an  isotropic 
material,  a  rotation  of  180  degrees  around  the  X3  axis  is  chosen  and  is  shown  in  Fig.  2.3. 
The  transformation  matrix  or  direction  cosine  matrix  [ajj]  is  found  by  taking  the  cosine 
of  the  angle  between  the  rotated  and  original  axes.  For  this  rotation  it  is  given  by 

r  T 


-1 

0 

0 


0  0 
-1  0 
0  1 


(2.20) 
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The  elect  of  this  rotation  on  the  strain  tensor  is  given  by 


Thus 


^ik^jl^kh 


(2.22) 


(2.21) 


The  rotated  strain  energy  density  function  can  be  written  after  substituting  equation 
(2.22)  as 


4  = 


(  2 
■{^11^1 


^22^2  +  <^33^3  “!■  ^44^1  “S”  ^55^5  ^66^1)  "i"  ^^12^1^2 


-1-  0136^63  —  0146264  —  C156265  0256265  +  O236263  —  C246264  —  O256265 

+  <^26^2^6  ~  C346364  —  O356365  -h  0356365  -j-  0456465  —  0456465  —  0556565).  (2.23) 


Evaluating  term  by  term  the  relation 


0'  =  0 


(2.24) 


yields 

0i4  =  —  C24 
^15  ~  ~  ^15 
^24  ~  “  ^24 
^25  ~  ~  ‘^25 

O34  =  —  C34  (2.25) 

C35  ~  —  C35 
*“46  ^46 

^56  ~  ~  <^56 

all  other  Cjj  =  Cjj. 


If  equations  (2.25)  are  to  be  true,  then 
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^14  =  <^15  =  <^24  ~  ^25  —  ^34  “  <^35  “  ^46  —  ^55  —  0.  (2.26) 

Application  of  further  symmetry  operations  will  show  that  for  an  isotropic  material,  the 
number  of  independent  elastic  coefficients  reduces  to  two.  The  matrix  then  appears  as 
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Ci2 

Ci2 
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0 

0  ~ 

Ci2 

Cll 

Ci2 

0 

0 

0 

(2.27) 
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Cll 

0 
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0 
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^11  ~  ^12 

0 

0 

t 

2 

0 

0 

0 

0 

^11  —  ^12 

0 

2 

0 

0 

0 

0 

0 

<^11  “  ^12 

2 

For  a  transversely  isotropic  material,  the  number  of  independent  coefficients  is  five  while 
for  an  orthotropic  material  it  is  nine.  The  derivation  of  the  independent  coefficients  as 
well  as  their  final  matrix  forms  for  these  two  elastic  symmetries  which  are  important  for 
unidirectional  fiber  reinforced  composites  is  shown  in  Appendix  A. 

The  relationship  between  strain  and  stress  is  given  by 

=  SijkiOki.  (2.28) 

The  moduli  Syi^i  form  the  linear  elastic  compliance  tensor.  The  conditions  of  symmetry 
for  the  stiffness  tensor  also  apply  to  the  compliance  tensor.  The  relations  between  the 
compliance  and  stiffness  moduli  for  a  given  symmetry  can  be  determined  by  inverting 
the  stiffness  matrix. 

The  value  of  the  elastic  moduli  are  also  dependent  upon  the  conditions  (isother¬ 
mal  or  adiabatic)  in  which  they  are  measured.  Isothermal  moduli  are  measured  at  con¬ 
stant  temperature  (T)  while  adiabatic  moduli  are  measured  at  constant  entropy  (S). 
The  work  by  Brugger  [12]  provides  a  precise  thermodynamic  definition  of  the  elastic 
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moduli  and  the  relationship  between  adiabatic  and  isothermal  coefficients.  In  terms  of 
the  free  energy  (F),  the  isothermal  stiffness  moduli  are  defined  by 


^ijkl  ^^0 


jr 


(2.29) 


where  is  the  unstrained  density. 
The  adiabatic  stiffness  moduli  are 

5 

^ijkl  ”  Po 


(2.30) 


where  U  is  the  internal  energ^n 

The  adiabatic  and  isothermal  compliances  are  defined  in  terms  of  the  enthalpy  (H)  and 
the  Gibbs  function  (G)  respectively  and  the  thermodynamic  tensions  (ty)  which  ai 

au 
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The  compliance  moduli 
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(2.32) 


and 


T  [  (9^G 

Sijkl  =  -  Po 


It 


(2.33) 


The  relationship  between  adiabatic  and  isothermal  compliance  moduli  was  given  in 
matrix  form  by  Hankey  and  Schule  [40]  as 


_„S  ^ 

Sab  —  Sab 


T  1 

lcoCpG“= 


(2.34) 


where  Cp  is  the  specific  heat  at  constant  pressure, 
and  a.^  is  the  linear  coefficient  of  thermal  expansion  defined  by 
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^Va 

O^A  =  (2-35) 

The  difference  between  adiabatic  and  isothermal  moduli  is  usually  quite  small  and  often 
less  than  experimental  uncertainty. 

The  determination  of  the  adiabatic  linear  elastic  stiffness  moduli  can  be  made 
by  measuring  the  density  and  the  ultrasonic  wave  velocity.  The  derivation  of  the  equa¬ 
tions  of  motion  for  the  propagation  of  linear  elastic  waves  in  an  anisotropic  medium 
provides  the  theoretical  basis  for  this  technique.  This  subject  has  been  treated  by  a 
number  of  authors  among  the  earliest  being  Love  [41]  and  it  is  reviewed  in  detail  by 
Green  [5].  In  deriving  the  wave  equations,  a  number  of  initial  assumptions  about  the 
propagation  medium  are  made.  It  is  assumed  to  be  unbounded,  homogeneous  and  con¬ 
tinuous.  Next,  an  infinitesimal  volume  element  of  the  medium  is  examined  and  the  net 
unbalanced  forces  acting  on  it  are  determined.  This  is  accomplished  by  finding  the 
resultant  of  the  variation  of  stress  across  the  element  for  each  direction  and  multiplying 
this  by  the  respective  cross  sectional  area.  Fig.  2.4  shows  the  variation  of  stress  along 
the  Xi  direction  for  the  infinitesimal  volume  element  with  sides  of  length  <5xi,  <5x2 
If  there  were  no  variation  of  stress  and  thus  no  net  unbalanced  forces,  there  would  be 
no  wave  motion.  The  net  unbalanced  force  along  x^  is  given  by 


[(^11  +  +  [(^21  +  ^^^2)  “ 

1  2 


dCTr 


-  -  3]^ 

+  [(<^31  +  ■^'^3)  -  <73i]&i&2- 


(2.36) 


This  can  be  reduced  to 


5cr, 


<9xi 


11  dCTc 
(9xo 
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5x2  <5X3. 


(2.37) 


Similar  expression  can  be  written  for  the  net  unbalanced  force  along  X2  and  X3  .  Accord¬ 
ing  to  Newton’s  second  law  these  unbalanced  forces  must  equal  miii  where  m  is  the  mass 
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of  the  element  and  \i[  are  the  second  time  derivatives  of  the  displacements  which  are  the 
acceleration  components.  Actually,  to  be  more  exact,  it  is  pointed  out  that  Newton  s 
second  law  states  that  the  force  is  equal  to  the  time  derivative  of  the  momentum  pj  . 
However, 

Pi  =  mUj  (2.38) 

and  therefore 

Pi  =  mill  +  mU|.  (2.39) 


Since  the  mass  of  the  element  is  assumed  to  be  constant  in  time  the  forces  are  simply 

given  by  mUj  .  Therefore  the  balance  of  forces  can  be  written  in  component  form  as 

/  \ 


mUi  = 


5Xi 


(2.40) 


It  is  noted  that  the  body  forces  are  neglected  in  this  equation.  The  mass  can  be  rewrit¬ 
ten  in  terms  of  the  density  Po  and  the  volume  (V)  which  is  given  by 

V  =  &Ci5x2&3  (2-41) 


as 


m  =  /3o<5xi&2^3- 


(2.42) 


Substituting  (2.42)  into  (2.40)  and  canceling  appropriate  terms  yields 

Poiii  = 


aCTy 


(9xi 


(2.43) 


To  rewrite  this  in  terms  of  only  density  and  displacement,  the  stress  is  first  rewritten  in 
terms  of  the  strain  using  Hooke’s  law  (2.10).  Then  the  strain  definition  (2.6)  is  used  to 
produce  the  expression  for  stress  in  terms  of  displacement  as 


(2.44) 
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<9xi  ^  5xk 


The  symmetry  of  the  elastic  coefficients  (2.12)  is  applied  to  further  simplify  (2.44)  to 
yield 


(J-,.  =  — e, 

ij  2  ‘ 


<7;; 


,  1 
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bjlk 
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[  ^ 

<9x,  1 

(2.45) 


(2.46) 


This  is  now  put  bade  into  equation  (2.43)  to  yield  the  final  form  of  the  equation  of 
motion  for  linear  elastic  wave  propagation  in  an  anisotropic  medium  which  is 


5x|(9xj 


(2.47) 


Now  that  the  equation  of  motion  has  been  derived,  the  next  step  is  to  choose  a 
solution  for  the  displacements  as  a  function  of  time  and  space.  For  a  plane  wave,  a 
solution  can  be  written  as 

u,(xj„t)  =  Aoa,e‘^‘"^  ~  (2.48) 

v/here  is  the  amplitude  of  the  v/ave, 

Ofj  are  the  direction  cosines  of  the  particle  displacement  vector, 
t  is  the  time, 

a;  is  the  angular  frequency 
and  kn,  are  the  components  of  the  v/ave  vector. 

The  wave  vector  components  (k^^)  are  related  to  the  direction  cosines  of  the  wave  nor¬ 
ma!  and  the  wavelength  (X)  by  the  expression 
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where  k  is  the  wave  number. 

The  solution  for  the  displacement  can  now  be  substituted  into  the  equation  of  motion. 
This  will  lead  to  the  following  expression 


ijkik^lilj^k  = 


which  can  be  rewritten  as 


where  v  is  the  phase  velocity  of  the  elastic  wave  and  is  given  by 


In  equation  (2.51),  it  is  noted  that  is  arbitrary  and  therefore  not  necessarily  equal  to 
zero.  In  order  for  this  equation  to  have  nontrivial  solutions,  it  must  be  true  that  the 
determinant  of  the  matrix  of  coefficients  must  equal  zero.  That  is 


CijkiMi  -  =  0- 


To  more  conveniently  write  this  equation,  a  matrix  is  defined  by 


\k  ^ijklMj* 


Therefore,  equation  (2.53)  becomes 


I  ^ik  P6^  ^ik  I 


This  is  a  characteristic  equation  which  when  expanded  forms  a  cubic  expression  in  terms 
of  PqV^  .  Since  is  a  symmetric  matrix,  the  three  solutions  for  PqV^  must  all  be  real. 
The  physical  significance  of  this  is  that  for  any  arbitrary  direction  in  an  anisotropic 
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medium,  three  plane  waves  may  be  propagated.  The  velocity  of  each  of  these  waves 
depends  on  the  density,  the  elastic  constants  and  the  direction  cosines  of  the  wave  nor¬ 
mal.  Thus,  if  the  ultrasonic  wave  velocities  are  measured  for  a  number  of  directions  in 
an  anisotropic  medium  with  known  density  and  elastic  symmetry,  the  linear  elastic 
stiffness  moduli  can  be  determined.  Although  not  formally  designated  as  such  in  the 
preceding  equations,  these  moduli  will  be  the  adiabatic  stiffnesses. 

The  particle  displacement  vector  direction  cosines  can  be  determined  for  each  of 
the  three  v/aves  propagating  along  a  given  direction.  This  is  accomplished  by  solving 
for  the  eigenvectors  for  each  eigenvalue  [pov"^)  using  the  equation 

PqV  Qj  (2.56) 

and  also  applying  the  fact  that 

=  1.  (2.57) 

These  direction  cosines  determine  the  mode  of  the  v/ave.  If  they  are  the  same  as  the 
wave  normal  direction  cosines,  that  is 


(2.58) 

then  the  particle  displacement  of  the  elastic  wave  lies  along  the  direction  of  propaga¬ 
tion  of  the  wave.  A  wave  of  this  type  is  referred  to  as  a  pure  mode  longitudinal 
(comipressional)  wave.  If  the  particle  displacement  is  perpendicular  to  the  propagation 
direction  and  therefore 


«il,  =  0,  (2.59) 

then  the  wave  is  called  a  pure  mode  transverse  (shear)  v^ave.  Illustrations  of  longitudi- 
nal  and  transverse  mode  v/aves  are  given  in  Figs.  2.5  a)  and  b). 


Figure  2 
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5  -  (a)  Illustration  of  a  longitudinal  (compressional) 

elastic  wave  demonstrating  directions  of  particle 
vibration  and  propagation. 
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Figure  2,5  -  (b)  Illustration  of  a  transverse  (shear)  wave 
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However,  in  general  the  wave  need  not  be  pure  mode.  It  can  be  either  quasilong- 
itudinal  or  qnasitransverse  depending  on  which  type  of  particle  motion  is  more  dom¬ 
inant.  Of  the  three  waves  propagating  along  a  given  direction,  one  is  longitudinal  or 
quasilongitudinal  while  the  other  two  are  transverse  or  quasitransverse. 

/mother  consideration  of  linear  elastic  wave  propagation  in  an  anisotropic 
medium  is  that  the  direction  of  the  energy  flux  of  the  v/ave  may  not  in  general  be  the 
same  as  the  direction  of  the  wave  front  normal.  That  is,  there  may  be  a  refraction  of 
the  wave  because  of  the  anisotropy  of  the  medium.  The  energy  flux  vector  is  defined  as 
the  direction  of  the  flow  of  energy  per  unit  time  per  unit  area.  A  diagram  of  the  energy 
flux  deviation  of  a  quasitransverse  wave  is  shov/n  in  Fig.  2.5.  Equations  for  expressing 
the  energy  flux  vector  and  its  deviation  from  the  direction  of  the  wave  normal  are 
reviewed  by  Green  [5]  and  are  presented  here. 

First,  the  expressions  for  the  kinetic  energy  (K)  and  the  potential  energy  (P) 
which  is  contained  in  the  wave  field  are  given  by 


(2.60) 


and 


P 


-  Jcr|j€ijd^ 


(2.61) 


where  dr  is  a  volume  elemient  of  the  v/ave  field. 

The  sum  of  K  and  P  yields  the  total  elastic  energy  E  contained  in  the  wave  field  which 
is 


1  (  o  1 

=  =  7  [Jpufdr  -h  /cTije,jdr  . 


(2.62) 


The  time  derivative  of  H  is  then  given  by 


=  /pu.iiidr  + 


dt 


(2.63) 
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To  simplify  this  expression  the  term  under  the  second  integral  is  now  considered.  It  can 
be  rewritten  as 


(2.64) 


Hooke’s  law  (2.10)  can  now  be  substituted  into  this  equation  to  express  — ((Tyfy)  in 

(J% 

terms  of  the  elastic  moduli  and  the  strains  as 


A 

di 


(2.65) 


Because  of  the  symmetry  of  the  elastic  moduli,  equation  (2.85)  can  be  rewritten  as 


(2.86) 


Again  using  Hooke’s  lav/  this  becomes 


A 

dt 


(^ij^ij)  —  ^kl^kl  “i"  ^11^ 


ij  1] 


(2.67) 


which  can  be  rewritten  as 


d  ,  \  c,  ■ 


(2.68) 


after  the  dummy  indices  k  and  1  are  summed  out  of  the  equation.  The  time  derivative 
of  the  strain  in  equation  (2.68)  can  be  rewritten  in  terms  of  the  displacements  using  the 
definition  of  strain  as 


^ij 


1 
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(^lij 

dx. 


(2.69) 


This  is  substituted  into  equation  (2.68)  resulting  in 


A 

dt 


3xj 


(2.70) 
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which  because  of  the  symmetry  of  the  stress  tensor  is 


d 


(o-jje.j)  =  2a-,j 


5Ui 


dt  ^  ax. 


(2.71) 


Equation  (2.71)  is  now  put  back  into  equation  (2.63)  to  give  the  time  derivative  of  the 
total  energy  of  the  elastic  wave  field  as 


as  f  •  ^  f  , 

—  =  /pu.Uidr  +  J  ^ij— 


(9xi 


(2.72) 


Since 


(9 

+ “i- 


aui  aoii 


axj  ’ 


(2.73) 


or 


(2.74) 


the  second  integral  in  equation  (2.72)  can  be  rewritten  as 


,  ailj  ,  S  .  aCT;; 


(2.75) 


Gauss’s  theorem  can  now  be  applied  to  the  first  integral  of  equation  (2.75)  to  rewrite  it 
as  a  surface  integral  of  the  form 

/^(<^,jUi)d7-  =  /(Tjp.dS.  (2.76) 


This  can  be  substituted  into  equation  (2.75)  which  is  then  inserted  into  equation  (2.72) 
which  gives 
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puj  - 


dxi 


Ujdr  +  J  (7jjUidS. 


(2.77) 
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The  first  integral  equals  zero  because  of  the  equation  of  motion  (2.47).  Therefore,  it 
remains  that 


/o-ijU,dS. 


(2.78) 


This  equation  represents  the  energy  flov/  into  the  volume  element  r  through  the  surface 
S.  The  energy  flux  vector  v/ith  components  Ej  can  be  formed  which  gives  the  energy  flux 
flowing  outward  from  the  volume  r  through  S  by 

Ej  =  -  (2.79) 


If  the  time  averaged  values  E  of  the  energy  flux  vector  over  one  cycle  of  the  elastic 
v/ave  are  computed  and  it  is  defined  that 

E|  =  Ef  -f  E|  -K  E|,  (2.80) 


then  the  direction  cosines  of  the  energy  flux  vector  (E()  can  be  formed  by 
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Tl 

Et  ' 


(2.81) 


If 

E,'  =  1„  (2.82) 

then  the  v/ave  propagates  along  the  direction  of  the  wave  normal  and  there  is  no  energy 
flux  deviation.  Otherv/ise,  the  angle  (<5)  of  deviation  of  the  energy  flux  from  the  wave 
normial  can  be  determined  by 

8  =  COs“^(Eilj).  (2.83) 


Therefore,  to  determine  the  elastic  moduli  by  using  ultrasonic  velocity  measure¬ 
ments,  a  number  of  considerations  must  be  examined.  First,  the  elastic  symmetry  and 


thus  the  form  of  the  elastic  stiffness  matrix  must  be  known.  Then  a  direction  of  propa¬ 
gation  must  be  chosen  and  the  solutions  for  the  velocities  of  the  three  waves  that  pro¬ 
pagate  along  this  direction  must  be  calculated  in  terms  of  the  density  and  the  elastic 
moduli.  Also,  any  energy  flux  deviation  of  the  three  wave  modes  must  be  determined. 
This  is  repeated  for  a  number  of  different  directions  until  enough  independent  equations 
are  derived  to  calculate  all  of  the  elastic  moduli  from  the  corresponding  velocity  meas¬ 
urements.  Usually,  an  attempt  is  usually  made  to  choose  directions  which  propagate 
waves  which  are  pure  mode  and  suffer  no  energy  flux  deviations  to  make  the  measure¬ 
ments  as  easy  as  possible.  Unfortunately,  this  can  not  always  be  accomplished. 
Brugger  [42]  published  an  equation  which  can  be  used  to  determine  the  pure  mode  direc¬ 
tions  for  the  various  elastic  symmetries.  In  general  form  this  equation  is 


^ijk^jlrs^kMr^s  ® 


(2.84) 


where  is  the  permutation  tensor. 

It  is  defined  by  ^ 

0  if  any  i,j,k  are  equal 
1  if  ij,k  are  in  cyclic  order 
1—1  if  i,j,k  are  in  acyclic  order 


^ijk 


(2.85) 


Also  presented  are  the  pure  mode  wave  directions  for  the  various  elastic  symmetries. 

The  specific  case  of  elastic  wave  propagation  in  a  unidirectionally  graphite  fiber 
reinforced  composite  will  now  be  discussed.  First,  the  applicability  of  the  previously 
derived  wave  equations  to  a  composite  material  must  be  considered  with  respect  to  the 
heterogeneous  nature  of  the  medium.  This  violates  one  of  the  initial  assumptions  that 
the  medium  be  homogeneous.  However,  it  is  generally  assumed  that  if  the  wavelength 
of  the  elastic  waves  is  much  larger  than  the  size  of  the  inhomogeneity,  then  the  material 
is  macroscopically  homogeneous  and  the  wave  equations  are  applicable.  In  the  materi¬ 
als  studied,  the  size  of  the  inhomogeneity  is  given  by  the  diameter  of  the  fibers  which  is 
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on  the  order  of  seven  fim  .  In  the  frequency  range  used  2.25IV/IHz)  ,  the  wavelengths 
for  the  different  v/ave  modes  and  directions  of  propagation  ranged  from  about  700/im  to 
almost  5000/im  .  Therefore,  the  previously  stated  assumptions  are  valid  as  are  the  wave 
equations. 

Another  important  question  is  the  use  of  Hooke’s  law  because  of  the  viscoelastic 
nature  of  most  polymeric  materials.  However,  epoxies  and  their  composites,  because  of 
the  high  degree  of  crosslinking  along  the  molecular  chains,  generally  exhibit  little  viscoe¬ 
lastic  behavior.  This  is  especially  true  at  the  temperatures  at  which  these  measure¬ 
ments  were  made  (room*  temperature).  Measurements  of  the  wave  speeds  and  thus  the 
moduli  as  a  function  of  frequency  demonstrated  no  viscoelastic  effects.  Therefore  the 
application  of  Hooke’s  law  appeared  to  be  valid. 

The  appropriate  elastic  symmetry  model  for  a  unidirectionally  fiber  reinforced 
graphite/epoxy  composite  must  now  be  chosen  in  order  to  derive  the  specific  solutions 
for  the  velocities  of  ultrasonic  wave  propagation.  The  most  common  symmetry  model 
used  for  unidirectional  composities  is  that  of  transverse  isotropy.  In  a  transversely  iso¬ 
tropic  miaterial,  there  exists  a  preferred  axis  perpendicular  to  which  there  exists  a  plane 
in  which  the  material  behaves  isotropically.  In  unidirectional  graphite/epoxy,  the  fiber 
direction  (X3)  is  the  preferred  axis  while  the  X1X2  plane  is  the  plane  of  isotropy  as  shown 
in  Fig.  2.7.  Transversely  isotropic  behavior  of  unidirectional  composites  is  based  on  the 
assumption  of  a  random  distribution  of  fibers  in  the  X1X2  plane  as  well  as  alignment  of 
the  fibers  along  the  X3  axis.  The  symmietry  conditions  for  such  a  material  are  that  it  be 
elastically  insensitive  to  twofold  (180  degree)  rotations  about  each  of  the  three  axes  as 
Vvel!  as  insensitive  to  any  rotation  around  the  X3  axis.  Application  of  these  symmetry 
conditions  as  shov/n  in  Appendix  A  demonstrates  that  there  are  only  five  independent 
elastic  moduli.  The  matrix  form  of  the  elastic  stiffness  moduli  for  a  transversely  isotro 
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However,  if  for  some  reason  the  fibers  are  not  randomly  spaced  in  the  X1X2  plane, 
this  plane  will  no  longer  be  elastically  isotropic.  In  this  case,  the  material  will  possess 
an  orthotropic  symmetry.  This  might  occur  if  because  of  the  lay  up  process  of  manufac¬ 
ture  there  v/as  a  resin  rich  layer  between  each  lamina  as  shown  in  Fig.  2.8  or  if  some  of 
the  laminae  were  misoriented.  The  resulting  orthotropic  material  would  now  have  only 
the  tv/o  fold  symmetries  about  each  of  the  three  axes  as  is  the  case  with  cross  ply  or 
angular  ply  composite  lay  ups.  For  this  symmetry  the  number  of  independent  linear 
elastic  moduli  increases  to  nine  as  shown  in  Appendix  A.  The  matrix  form  of  the  moduli 
is  then 
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Ultrasonic  velocity  data  can  be  compared  with  the  equations  derived  for  both 
models  to  determine  which  gives  the  best  fit  to  the  data.  The  derivations  of  the  equa- 
tions  for  the  velocity  as  a  function  of  density  and  elastic  moduli  for  a  number  of 
different  directions  and  pure  wave  m.odes  for  transversely  isotropic  and  orthotropic  sym¬ 
metries  are  presented  in  Appendix  B.  However  a  summary  of  the  results  is  presented  in 
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Figure  2. 


resin 
ri  ch 
areas 


-  Illustration  of  unidirectional  composite  with  resin 
rich  regions  between  laminae  causing  the  material  to 
have  orthotropic  elastic  symmetry* 
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Table  2.1.  To  determine  all  of  the  elastic  moduli  for  the  two  models,  several  nonpure 
mode  waves  must  be  used.  The  results  for  these  waves  are  presented  in  Table  2.2. 


UloC, 

In  most  of  the  ultrasonic  velocity  measurements  made  in  this  study,  a  pulse  echo 
overlap  velocity  measurement  system  was  used.  This  system  was  briefly  discussed  in  the 
Introduction.  Hov/ever,  the  details  of  its  operation  v/ill  nov/  be  discussed.  Figure  2.9 
shows  a  block  diagram  of  the  entire  pulse  echo  overlap  system.  A  Matec  model  6800 
pulse  modulator  and  receiver  v/as  used  to  generate  an  electrical  tone  burst  of  variable 
amplitude  and  frequency.  This  v/as  used  to  excite  an  ultrasonic  transducer  (typically 
undamped  crystals  of  quartz  or  PZT  with  fundamental  resonances  of  approximately 
2.25  IvIHz.)  which  excited  ultrasonic  waves  into  the  graphite/epoxy  samples.  The  result¬ 
ing  echoes  from  the  back  surface  of  the  specimen  were  detected  by  the  same  transducer 
and  converted  back  into  electrical  signals.  The  received  signals  were  input  into  the 
receiver  and  amplified.  After  amplification,  they  were  displayed  on  a  Hewlitt  Packard 
(H.P.)  1743A  dual  trace  oscilloscope.  A  typical  pulse  echo  pattern  is  shown  in  Fig.  2.10. 
The  echoes  appear  identical  to  the  generating  signal  (main  bang)  except  for  an  exponen¬ 
tial  decrease  in  amplitude  due  to  attenuation  in  the  sample. 

To  have  the  pulse  echo  pattern  appear  stable  on  the  oscilloscope,  the  Matec  gen¬ 
erator  and  the  PI.P.  i743A  oscilloscope  must  be  triggered  simultaneously.  This  was 
accomplished  by  using  the  triggering  components  of  the  pulse  echo  overlap  system.  The 
signal  from  a  continuous  wave  (c.w.)  oscillator  (Matec  model  110)  was  input  into  a 
Matec  122B  decade  divider  and  dual  delay  generator.  The  Matec  122B  divided  the  ori¬ 
ginal  signal  by  a  preset  amount  (typically  1000)  and  then  generated  a  trigger  signal  at 
the  divided  frequency.  This  is  demonstrated  for  frequency  division  by  five  in  Fig.  2.11. 
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Table  2.1  -  Equations  for  pure  mode  linear  elastic  wave  propagation  in 
transversely  isotropic  and  orthotropic  media. 


Propagation  Particle  Displacement  Wave  Transversely  Orthotropic 

Direction  Direction  Mode  Isotropic  Model 

Cosines  Cosines  Model 
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Table  2.2  =  Equations  for  non  pure  mode,  off  axis  linear  elastic  wave 
propagation  in  transversely  istotropic  and  orthotropic  media. 
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Table  2.2  -  (Continued) 


Output  of  Matec  122  B 


Figure  2.1 


1  -  Illustration  of  frequency  division  of  a  continuous 
wave  (c.w.)  source  by  the  Matec  model  122B  by  a 
factor  of  five. 


Figure  2ol2  -  Illustration  of  an  oscilloscope  trace  of  two  echoes 
during  the  pulse  echo  overlap  technique  where  the 
echoes  are  not  properly  overlappedo 


Figure  2ol3  -  Illustration  of  two  properly  overlapped  echoes = 
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uncertainty  in  the  operator's  visual  ability  to  detect  the  overlapped  condition  and  the 
uncertainty  in  the  measurement  of  the  c.w.  oscillator  frequency. 

Because  of  time  delays  of  the  ultrasonic  pulse  in  the  bond  and  phase  shifts  at  the 
specimen-transducer  interface,  corrections  had  to  be  made  to  determine  the  true  time  of 
flight  (<5)  from  the  measured  time  (t).  The  theory  for  these  corrections  was  presented  by 
McSkimmin  [19]  and  is  outlined  in  Appendix  C.  In  this  theory  it  is  necessary  to  deter¬ 
mine  the  frequency  of  the  ultrasonic  pulse.  This  was  accomplished  by  inputting  the  tone 
burst  into  a  mixer  (Anzac  MD-40)  along  with  the  signal  from  a  HP.  3325A  frequency 
synthesizer/function  generator.  The  output  of  the  mixer  was  then  input  into  a  Tek¬ 
tronix  53C  oscilloscope.  The  output  of  the  mixer  was  the  sum  and  difference  frequencies 
of  the  tone  burst  and  the  HP.  c.w.  source.  If  they  were  of  the  same  frequency,  the  out¬ 
put  of  the  mixer  was  a  constant  voltage  (0  Hz)  dependent  on  the  phase  difference 
between  the  two  signals.  This  phase  difference  was  variable  each  time  the  Matec  6600 
was  triggered  and  thus  a  series  of  parallel  lines  was  generated  on  the  scope  trace.  This 
is  shown  in  Fig.  2.14.  If  the  frequencies  were  different,  then  the  output  would  have  the 
difference  frequency  again  phase  shifted  on  each  new  trace.  Therefore,  the  frequency  of 
the  tone  burst  from  the  Matec  6600  could  be  monitored  by  the  frequency  of  the  H.P.  fre¬ 
quency  synthesizer. 

The  remaining  equipment  was  used  to  collect  data.  A  Tektronix  4051  computer 
was  used  to  read  the  trigger  frequency  values  and  calculate  the  time  of  flight  which  was 
used  to  calculate  the  velocities  and  the  moduli.  The  program  used  to  control  the  pulse 
echo  overlap  data  acquisition  and  analysis  is  shown  in  Appendix  D. 

There  were  a  number  of  factors  which  affected  the  size  and  geometry  of  the  sam¬ 
ples  used.  They  had  to  have  parallel  sides  to  be  used  in  the  pulse  echo  technique.  Also, 
the  thicker  the  specimen,  the  more  accurately  its  time  of  flight  could  be  measured. 
However,  because  of  the  high  ultrasonic  attenuation  of  graphite/epoxy  composites,  espe¬ 
cially  for  shear  waves,  there  were  limits  on  the  maximum  thickness.  Also,  limits  on 


sample  thickness  occurred  because  of  the  desire  to  have  the  wavefront  remain  planar. 
The  distance  (d)  over  which  the  wavefront  remains  planar  was  given  by  Mason  [43]  as 


where  R  is  the  radius  of  the  transducer 
and  X  is  the  wavelength  of  the  ultrasonic  wave. 

However  limitations  on  the  thickness  because  of  the  large  attenuation  usually  occurred 
before  this  requirement  became  important.  There  were  also  restrictions  on  the  lateral 
dimensions  of  the  samples.  The  equation  of  motion  are  based  on  the  assumption  of  a 
wave  propagating  through  a  unbounded  medium.  However,  ultrasonic  velocities  have 
been  shown  by  Tu  et  al.  [44]  to  be  independent  of  lateral  dimensions  if  the  sample  was 
at  least  two  and  one  half  times  the  acoustic  wavelength. 

The  above  constraints  provided  the  basis  for  the  choice  of  the  sample  geometries 
used.  Ideally,  the  measurements  of  the  elastic  moduli  would  be  made  using  only  one 
sample  to  avoid  any  problems  with  sample  to  sample  variations  in  the  elastic  properties. 
These  problems  are  especially  prevalent  in  polymer  composite  materials  where  elastic 
property  variations  may  occur  not  only  from  batch  to  batch  due  to  different  curing  pro¬ 
cedures  but  also  within  a  given  batch  because  of  regional  variation  of  fiber  content,  void 
content  and  other  important  parameters.  However,  because  of  the  above  constraints 
and  the  necessity  to  use  off  axis  non  pure  mode  wave  propagation  to  determine  some  of 
the  moduli,  several  samples  were  required.  These  samples  were  all  cut  from  a  150  ply 
laminate  of  T300/5208  unidirectional  graphite/epoxy  that  was  nominally  12  X  15  X  0.9 
in.  (30.5  X  38.1  X  2.3  cm.)  after  cure.  This  panel  was  designated  MCIS  1  and  layed  up 
and  cured  at  NASA,  Langley  Research  Center.  The  finished  panel  was  ultrasonically 
C-scanned  by  the  Materials  Characterization  and  Instrumentation  Section  at  NASA 
Langley  for  the  presence  of  gross  voids  or  delaminations.  None  were  found  to  be 
present,  A  number  of  tensile  and  compression  samples  were  then  cut  and  machined 
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from  this  panel.  Those  used  for  the  ultrasonic  testing  were  MCIS  l.lOA,  1.8A,  I.IOB, 
1.7A,  1.5  A,B,E,  and  F  and  l.oA^.  ¥/hose  dimensions  are  listed  in  Tables  2.3  and  2.4.  The 
samples  were  generally  flat  and  parallel  to  -h/-  O.OCOS  in.  {O.0D076  cm.)  or  better.  Sam¬ 
ples  MCIS  l.lOAj  i.8A,  l.lOB  and  1.7A-  ¥/ere  all  rectangular  parallelepipeds  >vith  their 
faces  oriented  perpendicular  to  the  three  coordinate  axes.  The  densities  for  each  are 
also  given  in  Table  2.3  v/hich  show  some  sample  to  sample  variation.  The  remaining 
samples  v/ere  cut  to  allow  wave  propagation  in  a  direction  that  v/as  45  degrees  betv/een 
two  of  the  axes  and  perpendicular  to  the  third.  These  samples  were  used  to  propagate 
the  waves  listed  in  Table  2.2.  The  sample  dimensions  are  given  in  Table  2.4.  The  densi¬ 
ties  of  these  samples  were  not  measured  but  the  value  of  1.54  g/cmT  was  used  for  all  cal¬ 
culations  involving  these  specimen. 

Measurem_ents  of  the  pure  miode  longitudinal  wave  speeds  along  each  of  the  three 
coordinate  axes  were  made  using  specimen  MCIS  I.IOA.  Specimen  MCIS  1.8A,  l.lOB 
and  1.7A  were  used  to  determine  the  six  diflerent  combinations  of  polarisation  and  pro¬ 
pagation  directions  for  the  pure  mode  shear  waves  along  the  three  axes.  The  velocities 
and  the  product  PqY^  for  each  of  these  waves  are  shown  in  Table  2.5.  The  uncertainties 
for  the  shear  wave  velocities  were  larger  than  those  of  the  longitudinal  velocities 
because  of  the  much  higher  attenuation  of  shear  waves.  Even  using  the  thinner  speci¬ 
men  for  the  shear  waves  did  not  overcoKLe  the  high  attenuation.  This  resulted  in  a 
much  smaller  amplitude  second  echo  which  was  difficult  to  reproducibly  overlap.  The 
uncertainty  in  the  longitudinal  velocity  along  Xg  was  also  large  because  of  the  high  velo¬ 
city  and  thus  the  short  time  of  flight  along  this  direction. 

The  measured  values  were  then  compared  with  theory  to  check  for  transversely 
isotropic  behavior.  The  products  PqY^  calculated  using  the  longitudinal  waves  along 
Xi  and  X2  in  sample  MCIS  l.lOA  should  be  equal  as  should  the  products  calculated  from 
the  two  shear  waves  propagating  along  X3  in  sample  MCIS  1.8A.  Comparisons  for  the 
other  similar  products  were  not  valid  because  they  ’were  obtained  on  different  samples. 
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Table  2.3  -  Sample  dimensions  and  densities  for  samples  used  for  pure 
mode  elastic  wave  propagation. 


Sample 

Number 

Dimensions  (in.) 

Xi  X2  X3 

Density 

(g/cm^) 

I.IOA 

0.8001 

0.7999 

0.8002 

1.5404  +/-  0.0004 

l.lOB 

0.7999 

0.5002 

0.7996 

1.5460  +/-  0.0006 

1.8A 

0.8002 

0.8000 

0.5002 

1.5384  +/-  0.0005 

1.7A 

0.5001 

0.8001 

0.7999 

1.5411  +/-  0.0005 

Table  2.4  -  Sample  dimensions  for  samples  used  for  off  axis,  non  pure 
mode  wave  propagation. 


Sample 

Number 

Dimension  (in.)  along  axis 

45  degrees  between  axes 

Dimension  (in.)  along  axis 

1.5A 

0.4499 

Xl,  X3 

0.9 

X2 

1.5B 

0.5000 

Xl,  X3 

0.9 

X2 

1.5E 

0.4491 

X21  X3 

0.9 

Xl 

1.5F 

0.4993 

X2,  X3 

0.9 

Xl 

1.6A 

0.5001  ' 

Xl,  X2 

0.9 

X3 

Table  2.5  -  Data  for  pure  mode  linear  elastic  wave  propagation. 


Sample 

Propagation 

Polarization 

Velocity 

Number 

Direction 

Direction 

(lO^cm/s) 

(GPa) 

l.lOA 

Xl 

Xl 

3.0467  +/-  0.0005 

14.295  +/-  0.005 

l.lOA 

X2 

X2 

3.0390  +/-  0.0008 

14.226  +/-  0.006 

l.lOA 

X3 

X3 

8.390  +/-  0.007 

108.4  +/-  0.1 

1.7A 

Xl 

X2 

1.557  +/-  0.002 

3.736  +/-  0.007 

1.7A 

Xl 

X3 

1.849  +/-  0.003 

5.27  +/-  0.01 

l.lOB 

X2 

Xl 

1.558  +/-  0.003 

3.75  +/-  0.01 

l.lOB 

X2 

X3 

1.850  +/-  0.003 

5.29  +/-  0.01 

1.8A 

X3 

Xl 

1.850  +/-  0.002 

5.265  +/-  0.008 

1.8A 

X3 

X2 

1.851  +/-  0.002 

5.271  +/-  0.008 
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multiple  results  were  averaged  except  for  C13  and  C23  where  there  was  a  wide  range  of 
values.  This  wide  range  was  due  to  the  large  uncertainty  in  the  measurements  which 
propagates  into  an  even  larger  uncertainty  in  the  calculated  values.  The  values 
presented  for  C13  and  C23  were  calculated  using  the  quasishear  data  because  the  quasi¬ 
longitudinal  data  yielded  complex  values  for  the  moduli. 

The  values  for  the  transversely  isotropic  model  can  be  compared  with  those 
reported  by  Kriz  and  Stinchcomb  [45].  They  reported  data  for  the  moduli  of  T300/5208 
graphite/epoxy  measured  ultrasonicaliy  as  well  as  the  values  theoretically  predicted 
from  the  properties  of  the  fibers  and  matrix.  The  theoretical  values  are  derived  for  an 
assumed  volume  fraction  of  fibers  of  0.67.  These  are  all  presented  in  Table  2.8.  The 
values  all  compare  favorably  except  for  C33  which  was  much  lower  in  the  present  study. 
The  lower  density  of  the  samples  used  in  the  present  experiments  implies  that  there  may 
be  a  lower  fiber  content  and/or  a  higher  void  content.  This  would  account  for  some  of 
the  discrepancy  in  the  value  of  C33  .  There  may  also  be  some  misalignment  of  the  fibers 
along  the  X3  axis  which  would  also  lower  the  value  of  C33  although  there  is  no  other 
confirming  evidence  of  this  fact.  The  values  of  Kriz  and  Stinchcomb  also  show  much 
more  scatter  demonstrating  sample  to  sample  variations. 


To  calculate  the  nonlinear  elastic  moduli^  the  isothermal  linear  elastic  compli¬ 
ance  moduli  are  also  needed.  The  adiabatic  compliance  moduli  can  be  obtained  by 
inverting  the  stiffness  matrix.  Corrections  for  differences  in  adiabatic  versus  isothermal 
moduli  can  then  be  applied  if  they  are  larger  than  the  uncertainties  of  the  inverted 
values.  For  the  case  of  transverse  isotropy  the  compliances  are  given  by  the  following 
equations, 

2 

^ _ *^11^33  ~  ^13 _ 

^33(Cn  ~  (^12)  +  20^3(0^2  —  C,,) 


S33  — 


^11  +  C12 

<^33(^11  +  C12)  —  20^3 


(2.92) 
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Table  2.8  -  Comparison  of  measured  transversely  isotropic  stiffness 
data  with  that  presented  by  Kriz  and  Stinchcomb  [45]. 
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Table  2.9  -  Calculated  linear  elastic  compliance  values  using  ultrasonic 
stiffness  data  for  transversely  isotropic  and  orthotropic  symmetries. 
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and 


~(CllC2:3  —  C13C12) 

S23  =  - ^ ^ ^ ^ - .  (2,104) 

^11 1^22^33  —  C23)  —  Ci2(Cr2C33  ““  2C13C23)  —  CI3C22 

The  calculated  values  for  these  moduli  are  also  presented  in  Table  2.9.  Again  there  was 
a  narrow  range  of  values  except  in  the  constants  S13  and  S23  because  of  the  large  uncer¬ 
tainty  in  Ci3  and  Cos  . 

The  isothermal  correction  for  these  adiabatic  compliance  moduli  can  be  calcu¬ 
lated  using  equation  (2.34).  The  temperature  of  these  measurements  was  nominally  24 
C  and  the  density  was  1.54  g/cm^  .  The  coefficients  of  thermal  expansion  for  similar 
unidirectional  carbon/epoxy  composites  were  given  by  Hull  [46]  as 
'  10“^/C  and  Q'3  =  —  0.2  ■  10“^/C  and  the  specific  heat  was  given  by  as 
Cp  =  2.1  J/gC  (unpublished  data  from  C.  Welsh,  1986).  Using  these  values,  the  max¬ 
imum  correction  term  was  only  about  7  *  10“®  GPa“^  which  was  much  smaller  than  the 
uncertainty  of  the  calculated  values.  Thus  the  difference  between  the  adiabatic  and  the 
isothermal  values  for  the  compliance  moduli  was  negligible. 

Some  of  the  elastic  compliance  moduli  w’ere  also  measured  using  static  compres¬ 
sive  loading  vvdth  bonded  strain  gauges  to  measure  the  strain.  The  theory  for  this  is 
based  on  Hooke’s  lanv.  If  a  unidirectional  stress  is  applied,  for  example,  along  the  Xj 
axis,  the  strain  can  be  written  as 


61  —  Sji  cri , 

(2.105) 

€0  ”  S  2  ; 

(2.106) 

63  =  Si.3(7i 

(2.107) 

and 


£4  —  t's  —  eg  =  0. 


(2.108) 


Therefore  the  moduli  can  be  rewritten  as 


SiA  — 


d<Ji 


or  in  more  general  form  as 


Sab  = 


dox’ 


(2.109) 


(2.110) 


for  a  unidirectional  stress  along  the  A  axis.  Therefore  the  stress  and  strain  can  be  meas¬ 
ured  and  the  slope  of  the  stress-strain  curve  will  provide  the  moduli. 

The  samples  used  in  these  measurements  were  cubes  nominally  0.7  in.  (1.8  cm.) 
on  a  side  oriented  with  faces  perpendicular  to  the  coordinate  axes.  They  were  also  cut 
from  the  panel  MCIS  1  and  were  designated  MCIS  1.8D  and  E.  The  load  was  applied 
with  a  MTS  810  load  frame  in  compression  and  the  calibrated  output  of  the  load  cell 
was  input  into  a  H.P.  3478A  digital  multimeter.  To  insure  that  the  load  remained  uni¬ 
directional,  a  fixture  was  designed  and  inserted  into  the  load  frame  that  placed  the 
upper  compression  plate  on  a  ball  bearing.  This  allowed  the  compression  plate  to  com¬ 
pensate  for  any  slight  lack  of  parallelism  in  the  specimen.  Fig  2.15  illustrates  this 
fixture.  The  strains  were  measured  with  Measurements  Group  Inc.  350n  strain  gauges. 
The  surfaces  of  the  samples  were  cleaned  according  to  manufacturer’s  instructions  and 
the  strain  gauges  were  then  bonded  onto  the  samples.  A  conventional  Wheatstone 
Bridge  circuit  was  used  with  a  H.P.  6002A  power  supply  providing  a  five  volt  input  and 
a  H.P.  3478A  voltmeter  measuring  the  output  voltage.  A  Tektronix  4051  computer  was 
used  to  read  the  data  from  the  voltmeters  and  calculate  stress,  strain  and  the  resulting 
modulus.  The  control  program  is  shown  in  Appendix  E.  Fig  2.16  shows  the  block 
diagram  for  the  experimental  set  up.  During  each  test,  a  unidirectional  stress  was 


applied  and  the  corresponding  strain  was  measured  in  two  directions.  One  direction 
was  parallel  to  the  direction  of  applied  stress  while  the  other  was  perpendicular.  Then 


Sail  Bsaring 


)p%r  Compression 
Plate 


Specimen 


Plats  ' 


Figure  2^15  -  Diagram  of  the  compression  fixture  used  during  the 
static  loading  for  the  measurement  of  the  elastic 
compliance  moduli  and  the  uniaxial  stress  acoustic 
constants  ^ 
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Figure  2.16  -  Block  diagram  of  the  apparatus  used  for  static 
measurement  of  the  linear  elastic  compliance 
moduli. 
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the  sample  was  reoriented  to  apply  stress  in  a  different  direction. 

The  resulting  values  for  the  compliance  moduli  are  given  in  Table  2.10.  The  shear 
moduli  Sd4,  S55  and  s^g  were  not  measured  by  this  technique.  Again  the  deviation  from 
transversely  isotropic  behavior  is  cpiite  small  as  demonstrated  by  the  good  agreement 
between  s^  and  S02  and  betv/een  S13  and  S23  for  the  orthotropic  moduli.  These  values 
compare  favorably  with  the  values  obtained  from  the  inversion  of  the  ultrasonic  data 
except  for  the  value  of  S12  •  This  discrepancy  could  be  caused  by  error  accumulation  in 
the  inversion  of  the  stiffness  moduli  of  it  could  be  due  to  inaccuracy  of  the  static  meas¬ 
urement.  Ledbetter  and  Read  [47]  also  compared  static  measurements  to  ultrasonic 
measurements  for  a  NbTi  liber  reinforced  Cu  composite  and  reported  discrepancies  in 
the  off  diagonal  moduli.  The  measured  values  can  also  be  compared  to  those  presented 
by  Garber  [48]  which  are  listed  in  Table  2.11.  His  measurements  were  made  in  tension 
on  thin  laminates  of  T300/5208.  Because  thin  laminates  w'ere  used,  Si2  was  not  meas¬ 
ured.  The  values  compare  favorably  except  for  S33  .  This  is  consistent  with  the  lower 
value  of  C33  reported  earlier  for  the  material  in  the  present  study. 


Il.Dc  Conclusion 

The  linear  elastic  stiffness  and  compliance  moduli  of  a  unidirectional 
graphite/epoxy  composite  have  been  measured  by  ultrasonic  and  static  measurements. 
These  measurements  demonstrate  a  measurable  deviation  from  transversely  isotropic 
behavior.  However,  the  deviation  is  very  small,  especially  wTen  compared  to  material 
variations  from  sample  to  sample,  and  thus  the  assumption  of  transverse  isotropy  for 
the  linear  elastic  behavior  of  this  material  is  valid. 

Comparisons  of  the  present  data  wdth  that  presented  in  the  literature  for  similar 
materials  are  favorable  except  for  the  values  of  C33  and  S33  .  The  difference  in  these 
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Table  2.10  -  Measured  linear  elastic  compliance  values  for  transversely 
isotropic  and  ortliotropic  symmetries. 


Transversely 

Isotropic 

Modulus 

Value 

(GPa)-* 

Orthotropic 

Modulus 

Value 

(GPa)-i 

Sll 

0.089 

Sll 

0.0890 

S12 

-0.063 

S12 

-0.0626 

Sl3 

-0.0021 

Sl3 

-0.00208 

S33 

0.00935 

S33 

0.00935 

S22 

0.0891 

S23 

-0.00215 

Table  2.11  -  Comparison  of  measured  compliance  values  with  those 
presented  by  Garber  [49]. 


Modulus 

Present 

Garber 

— 

Work 

(GPa)-^ 

Sll 

0.089 

0.0921 

S12 

-0.063 

not  measured 

Sl3 

-0.0021 

-0.0025 

S33 


0.00935 


0.0078 
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moduli  together  with  a  lower  measured  density  may  indicate  a  lower  fiber  volume  con¬ 
tent  or  a  higher  void  content  in  the  material  of  the  present  study.  It  also  may  indicate 
either  a  poor  alignment  of  the  fibers  or  a  waviness  of  the  fibers  along  the  fiber  axis. 

The  comparisons  of  the  ultrasonic  and  static  data  are  good  except  for  the 
modulus  Si2  .  The  difference  may  be  an  artifact  of  the  measurement  or  the  calculations 
used  in  the  comparisons.  Because  of  the  higher  accuracy  of  ultrasonic  measurements, 
the  values  used  in  the  computation  of  the  nonlinear  moduli  in  the  next  chapter  are 
those  obtained  by  ultrasonic  measurements  except  in  the  case  of  S13  where  large  uncer¬ 
tainties  exist. 
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m.  NONLINEAR  ELASTICTIY 


in.A*  Introduction 

Using  the  linear  elastic  stiffness  and  compliance  moduli  presented  in  the  previous 
chapter,  the  nonlinear  elastic  moduli  can  be  calculated  from  measurements  of  the  stress 
dependence  of  the  ultrasonic  **naturar  velocity.  In  this  chapter,  the  theory  of  how  this 
is  accomplished  will  be  presented  along  with  actual  measurements  and  calculations  for 
the  unidirectional  graphite/epoxy  composites  being  studied. 


in.B.  Theory 

The  theory  for  nonlinear  elasticity  differs  in  several  respects  from  the  theory  of 
linear  elasticity  which  was  presented  earlier.  The  first  major  difference  is  that  the 
assumption  of  infinitesimal  deformation  is  no  longer  valid.  Therefore,  the  strain  tensor 
must  be  defined  using  either  the  Lagrangian  or  Eulerian  definitions.  In  this  work  the 
Lagrangian  formulation  will  be  used.  The  constitutive  equation  (i.e.  relation  between 
stress  and  strain)  must  be  extended  to  include  higher  order  terms.  Also  the  strain 
energy  function  must  be  extended.  It  is  given  by  Brugger  [12]  as 

^  +  '  ’  '  •  (3.1) 

The  numerical  coefficient  for  each  of  the  terms  is  —  where  n  is  the  order  of  strain  for 

n! 

that  term.  The  Voigt  reduced  notation  can  be  used  to  rewrite  this  equation  as 

^  +  S  ^AB^A^  +  "rE^AAA^A 

^  A  A<B  A 

+  T  E  CaAB^A%  +  E  CABD^A^^d  +  •  •  •  (3-2) 

^  Af^B  A<B<D 
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¥/hers  no  sum  is  implied  except  where  indicated  and  the  capital  indices  range  from  1  to 

6. 

The  definitions  for  the  elastic  moduli  are  similar  to  those  of  linear  elastic  theor^n 
The  higher  order  adiabatic  and  isothermal  moduli  were  defined  by  Brugger  [12]  as 

(3.3) 


(3.4) 


„s 

—  Po 

,9^U 

^ijklmn... 

•  •  • 

=  Poj 

5“F 

^ijklmn... 

.S 

<9"H 

®ijklmn...  Po 


(3.5) 


and 


^ijklmn. 


=  -  Po 


<9"G 


(3.6) 


However,  in  this  research,  the  third  order  moduli  were  obtained  by  measuring  the  stress 
derivative  (isothermal)  of  the  ultrasonic  velocity  (adiabatic).  Therefore  the  resulting 
moduli  were  the  "mixed”  third  order  moduli  as  defined  by  Brugger  [12]  as 


.M 


dc 


ijkl 


jiT 


(3.7) 


or  in  reduced  notation 


.M 

^ABD 


(9c  s 


AB 


<9Pd 


(3.8) 


They  are  related  to  the  pure  adiabatic  third  order  moduli  by 

i^ijklqr^qr  ” 


^ijklmn  "  ^ijklmn 


T  s 

mnop^op 


Po^T 


dc 


ijkl 


dT 


fJ 


(3.9) 


or 
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r  ( 'I  ■ 

CABD  =  ^ABD  +  CabF^  “  ^ 

However,  again  the  differences  are  usually  small  and  therefore  less  than  experimental 
uncertainty. 

Using  the  same  arguments  based  on  symmetry  of  the  stress  and  strain  tensors 
and  the  existence  of  the  strain  energy  function,  the  number  of  independent  third  order 
elastic  moduli  can  be  reduced  to  56  for  the  most  general  anisotropic  material.  Further 
considerations  of  the  symmetry  of  the  material  reduces  the  number  to  20  for  an  ortho¬ 
tropic  material  and  9  for  a  transversely  isotropic  material  as  shown  in  Appendix  F .  For 
the  orthotropic  case  the  independent  moduli  are  Cm,  Cii2>  Ci22; 
^123?  ^133>  ^144>  ^166>  ^222)  ^223>  ^233?  ^244?  ^266>  ^266)  ^333>  ^344?  ^366?  ^368  ^460  •  Other 

moduli  are  equal  to  zero.  The  independent  moduli  for  transverse  isotropy  are 
Cm,  c„2,  CU3,  Ci23,  Ci33,  C144,  c,66,  C333  and  C344  with  the  following  relations  for  the  other 
moduli 

^222  =  ^111,  ^355  ==  ^344 

^223  “  ^113,  ^166  ”  ^266  “  “^(^111  ""  ^112) 

^255  ~  ^144,  ^233  ^  ^133  (3.11) 

^244  “  ^155,  ^366  ”  ""  ^^23) 

^122  ==  ^112,  ^456  =  “^(^155  ^144) 

with  all  remaining  c^bd  ==  ^  • 

As  mentioned  previously,  these  moduli  can  be  obtained  from  measurements  of  the 
stress  derivative  of  the  ultrasonic  "natural"  velocity.  The  theoretical  basis  for  this  was 
presented  by  Thurston  and  Brugger  [14].  They  rederived  the  equation  of  motion  for  a 
small  amplitude  elastic  wave  propagating  in  a  homogeneously  stressed  medium  using  the 
following  notation.  The  coordinates  aj  refer  to  the  position  of  a  point  of  material  in  the 
"natural"  or  unstressed  state.  Xj  are  coordinates  of  the  homogeneously  stressed  material 
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=  N,N3(5jk^rs  +  (V  +  2^qk)Cj?qs)-  (3-18) 

The  eigenvector  U  is  normalized  so  that  its  components  are  the  direction  cosines  of  the 
particle  displacement  with  respect  to  the  natural  state.  Differentiating  equation  (3.17) 
with  respect  to  pressure  (p)  and  evaluating  at  zero  pressure  yields 

(PoW“);  .  0  =  (U.Cc^'kUJp  .  „  (3.19) 

which  was  reduced  to  the  following  forms  for  the  conditions  of  hydrostatic  pressure  and 
uniaxial  compressive  stress.  For  hydrostatic  loading, 

— (Po^^)p  =  0  =  1+  2wFhc  +  GhC  (3-20) 

where  Fhc  =  sJaraUfUs  .  (3-21) 

Ghc  =  sTuvC„vprqaNpN,UA  ,  (3.22) 

and  w  =  (/>oW2)o  =  (/>v2)o  =  CpS  ,,NpN,UA  •  (3-23) 

For  uniaxial  loading  the  equation  reduced  to 

-(PoW2)p  =  0=2wF^e  +  Guc-  (3-24) 

where  and  G^5  are  given  by 

F..  =  sJ,„M.MpUA  (3-25) 

and 

^uc  “  S^^yC^ypJ.qsM3^MbNpNqUpU5  (3.26) 

and  M|  are  the  components  of  a  unit  vector  in  the  direction  of  stress. 

It  is  assumed  for  the  uniaxial  case  that  the  direction  of  propagation  is  always  perpen¬ 
dicular  to  the  direction  of  stress. 
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Tliese  equations  provide  relationships  between  the  linear  elastic  stiffness  and 
compliance  moduli,  the  third  order  elastic  moduli  and  the  derivative  of  the  "natural” 
velocitv  with  respect  to  stress.  However,  the  measurements  made  in  this  study  were  of 
the  Suress  Acoustic  Constants  (H)  which  were  delined  by  Heyman  [31]  and  Cantrell  [49] 
to  be 


H  = 


(3.27) 


Since  there  \vere  a  number  of  different  orientations  of  stress  and  propagation  and  polarH 
zation  directions,  the  different  SAC’s  were  designated  with  three  subscripts  as 
The  first  subscript  gave  the  direction  of  applied  stress  which  was  limited  to  one  of  the 
three  axes  for  these  experiments.  A  zero  for  the  first  subscript  implies  hydrostatic  load- 
ing.  The  second  and  third  subscripts  w-ere  the  directions  of  propagation  and  polariza¬ 
tion  of  the  ultrasonic  wave  respectively.  These  again  were  limited  to  being  pure  mode 
waves  along  the  three  axes. 

Using  this  notation,  the  equations  for  the  SAC’s  were  derived  for  both  the  models 
of  transverse  isotropy  and  orthotropy.  The  results  for  the  27  cases  considered  are  given 
in  table  3.1. 


nirC,  Experiment 

The  nieasui  e ments  of  the  stress  derivati\''es  of  the  natural"  velocitv  "were  earned 
out  using  some  oi  the  samples  that  were  used  in  the  previously  described  ultrasonic  velo¬ 
city  measurements.  Sample  MCIS  l.lOA  was  used  for  the  measurements  of  the  longitu¬ 
dinal  velocity  derivatives  while  samples  MCIS  1.8A,  l.lOB  and  1.7A  were  used  in  the 


shear  wave  measurements. 


-  81  - 


Table  3.1  -  Stress  Acoustic  Constant  Equations  for  transversely  isotro¬ 
pic  and  orthotropic  media. 


Transverse  Isotropy 


Hn„  =  - 


1  +  2cn(sii  +  Si2  +  Sis)  +  (sii  +  Si2  -f  Si3)(cni  +  C112)  +  (2Si3  +  833)01 


Hnoo  —  — 


H022  —  Hqh 

1  +  2c33(2Si3  H-  S33)  +  2{sii  +  Si2  +  813)0 133  +  (2Si3  +  833)0333 


Hnqi  —  — 


1  -f  2C44(Sii  H-  S12  +  S13)  +  (s^i  +  S12  +  Si3)(ci44  -f  C155)  +  (2Si3  +  833)0344 


Hnio  — 


1  -h(Cn  —  +  ^12  ^13)  +  ^((^11  +  S12  S]^3)(Cm  —  0^12) 

(cii  —  C12) 

_  (^^13  +  S33)(<^113  ^123) 

(Cii  “  C12) 


H021  —  H012 


=  - 


1^211  H122 

2CiiSi3  +  813(0111  +  O112)  +  S33O113 
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Table  3ol  -  Continued 


Orthotropy 

Hon  =  - 

1  4  2ciiSi  4  SjCiii  4  S2C112  +  S3C113 

2cii 

H022  = 

1  4  2C22S2  4  S1C122  +  S2C222  4  S3C223 

2c22 

Ho33  =  — 

1  4  2C33S3  4  S1C133  4  §2^233  +  §3^333 

2C33 

II 

1 

1  +  2C55S1  4  S1C155  4  S2C256  4  S3C355 

2^66 

Ho32  —  “ 

1  4  2C44S2  4  S1C144  4  S2C244  4  S3C344 

2C55 

^012  ”  “ 

1  4  2c0qS2  4  SiCjQQ  4  S2C268  4  S3C306 

2C68 

1 

II 

§ 

1  4  2ce8Si  4  SjCj88  4  S2C2C8  4  §30380 

2c00 

Hoi3  —  — 

1  4  2C55S3  4  S1C155  4  §2^256  +  §3^356 

2^56 

^023  =  “ 

1  4  2C44S3  4  S1C144  4  §2^244  +  §3^366 

2C44 

II 

1 

2C33S13  4  SnCi33  4  S12C233  4  Si3C333 

2C33 

H233  —  ” 

2C33S03  4  S12C133  4  S22C233  “t  ^23^333 

2^33 

H122  ^  ~ 

2^22^12  4  S11C122  4  S12C222  d-  S13C223 

2^22 

11 

1 

2ciiSj2  4  S12C111  4  S22Cn2  +  S23C113 

2cii 

where  Sj  =  +  Sjg  +  Sjs 

^2  ~  ®12  ^22  "h  ^23 

and  S3  =  Si3  4-  S23  +  S33 


Continued 


Orthotropy 


2^22^23  S13C122  So3^222  “J"  ^33^2 


2CiiSi3  H-  S13C211  -f  S23C212  H-  §33^113  1 


"^66^11  't-  SjjCigQ  -h  S12C266  "t-  Si3C3g0 


^^66^22  ^12^166  "t  §22^236  ^23^366 


SceeSos  ■!"  SisCjgg  H-  S23C236 


‘^^66^13  ^13^166  S23C263  ^1“  S33C3gg 


3C44S13  +  S22C244  -h  S22C244  4“  S23C344 


-^55^23  +  S22C253  +  S22C255  -h  S23C3^ 


"^44^33  4-  S23C244  "h  S23C244  -h  S33C344 


^^55^33  4-  S23C153  +  S23C255  -h  S33C355 


^^55^11  +  S22C255  4-  S12C255  4-  S23C353 

2CeE: 


2C44S22  4“  S20C244  -h  S22C244  4-  S23C344 


2C44S22  4-  S2iC244  4-  S22C244  4“  S23C344 


[  -^55^12  +  S22C2.55  4-  S22C206  4-  S23C355 


ifi' 


II 
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The  change  in  "natural"  velocity  was  measured  using  a  pulsed  phase  locked  loop 
(P2L2)  interferometer.  The  theory  of  the  measurement  of  the  change  in  resonant  fre¬ 
quency  and  it’s  relationship  to  the  change  in  "natural"  velocity  was  presented  earlier. 
Now  the  theory  of  the  operation  of  the  P2L2  which  is  shown  in  Fig,  3.1  will  be  dis¬ 
cussed. 

The  heart  of  the  P2L2  is  a  voltage  controlled  oscillator  (VCO)  which  generates  a 
continuous  wave  signal  of  a  frequency  that  is  controlled  by  a  D.C.  input  signal.  A  por¬ 
tion  of  the  signal  from  the  VCO  is  gated  out  into  a  tone  burst  and  is  then  used  to  excite 
the  ultrasonic  transducer.  The  resulting  elastic  wave  is  launched  into  the  sample  and  is 
then  detected  either  by  the  same  transducer  in  the  pulse  echo  mode  or  another  trans¬ 
ducer  in  the  through  transmission  mode.  The  received  signal  is  input  into  the  P2L2 
where  it  is  phase  compared  with  the  signal  from  the  VCO  at  a  preselected  phase  point 
using  a  sample  and  hold.  The  sampled  voltage  from  the  phase  detector  is  then  used  to 
drive  the  VCO  to  a  condition  of  quadrature.  The  acoustic  phase  shift  (^)  given  by 

B  =  27rft  (3.28) 

is  then  maintained  as  a  constant  by  the  feedback  loop  and  therefore 

AB  =  0  =  27r(t  Af  +  f  At)  (3.29) 


and 


A^  _  Af  At 

B  f  t  ■ 


Thus,  as  stated  earlier, 


Af  ^  AW  ^  At 
f  W  t 


(3.30) 


(3.31) 


Therefore,  by  monitoring  the  change  in  the  frequency  of  the  P2L2,  the  change  in 
"natural"  velocity  is  also  monitored. 
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A  block  diagram  of  the  experimental  apparatus  used  for  the  uniaxial  SAC  meas¬ 
urements  is  shown  in  Fig.  3.2.  The  uniaxial  compressive  load  was  applied  as  in  the 
strain  gauge  measurements  with  a  MTS  810  load  frame.  The  calibrated  output  of  the 
load  cell  was  monitored  by  the  H.P.  3478A  digital  multimeter.  The  tone  burst  from  the 
P2L2  was  used  to  excite  a  conventional  damped  broad  band  2.25  MHz  transducer. 
Because  broad  band  transducers  were  used,  phase  shifts  at  the  transducer  specimen 
interface  were  assumed  to  be  negligible.  The  P2L2  was  used  in  a  pulse  echo  mode  with 
the  reflected  echos  input  back  into  the  instrument  for  phase  comparison.  A  Tektronix 
2445  oscilloscope  was  used  to  view  the  ultrasonic  signal,  the  phase  signal  and  the  sample 
and  hold  position.  It  was  triggered  by  a  reference  sync  from  the  P2L2.  The  frequency 
was  monitored  by  a  H.P.  5316A  universal  counter.  The  voltage  and  frequency  values 
were  read  by  a  Tektronix  4051  computer  and  converted  to  stress  and  normalized  fre¬ 
quency  shift.  The  program  listed  in  Appendix  G  also  computed  the  slope  of  the  curve 
by  linear  regression  as  well  as  plotted  and  stored  the  data. 

The  value  for  each  uniaxial  SAC  was  determined  by  averaging  the  values  from 
ten  measurements.  A  spring  loaded  clamping  device  was  used  to  maintain  constant 
pressure  of  the  transducer  on  the  specimen  in  order  to  hold  the  bond  thickness  constant 
during  the  measurement.  Variations  in  the  bond  thickness  would  cause  phase  shifts 
which  introduced  errors  into  the  measurement.  These  errors  which  have  less  effect  on 
thick  samples  made  initial  measurements  on  thin  laminates  impossible  and  placed 
another  constraint  on  minimum  specimen  thickness.  The  transducer  was  also  rebonded 
onto  the  specimen  three  or  four  times  during  the  ten  measurements  to  determine  the 
effect  of  slight  variations  in  initial  bond  thicknesses.  After  each  rebonding,  the  specimen 
was  ramped  through  the  load  cycle  several  times  to  settle  the  bond.  Then,  the  sample 
was  allowed  to  remain  at  rest  approximately  15  minutes  to  bring  the  sample  to  equili¬ 
brium  temperature.  This  was  necessary  because  the  handling  of  the  sample  during 
rebonding  increased  the  temperature  slightly  which  could  significantly  affect  the  results. 
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Another  concern  in  these  measurements  was  again  the  presence  of  viscoelastic 
effects  in  these  materials.  The  loading  rate  was  varied  considerably  to  examine  this 
effect.  The  specimen  were  loaded  to  between  80  and  100  MPa  in  times  ranging  from  10 
sec.  to  over  three  minutes  with  no  observed  effects.  The  measurements  were  made  for 
both  increasing  and  decreasing  load  with  no  measurable  differences.  The  effect  of  tem¬ 
perature  variations  due  to  thermoelastic  heating  was  also  considered.  Theory  assumes 
these  measurements  to  be  made  under  isothermal  conditions.  At  these  loading  rates 
which  were  very  nearly  isothermal,  the  effect  of  thermoelastic  heating  was  calculated  to 
be  much  less  than  experimental  uncertainty.  The  basis  of  these  calculations  was  the 
temperature  derivative  of  velocity  data  to  be  presented  later. 

The  previously  mentioned  constraints  on  the  choices  of  stress,  propagation  and 
polarization  directions  allows  for  18  measurements  to  be  made.  The  average  values 
from  the  ten  measurements  of  each  and  their  standard  deviations  are  presented  in 
Table  3.2.  The  uncertainties  were  typically  less  than  five  percent  of  the  measured 
values  for  these  measurements.  The  values  were  consistently  positive  indicating  increas¬ 
ing  frequency  as  a  function  of  increasing  compressive  stress  except  for  H323  and  H313  .  A 
few  of  the  wave  modes  exhibited  nonlinear  (quadratic)  curves  indicating  higher  than 
third  order  nonlinear  effects.  These  are  designated  by  an  *  in  the  table  and  the  SAC 
values  for  these  are  the  linear  coefficients  of  a  quadratic  fit  to  the  data.  Typical  curves 
for  the  uniaxial  SAC  measurements  are  shown  in  Figs.  3.3  -  3.20. 

As  with  the  linear  elastic  measurements,  some  of  these  measurements  can  be  used 
to  check  for  consistency  with  transversely  isotropic  behavior.  The  SAC  pairs 
(H122,  H2117)  (H311,  H322)  ,  (Hi33,  H233)  ,  H232)  ,  and  (H132,  H231)  should  all  have  equal 

values.  Other  pairs  that  should  have  equal  values  were  not  compared  as  they  were 
measured  on  different  specimens  making  comparisons  invalid  because  of  sample  to  sam¬ 
ple  variations.  The  values  of  some  of  these  pairs  compare  well  within  experimental 
uncertainty.  However,  H133  and  H233  differ  by  more  than  the  uncertainty  of  the 
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Figure  3.7  -  Delta  F/F  versus  stress  for  a  longitudinal  wave 
propagating  along  X2  with  stress  applied  along 


Figure  3.8  -  Delta  F/F  versus  stress  for  a  longitudinal  wave 
propagating  along  xj^  with  stress  applied  along  x 
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_  Delta  F/F  versus  stress  for  a  shear  wave  propagating 
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Figure  3.12  -  Delta  F/F  versus  stress  for  a  shear  wave  propagating 
along  X2  polarized  along  Xj^  with  stress  applied  along 

x^ . 
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;ure  3.15  -  Delta  F/F  versus  stress  for  a  shear  wave  propagating 
along  X2  polarized  along  with  stress  applied  along 
x^. 
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Figure  3.16 


Delta  F/F  versus  stress  for  a  shear  wave  propagating 
along  x-j^  polarized  along  x^  with  stress  applied  along 
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.19  -  Delta  F/F  versus  stress  for  a  shear  wave  propagating 
along  polarized  along  X2  with  stress  applied  along 
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Figure  3.20  -  Delta  F/F  versus  stress  for  a  shear  wave  propagating 
along  X3  polarized  along  X]^  with  stress  applied  along 
X2. 
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measursmsnt  as  v/eli  as  H131  and  H2S2  ■  Thus,  again  measurable  deviations  from 
transverse  isotropy  were  demonstrated.  Hov/ever,  why  this  behavior  is  demonstrated 
only  in  some  msasurements  and  not  in  others  is  not  clear. 

The  hydrostatic  measurements  were  made  in  a  similar  fashion  as  those  made 
under  uniaxial  compression.  A  block  diagram  of  the  set  up  is  shown  in  Fig.  3.21.  The 
change  in  "natural”  velocity  v/as  again  measured  with  the  P2L2.  Argon  was  used  to 
increase  the  pressure  v/ithin  the  pressure  chamber.  A  FLP.  5316A  universal  counter 
monitored  the  frequency  while  a  Tektronix  465B  oscilloscope  v/as  used  to  view  the  ultra¬ 
sonic,  phase,  and  sample/hold  signals.  The  Tektronix  4051  computer  was  used  to  read 
the  frequency  data.  The  pressure  data  v/as  manually  read  from  a  Heise  C-57488  pres¬ 
sure  gauge  and  input  into  the  computer. 

There  v/ere  several  differences  from  the  uniaxial  measurements,  however.  While 
temperature  variations  during  the  uniaxial  SAC  measurements  were  insignificant,  they 
were  important  in  the  hydrostatic  measurements.  Increasing  the  pressure  25  psi.  sud¬ 
denly,  caused  an  increase  in  temperature  of  the  gas  and  thus  the  sample  of  several 
tenths  of  a  degree  C.  This  temperature  increase  was  significant.  Therefore,  the  pres¬ 
sure  had  to  be  increased  slowly  and  then  time  for  thermal  equilibrium  to  be  obtained 
had  to  be  given.  Typically,  there  was  a  20  to  30  minute  delay  between  data  points. 
Also,  O'/er  this  period  of  time  the  room  temperature  could  vary  significantly.  This 
caused  additional  problems.  To  overcome  this,  heating  elements  were  installed  in  the 
piessu^e  cnamber.  A  cemperature  controller  v/as  used  with  these  heaters  to  maintain 
the  temperature  several  degrees  C  above  room  temperature.  This  speeded  the  time 
required  to  attain  thermal  equilibrium  and  avoided  the  effect  of  room  temperature  vari¬ 
ations  provided  it  did  not  rise  more  than  the  several  degree  buffer.  A  Lake  Shore  Cryo- 
tronics  Inc.  model  DRC84C  temperature  controller  was  used.  This  instrument  had  a 
reported  ability  to  maintain  the  temperature  v/ithin  -l-/-  0.1  degree  C.  Although  this 
inscrument  drastically  improved  the  situation,  these  small  variations  in  temperature 
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Figure  3.21  -  Block  diagram  of  apparatus  used  in  the  hydrostatic 
pressure  stress  acoustic  constant  measurements. 


reduced  the  accuracy  of  the  measurement  significantly. 


Another  problem  that  increased  the  uncertainty  in  the  hydrostatic  measurements 
urns  the  limitation  in  pressure.  The  upper  limit  of  pressure  in  the  chamber  was  set  at 
only  250  psi.  This  was  about  a  factor  of  50  less  than  the  maximum  stress  applied  dur¬ 
ing  the  uniaxial  measurements.  Theoretically^  the  slope  of  the  SAC  should  be  eyaluated 
at  zero  stress  and  thus  this  limitation  in  pressure  should  not  be  a  problem.  However,  the 
low  pressure  produced  much  smaller  velocity  changes  making  them  more  uncertain. 
Also,  the  effect  of  temperature  variations  were  much  greater.  Thus,  the  uncertainty  of 
these  measurements  v/as  much  larger.  It  was  on  the  order  of  ten  percent  of  the  meas¬ 
ured  SAG  for  most  of  the  hydrostatic  measurements.  The  limitation  in  maximum,  pres¬ 
sure  also  made  the  measurement  of  ¥±q^  impossible  because  the  changes  involved  v/ere  to 
small  to  detect  v/itli  the  P2L2.  The  possibility  of  anomalous  behavior  at  low  pressures 
must  also  be  considered.  Therefore,  although  the  hydrostatically  measured  SAC  values 
v/ere  used  with  the  uniaxial  SA_C  measurements  to  determine  the  nonlinear  moduli  in 
this  work,  measurements  at  higher  pressures  are  needed  to  determine  the  validity  of 
these  calculations. 

Because  of  the  tremendous  amount  of  time  involved  in  making  a  single  hydros¬ 
tatic  SAC  measurement  (10  to  15  data  points  at  20  to  30  minutes  per  data  point),  only 
three  measurements  v/ere  used  in  the  calculation  of  each  hydrostatic  SAC.  This  also 
contributed  to  the  large  uncertainty  in  these  measurements.  Measurements  v/ere  again 
made  for  increasing  and  decreasing  pressures  with  no  observed  differences. 

The  hydrostatic  measurements  v/ere  all  made  using  undamped  crystal  transduc¬ 
ers  (PZT  5A)  because  of  size  limitations  v/ithin  the  pressure  chamber.  This  created  the 
necessity  of  making  corrections  because  of  phase  shifts  at  the  transducer-specimen  inter¬ 
face.  These  phase  shifts  are  a  function  of  frequency  and  can  be  calculated  using  the 
theory  presented  by  McSkimmin  [l9j.  The  correction  to  the  experimentally  determined 
normalized  frequency  shift  assuming  zero  bond  thickness  was  given  by 
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F.  ""  F 

true  exp 


exp.  ^trans. 


where  — - was  the  true  frequency  change, 

F^true 


was  the  experimentally  measured  change, 


and - was  the  change  in  time  of  flight  due  to  phase  shifts. 

^trans. 

The  term  was  given  by 

ttrjins.  r  *  .  1 


Fi(360)  Fo(360) 


where  Fq,  were  the  measured  frequencies  of  two  data  points, 
to  was  the  time  of  flight  in  the  unstressed  specimen, 
and  7o,  7i  were  the  phase  angles  of  the  two  data  points  in  degrees. 


They  were  given  by 


.  2Z.Ro 

lo  -  tan  ;  , 

^0  “ 


(3.34) 


=  tan 


_1 

Rf  -  zi 


where  Zg,  were  the  mechanical  impedances  of  the  specimen  and 
transducer  respectively, 


Rq  =  Zx  tan  -  •  180 

Fp 


and 
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Ri  =  tan 


180 


(3.37) 


where  Fr  was  the  resonant  frequency  of  the  transducer. 

The  effect  of  nonzero  bond  thickness  on  this  correction  term  ¥/as  also  analyzed  and 
found  to  be  less  than  experimental  uncertainty. 

The  original  and  corrected  values  for  the  hydrostatic  SAC’s  are  shown  in  Table 
3.3.  Typical  curves  are  shov/n  in  Figs.  3.22  -  3.29.  Comparisons  of  the  data  to  check  for 
transversely  isotropic  behavior  shov/  no  measurable  differences.  This  is  most  likely  due 
to  the  large  uncertainty  of  the  measurements. 

The  effect  of  temperature  on  the  ’’natural”  velocity  was  also  measured  for  small 
variations  about  room  temperature.  The  apparatus  used  for  these  measurements  was 
the  same  as  that  used  in  the  hydrostatic  measurements.  The  temperature  was  increased 
above  room,  temperature  several  degrees  C  in  0.3  C  steps  using  the  heaters  and  the  tem¬ 
perature  controller.  The  frequency  v/as  measured  at  each  point  after  thermal  equili¬ 
brium  was  attained.  Again  corrections  were  applied  to  account  for  phase  shifts  at  the 
transducer-specimen  interface.  Typical  curves  are  shov/n  in  Figs.  3.30  -  3.38.  The  slopes 
of  these  curves  were  designated  as  the  Thermal  Acoustic  Constants  and  are  given 

in  Table  3.4  where  the  second  and  third  subscripts  give  the  directions  of  propagation 
and  polarization  respectively.  The  same  samples  used  for  the  hydrostatic  and  uniaxial 
measurements  v/ere  also  used  for  these  measurements. 

The  curves  all  demonstrated  linear  relationships  between  temperature  change 
and  normalized  frequency  shift  for  the  limited  range  of  the  measurements.  ¥7hile  these 
measurements  were  not  used  in  the  data  reduction  to  determine  the  third  order  elastic 
coefficients,  they  were  important  in  determining  the  effect  of  temperature  changes  on 


the  SAC  measurements. 
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Table  3.3  -  Measured  and  corrected  values  for  the  Hydrostatic  Stress 
Acoustic  Constants. 


AF/F  ",  T  AF/F 
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stress  (Pq) 

5.24  -  Delta  F/F  versus  pressure  for  a  shear  wave  propa¬ 
gating  along  Xj^  polarized  along  X2. 


STRESS  (Pa) 

Figure  3.25  -  Delta  F/F  versus  pressure  for  a  shear  wave  propa¬ 
gating  along  X2  polarized  along  Xj^. 


Figure  3,26  -  Delta  F/F  versus  pressure  for  a  shear  wave  propa^ 
gating  along  X]_  polarized  along 
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Figure  3.27  -  Delta  F/F  versus  pressure  for  a  shear  wave  propa 
gating  along  Xt  polarized  along  X3. 
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Figure  3.28  -  Delta  F/F  versus  pressure  for  a  shear  wave  propa¬ 
gating  along  polarized  along  Xp 
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Figure  3.29  -  Delta  F/F  versus  pressure  for  a  shear  wave  propa¬ 
gating  along  x^  polarized  along  X2. 
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Constants,, 


Measured  and  corrected  values  for  the  Thermal  Acoustic 


Thermal  Acoustic  Constants  C"') 


Original  Value 


-7.58  -1-/-  0.02 


-7.55  0.03 


-9.9  +/-  0.1 


-9.9  0.1 


-11.9  -}-/-  0.2 


-11.8  -5/-  0.2 


-8.4  -f/-  0.5 


-9.0  +/-  0.7 


-0.84  4-/-  0.09 


Corrected  Value 


-6.36  -f/-  0.02 


-6.33  -h/-  0.03 


-8.5  -1-/-  0.1 


-8.5  -h/-  0.1 


-10.2  +/-  0.2 


-10.1  +/-  0.2 


-7.1  +/-  0.5 


-7.6  +/-  0.7 


-0.72  +/-  0.09 


AF/F 


Figure  3.30  -  Delta  F/F  versus  temperature  change  for  a  longitudinal 
wave  propagating  along  Xj^. 
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Figure  3.31  -  Delta  F/F  versus  temperature  change  for  a  longitudinal 
wave  propagating  along  X2. 


A  F/F  c  A  F/F 


-113- 


AT  (C) 


0.0  0.5  1.0  1.5  2.0  2.5 


re  3o34  -  Delta  F/F  versus  temperature  change  for  a  shear 
wave  propagating  along  polarized  along 
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Figure  3.35  -  Delta  F/F  versus  temperature  change  for  a  shear 
wave  propagating  along  X2  polarized  along  Xj^. 
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Figure  3.38  -  Delta  F/F  versus  temperature  change  for  a  longitudinal 
wave  propagating  along 
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The  SAC  data,  along  -with  the  linear  moduli  ¥/ere  used  to  determine  some  of  the 
third  order  elastic  moduli  for  this  material.  Since  there  ¥/ere  more  equations  and  data 
than  unkno¥/n  moduli,  the  data  ¥/as  reduced  by  using  a  least  squares  fit  to  determine 
the  best  values  for  the  moduli.  The  data  was  fit  using  both  the  transversely  isotropic 
and  orthotropic  models.  Hov/ever,  the  lack  of  data  for  H033  prevented  the  calculation  of 
the  moduli  and  C33S  for  the  transversely  isotropic  model  and  the  moduli 
^133;  ^233)  ^^33?  <^123;  ^453  for  the  orthotropic  m>odel. 

The  least  squares  procedure  for  the  determination  of  third  order  moduli  was  out¬ 
lined  by  Hankey  and  Schuele  [40].  The  formulation  is  as  follows.  The  measured  values 
for  the  SAC’s  are  added  to  the  constant  term.s  in  the  equations  of  Table  3.1.  These  are 
designated  Ivl^  .  The  coefficients  of  the  equations  are  designated  A^j  and  the  moduli  to 
be  determined  are  Dj.  Thus  the  equations  of  Table  3.1  are  rewritten  as 

Mj  =  AjjDj  (3.38) 

wffiere  i  is  summed  over  the  number  of  data  used  and  j  is  summed  over  the  number  of 
moduli  to  be  determined.  This  can  be  rev/ritten  as 


A,kMi  =  (A,kAij)Dj  (3.39) 

where  k  also  runs  from  one  to  the  number  of  coefficients  to  be  determined.  The  best  fit 
for  the  moduli  Dj  are  then  given  by 

=  (AikAjj)  ^Apj^Mp  (3.40) 

v/here  (Ajj^Ajj)"^  is  the  inverse  matrix  element  and  p  is  summed  over  the  number  of 
measured  values. 

If  it  is  assumed  that  the  error  in  the  measured  values  is  much  greater  than  the 
error  in  the  coefficients,  then  the  error  in  the  coefficients  can  be  dropped.  The  error  in 


the  fitted  values  is  then  given  by 


AD|  =  |(A|kA,jr‘Apfc|AMp 


where  ADj  are  the  maximum  limits  of  error, 
and  AMp  are  the  errors  in  the  measured  values. 
The  probable  limit  of  error  (ADj)  is  given  by 


(N  -  if 


(3.41) 


(3.42) 


where  N  is  the  number  of  measurements. 

The  above  stated  assumption  was  valid  in  these  measurements  as  the  errors  in  the 
measured  SAC’s  ranged  from  one  to  ten  percent  while  the  errors  in  the  linear  coefficients 
were  typically  less  (parts  in  10^) 

Programs  were  written  to  perform  the  data  reduction  for  the  two  symmetry 
models  and  are  listed  in  Appendix  H.  The  fitted  values  and  their  uncertainties  are 
given  in  Table  3.5. 


in.E.  Conclusion 

Using  a  pulsed  phase  locked  loop  interferometer,  the  normalized  change  in  the 
ultrasonic  "natural"  phase  velocity  as  a  function  of  both  stress  and  temperature  have 
been  measured.  The  stress  measurements  were  made  under  conditions  of  hydrostatic 
and  uniaxial  stress.  The  measured  curves  were  predominantly  linear  although  several 
were  nonlinear  (quadratic)  indicating  higher  than  third  order  nonlinear  elastic  effects. 
Because  of  equipment  limitations,  the  hydrostatic  measurements  were  made  at  much 
lower  pressures  in  comparison  to  the  uniaxial  measurements.  Thus  the  uncertainties  of 
the  hydrostatic  measurements  were  much  larger  making  the  validity  of  using  them  with 
the  uniaxial  measurements  to  calculate  the  third  order  moduli  questionable.  This 
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Table  3.5  -  Continued 


Orthotropic  Nonlinear  Moduli  (Continued) 

Modulus 

Value 

Max.  Error 

Probable  Error 

^222 

-186 

21 

4 

^223 

-60 

93 

19 

^144 

-33.0 

2.4 

0.5 

^244 

-47.8 

3.3 

0.7 

<^344 

-46 

27 

6 

^156 

-50.1 

3.9 

0.8 

^256 

-33.5 

2.8 

0.6 

^366 

-49 

32 

7 

^166 

-33.9 

0.6 

0.1 

<^266 

-33,1 

0.6 

0.1 

^366 

-28.5 

4.0 

0.8 
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limitation  in  maximum  hydrostatic  pressure  also  prevented  the  measurement  of  H033 
v/hich  prevented  the  calculation  of  several  of  the  nonlinear  moduli.  The  hydrostatic 
measurements  need  to  be  evaluated  at  higher  pressures  to  reduce  the  uncertainties  and 
to  evaluate  H033  .  This  would  enable  the  determination  of  the  remaining  third  order 
moduli. 

The  effect  of  temperature  changes  on  "naturaF’  velocity  was  also  evaluated  for 
nine  combinations  of  propagation  and  polarization  directions.  These  measurements 
v/ere  made  over  a  narrow  range  of  temperature  change  slightly  above  room  tempera¬ 
ture.  The  curves  all  demonstrated  linear  relationships  betv/een  velocity  and  tempera¬ 
ture.  Although  not  directly  used  for  the  evaluation  of  third  order  moduli,  these  meas¬ 
urements  v/ere  important  in  the  evaluation  of  the  effects  of  temperature  variations  dur¬ 
ing  the  SAC  measurements. 

The  SAC  values  and  the  second  order  moduli  were  used  to  calculate  some  of  the 
third  order  elastic  moduli.  A  least  squares  fit  was  applied  for  both  the  transversely  iso¬ 
tropic  and  the  orthotropic  models  of  elastic  symmetry.  The  values  for  these  moduli 
were  all  negative  with  the  exception  of  the  Ci2s  •  This  is  in  agreement  with  the  non¬ 
linear  behavior  of  conventional  materials.  Although  the  linear  elastic  moduli  and  the 
SAC  data  demonstrated  small  deviations  from  transverse  isotropy,  large  uncertainties  in 
the  third  order  moduli  masked  any  su'ch  deviations.  The  orthotropic  moduli  agree 
v/ithin  experimental  error  to  the  conditions  of  transverse  isotropy. 
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IV.  SUMMARY  AND  CONCLUSIONS 


The  mechanical  behavior  of  a  unidirectional  T300/5208  graphite/epoxy  compo¬ 
site  material  has  been  evaluated  in  this  research.  In  particular,  the  linear  (second 
order)  and  nonlinear  (third  order)  elastic  properties  have  been  measured.  The  linear 
moduli  were  measured  by  ultrasonic  velocity  measurements  as  well  as  by  static  loading 
strain  gauge  measurements.  The  nonlinear  moduli  were  calculated  from  the  linear 
moduli  and  measurements  of  the  change  in  ultrasonic  "natural"  phase  velocity  as  a 
function  of  stress. 

The  measured  linear  elastic  properties  of  this  material  agreed  well  with  previ¬ 
ously  reported  values  for  a  similar  material  except  for  the  moduli  C33  and  S33  .  These 
were  the  moduli  that  reveal  the  behavior  of  the  composite  under  loading  along  the  fiber 
direction.  The  low  value  for  C33  and  corresponding  high  value  for  S33  may  have  been  a 
result  of  a  lower  fiber  content  or  higher  void  content  in  this  material.  It  could  also  have 
been  the  effect  of  misalignment  of  fibers  in  some  of  the  laminae.  However,  the  values  of 
the  ultrasonic  and  static  measurements  compared  well  with  the  exception  of  the  off 
diagonal  modulus  C13  .  This  discrepancy  was  probably  the  result  of  either  large  uncer¬ 
tainties  of  the  static  measurements  or  error  propagation  in  the  inversion  of  the  ultra¬ 
sonic  moduli.  These  measurements  also  pointed  out  measurable  deviations  from 
transversely  isotropic  elastic  behavior  which  is  usually  assumed  for  unidirectional  com¬ 
posites.  The  material  more  correctly  possessed  orthotropic  symmetry  although  the  devi¬ 
ations  from  transverse  isotropy  were  small. 

The  change  in  ultrasonic  "natural"  velocity  as  a  function  of  stress  and  tempera¬ 
ture  was  measured  with  a  pulsed  phase  locked  loop  ultrasonic  interferometer  for 
numerous  wave  modes.  These  measurements  again  demonstrated  deviations  from 
transversely  isotropic  behavior.  Also,  some  of  the  stress-velocity  curves  were  nonlinear 
(quadratic)  which  indicated  higher  than  third  order  nonlinear  effects.  The 
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temperature-yelocit}’  curyes  wsre  a!l  linear  oyer  the  small  range  of  the  measurements. 
Some  of  the  third  order  elastic  stiffness  moduli  vvere  calculated  from  the  stress  depen¬ 
dence  of  velocitj^  These  lyere  all  negatiye  in  agreement  with  nonlinear  properties  of  con¬ 
ventional  materials.  The  uncertainties  of  some  of  these  moduli  were  large  as  a  result  of 
the  large  uncertainties  of  the  hydrostatic  pressure  measurements  propagating  in  the 
least  squares  data  reduction.  These  large  uncertainties  masked  any  deviations  from 
transverse  isotropy  that  might  have  been  present. 

However,  it  should  be  noted  here  that  the  measured  SAC’s  are  probably  more  use¬ 
ful  than  the  calculated  nonlinear  moduli.  This  is  because  they  are  more  easily  obtained 
v/ith  snaaller  experimental  uncertainty  and  can  be  measured  in  an  individual  specimen 
instead  of  calculated  from  meas?urernents  made  in  several  specimens.  AAso,  the  SAC’s  are 
the  quantities  that  are  more  useful  in  experimental  attempts  to  measure  applied  and 
residual  stress  and  may  be  important  in  nondestructive  evaluation  of  ultimate  strength 
of  materials. 


It  v/ouid  be  remiss  to  fail  to  point  out  several  factors  that  could  be  improved  to 
yield  better  measurements.  First,  the  hydrostatic  measurements  should  be  carried  out 
to  higher  pressures  v/ith  improved  tem.perature  controller  capabilities.  This  would 
reduce  the  uncertainties  of  these  measurem.eDts  and  thus  those  of  the  nonlinear  moduli. 
Also,  the  effect  of  sam.ple  to  sam.pie  variations  on  both  the  linear  and  nonlinear  proper¬ 
ties  needs  to  be  evEiluated  further.  These  measurements  v/ere  made  on  as  few  samples 
as  possible  to  reduce  this  effect  in  the  calculations.  Future  measurements  should  also 
evaluEdus  the  changes  in  velocity  over  v/ider  temperature  variations  and  v/ith  much 
Higher  stresses  (i.e.  near  ultimate  strength  stresses)  as  well  as  evaluate  the  nonlinear 
properties  hj  other  techniques  such  as  harmonic  generation. 

'chough  the  primary  goal  of  this  investigation  v/a.s  a  more  thorough  characteri- 
zacion  of  tne  mechcmical  properties  of  graphite/epoxy  composites,  these  measurements 
serve  as  a  basis  for  further  study  and  potential  applications  in  the  area  of 
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nondestructive  testing  of  these  materials.  Also,  the  effects  of  fiber  volume  ratio  and 
more  complex  lay  up  orientations  on  nonlinear  properties  of  composites  need  to  be 
evaluated.  The  relationships  between  nonlinear  properties  and  important  engineering 
properties  such  as  strength,  residual  strength  after  impact  and  fatigue  loading,  and 
fiber-matrix  interfacial  strength  need  to  be  studied.  These  may  lead  to  the  development 
of  needed  nondestructive  evaluation  techniques  for  composite  materials  which  will  play 
an  important  role  in  the  future  of  the  aerospace  industry. 

The  problem  or  residual  strain  in  composites  which  arises  because  of  the 
mismatch  of  coefficients  of  thermal  expansion  between  fiber  and  matrix  is  another 
potential  application  for  these  measurements.  When  the  composite  material  is  cooled 
from  the  higher  temperatures  of  polymerization  during  cure,  large  strains  develop  which 
may  degrade  the  ultimate  strength  of  the  material.  Although  special  curing  cycles  have 
been  developed  to  reduce  this  effect,  it  still  remains  as  an  important  problem  in  the  field 
of  composites.  And  as  nonlinear  ultrasonics  have  been  used  in  attempts  to  nondestruc- 
tively  evaluate  residual  strains  in  conventional  materials,  they  may  well  be  useful  in  the 
measurement  of  residual  stains  in  composites.  Especially,  since  it  may  be  possible  to 
characterize  the  elastic  symmetry  of  the  composite  better  than  in  conventional  materi¬ 
als  because  the  symmetry  of  the  material  is  determined  by  the  lay  up  orientation.  In 
conventional  materials  the  problems  associated  with  anisotropic  textures  of  the  material 
have  limited  the  applicability  of  this  technique. 
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APPENDIX  A 


The  derivation  of  the  independent  linear  (second  order)  elastic  moduli  for  the 
symmetries  of  transverse  isotropy  and  orthotropy  are  now  presented.  The  case  of  ortho- 
tropic  elastic  symm_etry  is  presented  first.  In  an  orthotropic  material,  there  are  three 
orthogonal  planes  of  two  fold  symmetry.  Therefore  if  the  three  coordinate  axes  are 
chosen  such  that  they  are  the  normals  to  these  planes,  then  the  independent  elastic 
moduli  can  be  determined  by  applying  the  tv/o  fold  rotations  about  each  of  these  axes 
on  the  strain  energy  function.  The  strain  energy  function  is  initially  given  by 


J- 


^22^2 


+  C 


33^3 


“r 


-55^5 


-1-  C 


68^6, 


-1-  Ci3ei€3  -i-  0146184  +  0156165  +  0156165  +  0336263  +  0346064  -h  0356365 

-1-  0356365  4-  0346364  -f  0356365  -1-  0356365  -f  0456465  -I-  0456465  -h  0556565.  (A.l) 


The  effect  of  a  180  degree  rotation  about  the  X2  axis  is  nov/  considered.  The 
matrix  of  transformation  is  given  by 


-10  0 
0  1  0 
0  0-1 


(A.2) 


The  efect  on  this  rotation  on  the  strains  is  given  by 


(A.3) 


which  yields  the  following  relations  between  the  rotated  and  original  strains  expressed 
in  the  Voigt  reduced  notation 
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The  strain  energy  function  can  be  written  for  the  rotated  state  and  the  above  equations 
substituted  to  yield 

=  — (Ciief  +  022^2  +  <^33^3  ^44^4  +  ^55^5  +  ^66^1)  +  ^12^1^2 

2 

+  ^13^1^3  “  +  *^15^1^5  “  ^166165  +  023^2^3  —  C24e2e4  +  ^2562^5 

—  C26e2^6  “  <^34^3^4  '^35^3^5  —  ^36^3^6  “  ^456465  +  C46e4e6  “*  ^56^5^6-  (■^•^) 

Equating  ^  and  term  by  term  with  respect  to  strain  yields  the  following 


Ci4  —  Ci4 

C16  =  -C16 

C24  =  “^24 
^26  ~  ~^26 
C34  ==  “^34 
C36  ~  "“^36 
C45  =  ~‘^45 
C56  ~  "“<^56 


(A.6) 


For  these  relations  to  be  true,  then 

<^14  ~  *^16  ~  ^24  “  ^26  “  ^34  ~  ^36  “  *^45  “  ^56  “ 


(A.7) 


The  matrix  of  elastic  stiffness  moduli  is  then  given  by 
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Cab]  “ 


1 C 1  -j 

Ci2 

Cl3 

,0 

'O 

Cl5 

0  ' 

CT2 

£'22 

C23 

0 

Cos 

0  i 

Cl3 

C23 

C33 

■n.. 

U 

C35 

0 

0 

0 

0 

C44 

0 

C46: 

Clo 

C25 

C35 

0 

C55 

0  ^ 

0 

0 

0 

C46 

0 

Cse; 

(A.1 


The  strain  energ;/-  function  is  non/  given  by 

”  ylCiUl"  <^22^2  C33  -I”  £3  -I-  c|  -1-  C55£|  -j-  e|)  V  Ci2eie2 

-]-  G-iS^i^s  -I-  -'r  C23C2e3  -r  C25£2C5  “1“  C35C3e5  -j- 


(A.8) 


ihs  enect  01  a  180  degres  rotation  about  the  Xj  axis  can  now  be  considered.  The 
transformation  matrix  is  given  by 


a,i  = 


10  0  S 

0-1  0 

0  0  -l! 


(A.9) 


The  rotated  and  original  strains  are  related  by 


(A.IO) 


Again  these  are  substituted  into  the  equation  for  the  strain  energy  function  resulting  in 
1 


n  c 


C02C2  -r  €3363  C44C4  +  C55C5  -b  CqqCq)  -b  01262^2 
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+  ^13^1^3  ~  <^15^1^5  +  ^23^2^3  ~  <^25^2^5  ~  ^356365  —  0456465 

Equating  term  by  term  yields 

^15  =  <^25  =  C35  =  C45  =  0 
and  reduces  the  elastic  stiffness  matrix  to 

^11  ^12  ^13  0  0  0 

^12  ^22  ^23  0  0  0 

^13  ^^23  ^33  0  0  0 

I'abI  =0  0  0  c„  0  0 

0  0  0  0  C55  0 

0  0  0  0  0  Cge 

The  third  symmetry  rotation  about  X3  can  be  shown  to  have  no  further  effect  on  the 
elastic  stiffness  tensor  and  thus  the  final  form  for  an  orthotropic  body  is  given  by  equa¬ 
tion  A. 13.  There  are  nine  independent  moduli  for  an  orthtropic  material. 

The  reduction  of  moduli  for  a  transversely  isotropic  material  can  be  carried  out 
beginning  with  the  matrix  form  for  orthotropy.  This  is  because  it  possesses  the  same 
symmetry  conditions  as  an  orthotropic  material  in  addition  to  having  an  isotropic  plane 
perpendicular  to  X3  .  Thus  the  strain  energy  function  can  be  written  out  with  the  ortho- 
tropic  moduli  and  the  effect  of  an  arbitrary  rotation  around  the  X3  axis  can  be  con¬ 
sidered.  ^  is  given  by 

^  ~  +  ^22^1  +  63363  +  C446I  +  C5565  -1-  65565) 

+  +  Ci3ei63  +  6236263.  (A.14) 

The  effect  of  an  arbitrary  rotation  around  X3  by  an  angle  9  on  the  strains  is  given  by 


(A.11) 


(A.12) 
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£i  =  cos^^Ei  sin^cos^Cg  -1-  sin^0e2 

4  =  sin^^'ci  4-  sin^'cos^eg  -F  003^6*62 

/ 

^3  =  ^3 

64  =  — sin^Cg  4”  0036^64 
£5  =  cos^fig  —  sin^e4 

€5  =  — 23in^co36*ej  -1-  (cos^6*  —  3iii^6*)e5  -F  23in^co3^e2- 


(A.15) 


These  equations  can  be  substituted  into  the  strain  energy  function  as  before.  Then  com¬ 
paring  <P  and  term  by  term  will  yield 


^11  —  <^22 


^13  “  ^23 
C44  =  Cgg 

^66  =  ~  ‘^12) 


(A.16) 


which  reduces  the  number  of  independent  elastic  moduli  to  five  for  a  transversely  isotro¬ 
pic  body.  The  elastic  stiffness  matrix  can  be  written  out  as 


.^ABi 


- 

Cll 

C12 

<^13 

0 

0 

0 

^12 

^11 

Cl3 

0 

0 

0 

^13 

<^13 

C33 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

C44 

0 

0 

0 

0 

0 

0 

— ■ . 

0 

1 

vhh 

T - 1  j 

(A.17) 
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APPENDIX  B 

Using  the  equations  derived  for  linear  elastic  wave  propagation  presented  in 
Chapter  II.,  the  expressions  for  the  velocities  of  ultrasonic  waves  in  terms  of  elastic 
moduli  and  density  are  now  derived.  The  expressions  are  derived  for  both  transverse 
isotropy  and  orthotropy.  The  case  of  transversely  isotropic  elastic  symmetry  is  con¬ 
sidered  first. 

For  the  case  of  transverse  isotropy,  equation  (2.84)  can  be  used  to  show  that 
there  are  several  directions  of  propagation  for  pure  mode  elastic  waves.  Pure  mode 
elastic  waves  may  be  propagated  along  the  X3  axis  and  along  any  direction  perpendicu¬ 
lar  to  the  X3  axis  in  the  isotropic  plane  which  of  course  includes  the  Xj  and  X2  axes. 
They  may  also  be  propagated  in  a  direction  which  makes  an  angle  0  with  the  X3  axis 
given  by 

r  -,1 


tan  0  = 


<^13  +  ^^44  ~  <^33  2 
Ci3  +  2C44  —  C41 


(B.l) 


However,  for  unidirectional  graphite  epoxy  composites  where  the  modulus  C33  is  usually 
about  two  orders  of  magnitude  larger  than  the  other  moduli,  this  direction  does  not 
exist.  Thus  the  pure  mode  directions  to  be  considered  are  along  each  of  the  three  coor¬ 
dinate  axes. 

The  components  of  the  matrix  Xjj^  for  a  transversely  isotropic  material  are  given 
by 

^11  ~  “  ‘^12)11  +  ^44^3- 

^22  =  ~  +  ^44l|> 


^33  ~  ^44^2  <^33^3  > 


(B.2) 
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^12  ~  ^21  ~  2  '^12)^1^2> 
^13  =  ^3i  ~  (<^13  +  ‘^44)h^3> 


and 


^23 


X. 


32 


(c 


44 


is)^ 


13;^2^3- 


The  solutions  for  the  y/ave  modes  propagating  along  the  axis  are  novy  derived.  For 
vvaves  propagating  in  this  direction,  the  direction  cosines  of  the  wave  normal  are  given 
by  li  =  1,  and  I2  =  I3  =  0  .  The  components  of  the  Xy  matrix  reduce  to 

^11  ~  ‘^ii> 

^22  ~  ~ 


X 


33  —  C44, 


and  all  other  Xy  =  0  .  These  are  substituted  into  equation  (2.55)  to  determine  the  solu¬ 
tions  for  the  three  modes  of  propagation  in  this  direction  by  solving  the  equation 


DET 


cii  -  PqV 
0 
0 


0 


•(cii  —  C12)  “  Po^“ 


0 


0 

0 

C44  -  PqV^ 


0. 


(B.4) 


This  can  be  expanded  to 

(cii  -  Pov^)(y(cii  -  C12)  -  Pov^)(c44  -  =  0-  (B.5) 

2 


The  solutions  are  then  found  to  be 

2. 
2 
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Solving  for  the  eigenvectors  for  each  of  these  eigenvalues  yields  the  following  particle 
displacements.  For  the  first  wave,  the  direction  cosines  of  the  particle  displacements  are 
^2  ==  0  =  1  .  Therefore,  since  ^1  =  1^,  the  first  wave  is  a  pure  mode  longitu¬ 

dinal  wave.  For  the  second  wave,  the  direction  cosines  are  =  0^3  =  0  and  0^2  =  1  . 
Thus,  since  ajli  =  0  ,  it  is  a  pure  mode  transverse  wave.  The  particle  displacement  vec¬ 
tor  direction  cosines  for  the  third  wave  are  found  to  be  =  ^2  =  0  and  ^3  =  1  ,  which 
means  it  too  is  a  pure  mode  shear  wave.  The  energy  flux  equations  can  be  used  to  show 
that  there  is  no  energy  flux  deviation  for  these  waves. 

The  solutions  for  propagation  along  the  X2  axis  are  similar.  The  direction  cosines 
of  the  wave  normal  are  l^  =  I3  =  0,  and  I2  =  1  .  The  Xjj  components  are  then 

^11  ~  ~ 

^22  =  C,,,  (B.6) 

^33  =  ^44’ 

and  all  other  Xjj  =  0  .  The  solutions  for  these  wave  modes  are  of  course  the  same  as  for 
propagation  along  x^  except  the  first  and  second  waves  are  reversed.  Waves  propagat¬ 
ing  in  any  direction  in  the  Xj,  X2  plane  with  I3  =  0  can  be  shown  to  have  similar  solu¬ 
tions  to  the  waves  propagating  along  Xj  and  X2  . 

Propagation  along  X3  where  1^  :=  I2  =  0,  and  I3  =  1  is  now  considered.  The  com¬ 


ponents  of  Xjj  are  now  given  by 
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^11  ~  ^44> 

^22  ~  ^44> 

^33  =  ^33’ 

and  all  other  =  0  .  Substitution  into  equation  (2.55)  yields 

(C44  -  P0'^^)^(C33  -  P0^~)  =  0- 


(B,7) 


{B.8) 


The  solutions  for  this  equation  are  given  by 


/  \ 


Vi  =  Vn 


-44 


Po  J 


(B.9) 


andj 


Vo  = 


^33 

Po 


Solutions  for  the  particle  displacements  for  these  wave  modes  show  that  for  the  first  and 
second  wave  modes,  0^3  =  0  while  and  are  arbitrary  as  long  as  of  -f  a|  =  1  .  Thus 
these  v/aves  are  pure  mode  shear  waves  with  particle  displacements  in  the  X2  plane. 
For  the  third  wave,  they  are  ^  ~  03  =  1  ,  which  means  it  is  a  pure  mode 

compressional  wave.  Again,  it  can  be  shown  that  these  waves  suffer  no  energy  flux  devi¬ 
ation. 

Measurements  of  the  velocities  of  the  waves  in  these  three  directions  can  be  used 
to  determine  four  of  the  five  independent  elastic  stiffness  moduli.  However,  the  modulus 
Ci3  can  not  be  determined.  To  evaluate  this  modulus,  an  off  axis  non-pure  mode  wave 
must  be  used.  The  wave  chosen  for  this  purpose  was  a  v/ave  propagating  in  the  Xj,  X3 
plane  45  degrees  between  the  two  axes.  Thus,  l^  =  I3  =  ,  and  I2  —  0  .  The  com¬ 


ponents  of  the  Xjj  matrix  are  then 
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^11  =  +  ‘^44). 

^22  ~  "*”  ~ 

>^33  =  -^(^33  +  C44),  (B.IO) 

^13  =  ■j(‘^13  +  C44), 


and  all  other  X^j  =  0  .  These  are  substituted  into  the  eigenvalue  equation  and  the  solu¬ 
tions  are  given  by 

Po^l  =  [^11  +  2C44  +  C33  +  'V(cii  —  033)^  +  4(ci3  +  044)^1 


2  1 
P0V2  =  - 


C44  +  ^(<^11  “  Ci2)|, 


(B.ll) 


and 


P0V3 


~  4  ~  V^(^ll  ~  *^33)^  +  4(Ci3  +  044)^1 


Solutions  for  the  particle  displacements  of  these  waves  show  that  the  second  wave  is  a 
pure  mode  shear  wave.  However,  the  first  wave  is  a  quasilongitudinal  wave  and  the 
third  wave  is  quasishear.  Energy  flux  equations  show  that  these  wave  modes  also  suffer 
energy  flux  deviation.  Kriz  and  Stinchcomb  [45]  reported  the  energy  flux  deviation 
could  be  as  high  as  43.4  degrees  for  the  quasilongitudinal  wave  and  28,7  degrees  for  the 
quasitransverse  wave.  The  energy  flux  deviation  was  reported  to  be  about  17,8  degrees 
for  the  pure  mode  transverse  wave.  Waves  propagating  in  the  X2,  X3  plane  45  degrees 
between  the  axes  with  1^  =  0  have  the  same  solutions. 

The  solutions  for  the  wave  modes  in  an  orthotropic  medium  will  now  be  derived. 
The  matrix  components  are  now  given  by 
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Xji  —  ciilf  +  C55II  -1-  C55II, 
^22  ~  ^66^  f  “i”  ‘^22-2  ^44^3’ 
^33  “  ^55^ f  "i"  C44I2  ^33^1' 
^12  ==  ^21  ~  (‘^12  ■!■  ^66)^4’ 

^13  ~  ^31  ~  (^55  "i"  ‘^13)ll*3' 

and 

^23  =  ^32  ~  (^44  C23)l2l3- 

For  propagation  along  the  axis,  these  equations  reduce  to 

^11  “ 

^22  =  ^66^ 

^33  =  ^55^ 


and  all  other  components  are  zero.  When  these  are  substituted  into 


wave  solutions  can  be  found  to  be 


■^1 


{ cii 


2 


Po 


/ 


V2  = 


(Po  j 


and 


Vo  = 


^55 


Po 


(B.12) 


(B.13) 


equation  (2.55),  the 


(B.14) 
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with  the  first  wave  being  a  pure  mode  longitudinal  wave.  The  second  wave  is  a  pure 
mode  shear  wave  with  particle  displacements  along  the  Xg  axis  and  the  third  wave  is 
pure  mode  shear  polarized  along  the  Xj  axis.  There  is  no  energy  flux  deviation  for  these 
waves. 

For  propagation  along  X2  ,  the  components  of  Xy  reduce  to 

^11  =  ‘^66> 

^22  =  C22,  (B.15) 

^33  *=  ^44’ 


and  all  other  Xy  —  0  .  When  substituted,  these  yield  the  equation 


(C66  -  P0V^)(C22  -  P0V^)(C44  “  Po^^)  =  0 


(B.16) 


which  has  the  solutions 


Vi  = 


/  \ 


^66 


Po 


Vo  = 


-22 


Po 


(B.17) 


and 


V3 


Po 


In  this  case,  the  first  and  third  waves  turn  out  to  be  pure  mode  shear  waves  polarized 
along  the  Xj  and  X3  axes  respectively  while  the  second  wave  is  a  pure  mode  compres- 
sional  wave.  Again,  there  is  no  energy  flux  deviation  for  these  waves. 


-  136  - 


Now  the  case  of  wave  propagation  along  the  X3  axis  in  an  orthotropic  material 
is  presented.  With  =  I2  =  0,  and  I3  =  1  ^  the  components  of  Xy  reduce  to 

^11  “  ^55^ 

^22  “  ^44’  (B.18) 

^33  ^33^ 


and  all  other  Xy  components  equal  to  zero.  These  yield  the  solutions  for  the  three  modes 
of  propagation  as 


^44 


(B.19) 


and 


V3  = 


^33  2 

Po 

/ 


In  this  case,  the  first  and  second  waves  are  pure  mode  shear  waves.  The  first  wave  is 
polarized  along  the  x^  axis.  The  particle  displacements  for  the  second  wave  are  along 
the  X2  axis.  The  third  wave  is  a  pure  mode  compressional  wave.  There  is  no  energy  flux 
deviation  for  these  waves. 

As  with  the  case  of  transverse  isotropy,  the  pure  mode  waves  do  not  yield 
sufficient  data  to  determine  all  of  the  orthotropic  moduli.  Waves  must  be  propagated 
off  axis  to  determine  the  moduli  Ci2;  and  C23  .  Waves  propagating  in  the  x^,  x^  plane 
45  degrees  from  each  axis  with  13  =  0  are  considered  first.  The  direction  cosines  of  the 
wave  normal  are  given  by  ,  and  I3  =  0  .  Thus  the  Xy  matrix  components 


-  137- 


are 

^11  ~  ^66)> 

^22  =  ^(^66  +  ^22)’ 

^33  =  "^(^55  +  C44),  (B.20) 

^12  ~  ^21  =  +  ^55), 

and  all  other  Xij  =  0  .  The  solutions  for  the  waves  propagating  along  this  direction  can 
then  be  determined  to  be 

Po^l  =  ~  -f  2C56  -t-  C22  +  \/ (Cji  —  022)“^  +  4(c42  +  Cse)^! 

Pov|  =  y(c44  +  C55),  (B.21) 

and 

Po^l  =  ~  [^11  +  2C65  +  C22  —  \/(Cii  —  022)^  +  4(Ci2  +  055)^  . 

The  second  wave  is  a  pure  mode  shear  wave  while  the  first  is  quasilongitudinal  and  the 
third  is  a  quasishear  wave.  The  energy  flux  deviation  for  these  waves  is  small  because 
there  is  only  a  small  deviation  from  isotropy  in  this  plane. 

The  next  off  axis  direction  considered  is  that  in  the  x^,  X3  plane  45  degrees  from 
either  axis  where  the  direction  cosines  for  the  wave  normal  are  given  by 
li  =  I3  =  ,  and  I3  =  0  .  Therefore, 

^55)- 
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^22  “  t(‘^66  “!■  ^44)- 

2 

^33  =  ■^('^55  +  C33),  (B.22) 

Xi3  =  X31  =  ~(Ci3  4-  C55), 


and  all  other  Xjj  =  0  .  The  solutions  for  these  waves  are  then  found  by  solving  the  equa¬ 


tion 


(^22  -  PoV^)[(^ll  -  PoV^)(X33  -  PqY-)  -  Xfg]  =  0.  (B.23) 


The  three  roots  of  this  equation  are  given  by 


2C55  -1-  C33  -h  yj (cji  —  033)“  -|-  4(ci3 


PQ^2 


-44  i 


"66 


), 


(B.24) 


and 


+  2c 


55 


"33 


The  first  and  third  waves  are  quasilongitudinai  and  quasishear  mode  v/aves  respectively 
v/hile  the  second  is  pure  mode  shear.  The  energ>^  flux  deviations  are  again  large  as  was 
the  case  for  the  similar  waves  for  transverse  isotropy. 

The  final  direction  of  propagation  considered  is  that  in  the  X2,  X3  plane  45 
degrees  from  either  axis.  This  is  similar  to  the  previous  case  with  the  final  solutions 
reducing  to 


9  1 

=  - 


C22  +  2C44  +  C33  +  V  (C22  “■  C33)  -h  4(c23  -f  C44) 


PqW  “  (y  (^55  + 


(B.25) 
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and 

Pov|  =  J  [c22  +  2C44  +  C33  - 

Again,  the  second  wave  is  pure  mode  shear  while  the  first  is  quasilongitudinal  and  the 
third  is  quasishear.  These  equations  reduce  to  those  given  for  transverse  isotropy  if  the 
appropriate  substitutions  are  made.  These  equations  allow  the  determination  of  all  of 
the  second  order  moduli  for  an  orthotropic  material. 


The  relationship  between  the  measured  time  betv/een  echos  (t)  and  the  true  time 


of  Bight  (<5)  was  given  by  McSkimmin  [l9j  as 


t  =  p6 


P7 


/  ^ 
n 


360F 


(C.l) 


In  this  equation,  p  is  an  integer  expressing  the  difference  between  the  number  of  echos 
used  for  the  overlap  measurement.  Since  the  first  and  second  echos  were  used  for  these 
measurements,  it  was  always  true  that  p  =  1  .  The  variable  n  in  this  expression  is  an 
integer  v/hich  gives  the  number  of  cycles  of  incorrect  overlap  between  the  two  echos.  If 
the  echos  are  overlapped  correctly,  then  n  =  0  .  If  the  second  echo  is  overlapped  one 
cycle  before  the  first  then  n  =  — 1  and  if  the  reverse  is  true  then  n  =  1  .  7  is  the  phase 
angle  {in  degrees)  accounting  for  phase  shifts  at  the  transducer-specimen  interface.  F  is 
the  frequency  of  the  tone  burst  and  not  the  repetition  rate  frequency. 

Substituting  the  fact  that  p  =  1  simplifies  this  equation  to 

(  ^  \ 


t  =  s  — 


360F 


TL 


(C.2) 


Thus  botn  7  and  n  must  be  determined  to  evaluate  the  true  round  trip  time.  Therefore 
a  second  equation  is  introduced  v/hich  is 


At..=  ^ 
Fi 


( 


n  — 


_2l 

360 


1 

Fh 


n 


7h 

360 


(C.3) 


This  equation  expresses  the  change  in  measured  time  between  echos  as  a  function  of  a 
change  in  the  frequency  of  the  pulse  from  a  higher  frequency  (Fj,)  to  a  lower  frequency 
(F))  .  The  phase  shifts  at  these  two  frequencies  are  designated  and  yj  .  For  these 
measurements  F],  was  chosen  to  be  that  of  the  resonant  frequency  of  the  transducer  and 
F]  w’as  chosen  to  be  0,9  times  that  value.  Thus,  this  equation  can  be  rewritten  as 
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At 


1 

/  N 

7i 

1 

/ 

0.9Fr 

“  360 

Fr 

n  — 

\ 

360 

(C.4) 


where  is  the  resonant  frequency  of  the  transducer. 

The  phase  angles  7i  and  7,.  can  now  be  calculated  from  standard  electrical 
transmission  line  theory.  It  is  defined  that  Zj,  Z2;  and  Z3  are  the  mechanical  impedances 
of  the  transducer,  bond  material,  and  sample  respectively.  Z4  is  the  mechanical 
impedance  of  the  transducer-bond  composite  resonator  and  it  is  defined  that 


r 


(C.5) 


Therefore, 


Z4  —  iZ2 


r  tan  k2l2  +  tan 
r  tan  kj^lj  tan  k2l2 


(C.6) 


where  Ij,  I2  are  the  thicknesses  of  the  transducer  and  bond  respectively, 
and  k^,  k2  are  the  wave  numbers  for  the  transducer  and  bond  given  by 


ki 


OJ 

Vi 


27r 


(C.7) 


where  Vj  are  the  respective  velocities  of  sound  in  the  transducer  and  bond 
and  Xj  are  the  wavelengths  of  sound  in  the  transducer  and  bond. 

Now  it  is  assumed  that  the  bond  thickness  is  zero,  and  thus  equation  (C.6)  can  be 
reduced  to 


Z4  —  tan  kj^lj. 

Since  at  the  resonant  frequency  of  the  transducer. 


(C.8) 


(C.9) 


it  is  true  that 
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Icjii  =  180  degrees. 


(C.IO) 


At  i"',  , 


kiii  =  (0.9  •  180)  degrees. 


(C.ll) 


Therefore  at 


La  0 


{C.U 


and  at  Fj  j 


Z4  =  iZi  tan  (0.9  •  180) 


(C.13) 


Jbor  an  elastic  wave  reflecting  from  the  transducer-specimen  interface,  the  ratio 
of  the  refiected  pressure  (P,)  to  the  incident  pressure  of  the  wave  (Pj)  is  given  by 

Pr  Z..  -  Zo 


Pj  Z4  -h  Zj 


{C.14) 


and  the  phase  angle  is  given  by 


tan  7  = 


Pr  ' 

Im  : 

PJ 

Re 


"P.  ' 


(C.lc 


At  F  =  F. 


P.  0  -  Zc 


pj  0  -i-  Z3 


(C.16) 


and  therefor 


7r  =  tan  -0=0. 


(C.17) 


At  F  =  F 


where  R4  =  Zi  tan  (0.9  *  180)  . 
Thus  7]  can  be  calculated  to  be 

71  =  tan 


-1 


2Z3R^ 


p  2  nl 

(^K4  —  Z,3  j 


(C.18) 


Therefore  7^  can  be  calculated  from  knowledge  of  the  mechanical  properties  of  the  sam¬ 
ple  and  the  transducer. 

Since  7^.  =  0  ,  the  equation  for  the  change  in  measured  time  of  flight  as  a  func¬ 
tion  of  frequency  can  be  reduced  to 

(O.llln  -  7i/324) 


At(o)  = 


Fr 


(C.19) 


where  the  superscript  (0)  implies  that  this  equation  is  true  for  zero  bond  thickness. 
If  n  =  0  ,  this  equation  reduces  to 


324Fr 


(C.20) 


where  the  second  superscript  0  implies  that  n  =  0  . 

If  measurements  are  made  at  Fj.  and  for  a  number  of  different  overlap  numbers 
(n)  and  the  At's  calculated  for  each  n,  the  correct  choice  of  n  =  0  can  be  determined  by 
comparing  the  experimental  results  with  those  predicted  by  the  previous  equation.  The 
agreement  will  not  be  exact,  but  the  one  closest  to  the  calculated  value  of  and 

obeying  the  stipulation  that 

At„eas„r.<l  >  (C'21) 

which  implies  a  positive  bond  thickness  will  be  the  correct  choice  for  proper  overlap. 
This  theory  can  be  extended  to  account  for  nonzero  bond  thickness.  However,  the 
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corrections  were  determined  to  be  less  than  experimental  error  for  these  measurements. 


-145- 

APPENDIX  D 


108  REN 

PROGRAH 

TO  COLLECT  AND  ANALYZE  DATA  FRQN  PUi 

.SE  ECHO 

113  REff 
120  REN 
150  REN 
168  PEN 

170  prw 

OVERLAP  ULTRASONIC  TECHNIQUE.  WRITTEN  BY  PRI 
9/25/85. 

variable  list 

jSSER 

132  REH 

variabl 

e  description  type 

198  REH 

ii- 

naae  or  tile  character 

(32! 

pCH 

transducer  siatsrial  " 

(771 

212  m 

b$ 

bond  aaterial 

7321 

222  REfI 

5l 

saaple  ID 

(32) 

230  REN 

dl 

date  data  taken 

(281 

248  ^ 

data  pts/Qverlap  nuaeric 

252  REH 

ffi2 

series  toraat  of  at  " 

PF^ 

17 

AVE.  EXP.  TRANSIT  TIHE  ” 

PPM 

T"' 

m.  EXP.  DELTA  T(8i  " 

PPM 

P  i  T )  _ 

,F'ii  TRANSIT  TINE,  CORR.  DELTA  T 

H 

ncp  ppM 

F2 

XDUCER  RES.  FREQ.  IN  HZ  ' 

263  REH 

Ji 

iQHAXI 

-L\  DPM 

i.-J  L  l  -.WM 

F3 

DATA  FILE  NO. 

262  REN 

FILE  LENGTH  COUNTER 

DiT  DPM 

Cl 

FRES.  DIVISION  BY  HATEC  110 

S 

264  REH 

Pd), 

,Q(I)  AVS.  DELTA  T  AT  RES., LOW  FREQ.  " 

74^  pPM 

S(I) 

iR'I)  ST.  DEV.  OF  TINE  AT  RES., LOW 

FREQ.  ■ 

266  REH 

Dili 

DELTA  T'S 

n 

267  REH 

U(I1 

UNCERTAINTIES  IN  DELTA  T 

B 

268  REH 

W,X 

SUN  TOTALERS  FOR  AVS.  S 

H 

278  REH 

t 

saaple  tesip.  " 

232  REH 

tl 

Kducer  resonant  freq.  " 

290  REH 

I,K 

INTEGER  LOOP  COUNTERS  " 

300  REH 

l(i,i 

k)  data  at  low  freq.  (.9fl)  " 

318  REH 
320  REH 

h  I  i , 

k)  data  at  res.  freq. 

338  REH 
348  REH 

character  variable  diaensiDning 

350  Din  Fli3:S,):f(32i,B$(32),!)$(20),Sii32i 
360  m 

370  REH  file  naaing  and  variables 
380  RES', 

390  REH 

408  PRINT  "Mhat  is  the  naae  of  the  tile';’" 

410  INPUT  FT 

120  PRINT  "what  is  the  satple  ID?" 

430  INPUT  SI 

448  PRINT  "nhat  is  the  t.ducer  aaterial?" 

450  INPUT  Xt 

440  PRINT  "what  is  the  bond  Baterial";" 

470  INPUT  Bl 

480  PRINT  "what  was  the  date  data  was  taken?" 
498  INPUT  D$ 
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508  PRINT  "inihat  yas  the  ttop.  of  the  saaple  in  dog, 

513  INPUT  T 

520  PRINT  "whit  was  the  '^ducer  resonant  freq,  in  rlHz^" 

538  INPUT  FI 

540  PRINT  "whit  was  lONAK!^" 

551  PRINT  °=  oax,  0  of  apparent  echos  in  front  of" 

569  PRINT  "or  behind  the  reference  burst" 

570  PRINT  "(ref.  burst  is  taken  to  be  the  " 

5S0  PRINT  "larger  aoplitude  burst)" 

593  INPUT  J1 

600  PRINT  "what  is  the  I)  of  data  pts  per  overUp" 

613  INPUT  HI 

620  PRINT  "if  series  foraat  of  g!  was 
633  PRINT  "then  input  T 

640  PRINT  "if  series  foraat  of  g2  was  h2!4i8,,/ 

650  PRIr^T  "then  input  2" 

660  INPUT  H2 

673  IP  H2=l  THEN  7iS 

630  IP  H2=2  THEN  710 

699  PRINT  "expecting  input  of  1  or  2" 

700  SG  TO  629 
718  REH 

720  REN  description  of  echo-overlap  data 

733  REH  organization^  each  data  set  has  2 

740  REH  arrays  (1  for  res,  freq,  and  1  for  .9^:^rs5  freqJ 

750  REH  each  with  datapoints, 

760  REH  a  final  data  set  contains  xducer  oaterial, 

770  REH  bond  natenaij  teap,,  etc, 

798  REH  dinensioning  arrays 

800  REH 
819  LET 
823  LET  K=H1 

830  DIH  H(I,I(),P(I),Q(I),F{I) 

840  DIH  L(I,K),R(I),S(I),D(I],U(I),E(I) 

853  REH 

860  REH  input  data 

873  FOR  i(=l  TO  HI 

S89  PRINT  "INPUT  FREQ,  DIVISION  FACTOR  FROH  HATEC  Hi" 

899  INPUT  Cl 
909  FOR  1=1  TO 

919  PRINT  ”freq.  =  resonant  freq/ 

929  Z=I-(J1-H] 

929  PRINT  "PRESS  RETURN  TO" 

933  PRINT  "enter  array  GeiberHZ^K 

931  PRINT  "INPUT  YES  TO  CHANSE  FREQ.  DIVISION  FACTOR  FRGH;  ";C1 

932  INPUT  A$ 

933  IF  AJFO^YES"  THEN  943 

934  PRINT  "INPUT  m  FREQ,  DIVISION  FACTOR  " 

935  INPUT  Cl 

940  INPUT  i20;H(I,K) 

941  H(I,K)=H(I,h:)/13"Cl 

942  PRINT  H(IjK) 

943  PRINT  "  " 

950  NEH  I 
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?60  NEXT  K 

970  FOR  X-1  TO  Ml 

971  PRINT  "INPUT  FRE2.  DIVISION  FACTOR  FROM  MATED  110' 

972  INPUT  Cl 

980  FOR  1=1  TO  2»{JU1)-1 

990  PRINT  "Treq.  =  .9«res.  freq.' 

1000  I=I-(J1+1) 

1009  PRINT  "PRESS  RETURN  TO  ' 

1010  PRINT  "enter  array  ®eiiber";Z,K 

1011  PRINT  "INPUT  YES  TO  CHANGE  FREQ.  DIVISION  FACTOR  FROM: 

1012  INPUT  AJ 

1013  IF  AlO'YES’  THEN  1020 

1014  PRINT  "INPUT  NEH  FREE.  DIVISION  FACTOR" 

1015  INPUT  Cl 

1020  INPUT  @20:L(I,K) 

1021  L!I,K)=L!I,K1/10"C1 

1022  PRINT  L(!,K! 

1023  PRINT  '  ■ 

1030  NEH  I 
1040  NEXT  K 

1050  REM  COMPUTE  AVG.'S  AND  ST.  DEV.'S 
1060  REM 

1070  FOR.  1=1  TO  2MJl+l)-l 
1080  W=8 
1090  X=0 

1100  FOR  K=l  TO  Ml 
1110  IF  M2<>1  THEN  1140 
1120  W=H+K/H(I,K) 

1130  BO  TO  1160 
1140  Z=2'-(K-1! 

1150  H=H+I/H':i,K! 

1160  NEXT  K 
1170  P(I)=W/M1 
1180  FOR  K=1  TO  Ml 
1190  IF  M2')l  THEN  1220 
1200  X=X+(K/H(I,X)-P(I)!‘-2 
1210  60  TO  1240 
1220  2=2'- (K-1! 

1230  X=X+(2/H(I,t:;!-P(I)l-2 
1240  NEXT  K 
1250  S(I)  =  iX/Mir'0.5 
1260  NEXT  I 

1270  FOR  1=1  TO  2MJ1+1)-1 
1280  W=0 
1290  X=0 

1300  FOR  K=1  TO  Ml 
1310  IF  M2<>1  THEN  1340 
1320  «=H+K./L(I,K) 

1330  SO  TO  1360 
1340  2=2'- (K-1 1 
1350  «=H+2,'LiI,Kl 
1360  NEXT  K 
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n/yi-M/Ht 

1380  FOR  K=1  TO  HI 
1391  IF  ri2<)i  THEN  1428 
i4ii  ^^HK/L(I,K)-Q(I))"2 
I4ie  GO  TQ  1440 
1428  Z=2"(K'I) 

1430  ^=)(r(Z/L(I,K;-S(n)"2 

1448  Htll  K 

‘453  - ^ I ) = i K/Hl ) "i, 5 

1460  NEaT  I 
1470  REH 

Mee  REH  calculhteb  delta  t  and  uncertainties  in  delta  T 

1490  REH 
1300  REH 

1510  FOR  1=1  TO  2v(JU15-i 
1320  D(I)=0(I)-R{I) 

1330  LI(n=S(Ii^R(I) 

1348  NEm  ! 

1351  REH 

1568  REH  CALCULATES  DELTA  T(ei 


1370 

REH 

1580 

PRINT 

"NHAT 

IS  THE 

IHPEDhNCE 

'LHE5 

HECH, 

OHRS!  Cf" 

1398 

PRINT 

"THE  " 

,  Y*.  i> 

j  -r  j 

TRANSDUCERS 

n 

1608 

INPUT 

Z1 

1618 

PRIN*^ 

"NHhT 

IS  THE 

IHPESANCE 

(ME5 

HECH. 

OHRSl  OF  SARPLE  "sS'l 

1620 

INPUT 

L  ’J 

1630 

PRINT 

"NHAT 

IS  THE 

IHPEDANCE 

(^1E3 

HECH. 

OHRS)  OF  THE  BOND 

1648 

INPUT 

7  7 

1658 

Z4=:U 

‘TANi0. 

9^FI) 

S=ATN^ 

l2iUn 

3/(14" 

7_77''-  ‘ 

i.  L  .■  ; 

i  i.  i,  i 

i  U  i 

IF  S>i 

I  THEN 

1671 

1662  6=8+PI 

1670  F2=FlMO0ii00 

1680  Cl=-*l/(2^PUF2)^(G/e,9) 

1690  REH 

1780  REH  IDEHTIFV  CORRECT  N  AND  FIND  DELTA  T(i)  FROH  UP.  DATA 
1710  REH 
1^20  REH 

1730  PRINT  'TOHPARE  THE  CALCULATED  DELTA  T(i)  =  ^;01r  (SECJ^ 
1740  PRINT  ^UITH  THE  FOLLOWING  VALUES^ 

1750  PRINT 
1760  FOR  1=1  TO 

17:^0  PRINT  "DELTA  T  ( I-(JIH) ^ )  =  ^0(1)5“  iSECJ^ 

1788  PRINT 
r9i  NEXT  I 

iseo  PRINT  "ENTER  THE  #  OF  DELTA  T  THAT  HOST  CLOSELY  HATCHES  THE" 
1810  PRINT  "CALCULATED  DELTA  T(i)^ 

1820  INPUT  B 
1830  0=8 
IG48  REH 

135S  REH  CALCULATE  EKPERIHEfiTAL  DELTA  T(8)  AMD  AVG,  TRANSIT  TIRE 
1S6B  REH  AND  THEIR  UNCERTAINTIES  USING  A  HEISHTED  AV6. 

1872  RER 
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1830  n^0 
1390  X2=0 
1900  Y1=0 
1910  Y2=0 

1920  FOR  M  TO  2*(J1+1!-1 
1938  2=I-(J1+1) 

1931  N=Z-E 

1940  E(I)=D(I!+N/F2t(l/0.9-l) 

1950  F(!)=P(!)+N/F2 
I960  ;<l=n+E(I)*(l/UlI)“2) 

1978  Yl^il+FiDMl/S'D-'Z) 

1930  X2=l(2+l/y(I!''2 
1990  Y2=Y2+1/S(I)''2 
2008  NEXT  I 
2010  T2=:X1/X2 
2020  T3=Y1/Y2 
2030  Hl=X2'--0.5 
2040  H2=Y2''-0.5 
2080  60SUB  3430 
2090  REN  STORE  DATA 
2100  REtl 
2110  TLIST 
2120  PRINT  ■■ 

2130  PRINT  "ENTER  «  OF  (LAST)  FILE  TO  STORE  DATA  ON" 
2140  INPUT  F3 
2150  Z=30 

2160  Z=Z+LEN i S4 ) +2+LEN(Dti +2+LEN (Bt 1 +2+LEN (Ft! +2 

2170  Z=Z+LEN(X$l+2 

2180  Z=Z+10*(9+2K2*(Jl+l)tHl)) 

2190  2=Z+10t(6»2f(Jl+l)l 
2208  FIND  F3 
2210  MARK  1,Z 
2220  FIND  F3 

2230  WRITE  Ff,St,Bt,Dt,Xt,Jl 

2240  WRITE  T-,F1,F2,M1,H2,Z1,Z2 

2250  WRITE  Z3,Z4,S,Q1,D,C6,C7 

2260  FOR  1=1  TO  2t(Jl+l)-l 

2270  WRITE  P(I),Q(I1,S(I),R(I1,D(!),U(II 

2230  NEXT  I 

2290  FOR  1=1  TO  2t(Jl+ll-l 
2300  FOR  K=1  TO  Ml 
2310  WRITE  H(I,Ki,L(I,K! 

2320  NEXT  K 
2330  NEXT  I 

2340  WRITE  T2,T3,H1,H2 
2350  CLOSE 

2360  PRINT  "DATA  STORED" 

2370  PRINT  "' 

2410  PRINT  @2: "FILE  NAME  IS  ";F$ 

2420  PRINT  §2: "FILE  I  IS  ';F3 

2430  PRINT  02: "TRANSDUCER  MATERIAL  IS  *:X$ 

2440  PRINT  02: "BOND  MATERIAL  IS  ";Bt 
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2458  PRINT  i2:“'SAflPL£  ID  IS  'iSS 

2460  PRINT  t2:“DATE  BATA  MAS  TAKEN  IS  ";DS 

2461  PRINT  12:  “SAMPLE  TEflP  WAS  ";T;"  C" 

2478  PRINT  02: "DATA  PTS/GVERLAP  IS  ";fl! 

2480  IF  fl2<:>l  THEN  2528 

2493  PRINT  12; "SERIES  FGRHAT  OF  HI  VALUES  IS;' 
2530  PRINT  @2:'l:2-.3;4..." 

2518  6Q  TO  2548 

2523  PRINT  §2: "SERIES  FQRHAT  OF  HI  VALUES  IS;" 
2530  PRINT  92:”1:2:4;S...’'‘ 

2543  PRINT  12:"  ” 

2553  FOR  Z^l  TG  2 

2563  IF  201  THEN  259S 

2570  PRINT  12; “DATA  AT  RESONANT  FRE8." 

2583  30  TO  2610 
2593  PRINT  §2:“  ' 

2683  PRINT  §2;  "DATA  AT  .?-}RESONANT  FRESO 
2618  PRINT  S2:"  " 

2628  IHA6E  3A,17,S 
2638  IRASE  2D,3X,S 
2648  PRINT  §2;  L'SINS  2628;  "Q/R" 

2658  FOR  1=1  TO  Hl-l 
2668  IF  R202  THEN  2698 
2678  tl=2'"(I-Il 
2688  SO  TO  2781 
2698  W=I 

2788  PRINT  82;  USING  2638; h 

2718  NEXT  I 

2728  IMA6E  2D 

2730  IF  H2<>2  THEN  2768 

2748  li=20Rl-l) 

2758  60  TO  2771 
2760  W=H1 

2778  PRINT  12:  USING  2723:8 

2780  FOR  I=I  TO  20J10S-1 

2798  IHA6E  2D, IK, 3 

2888  PRINT  82:  USING  2790: I-'JI+ll 

2318  FOR;  K=1  TO  Hl-l 

2828  IP  ZOl  THEN  2868 

2830  IRAGE  7D.2D,1X,S 

2843  PRINT  §2;  USING  283e:H':i,Ki 

2858  GO  TO  2870 

2868  PRINT  §2:  USING  2838;L(!,K;; 

2878  NEXT  K 

2888  IP  ZOl  THEN  2920 

2898  PRINT  82;  USING  2988:8(1, HI! 

2988  IRAGE  7D.2D 
2918  GO  TO  2938 

2923  PRINT  92:  USING  2980; L( I, HI) 

2933  NEXT  I 
2948  NEXT  2 
2958  PRINT  12;"  " 
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2979  IF  Z=2  THEN  3910 

2989  PRINT  §2: 'DATA  AT  RESONANT  FREQ.' 

2999  60  TO  3010 

3000  PRINT  82: 'DATA  AT  .9«R£S0NAHT  FREQ.’ 

3010  PRINT  82:'  ' 

3020  PRINT  82:  USIN6  3030: 'Q' 

3030  IHASE  2A,1)(,S 

3040  PRINT  82:  USIN6  3050: 'AV6.  TINE  (SEC.)' 

3050  IHASE  28A,IX,S 

3060  PRINT  82:  USIN5  3079; 'STD.  DEV.  (SEC.)" 

3070  1HA6E  20A 
3090  PRINT  82:'  ' 

3090  FOR  1=1  TO  2t(Jl+l)-l 

3100  IF  1=2  THEN  3180 

3110  PRINT  82:  USING  3120: I-(J1+1) 

3120  IHAGE  2D,1X,S 

3130  PRINT  82:  USING  3140:P(I) 

3140  IHAGE  10E,10X,S 

3150  PRINT  82:  USING  3160:3(1) 

3160  IHASE  18E 
3170  GO  TO  3210 

3180  PRINT  82:  USING  3120: I-(J1+1) 

3190  PRINT  82:  USING  3140:0(1) 

3200  PRINT  82:  USING  3160:R(II 
3210  PRINT  82:"  ’ 

3220  NEXT  I 
3230  PRINT  82:"  * 

3240  NEXT  Z 
3250  PRINT  82:'  ' 

3260  PRINT  82:'IHPEDANCE  OF  XDUCER  IS  ';Z1;"  (tlE5  HECH.  OHHS)' 

3270  PRINT  92:"IHPEDANCE  OF  SAHPLE  IS  ';Z3;'  (HE5  HECH.  OHHS)' 

3280  PRINT  82:'IHPEDANCE  OF  BOND  IS  ';Z2;'  (HE5  HECH.  OHHS)' 

3290  PRINT  82: "PHASE  SHIFT  AT  LON  FREQ.  IS  ■;G»360/(2tPI);'  (DEG.)' 
3300  PRINT  82: "CALCULATED  DELTA  T(0)  IS  ';01;'  (SEC.)' 

3310  FOR  1=1  TO  2t(Jl+l)-l 
3320  PRINT  82:'  ■' 

3330  PRINT  82:'DELTA  T(";I-(J1+1);')  IS  'jDdl;'  (SEC.)' 

3340  PRINT  82: 'UNCERTAINTY  ("|I-(Jl+l)i')  IS  "lUIDi'ISEC.)' 

3350  NEXT  I 
3360  PRINT  82:'  ' 

3370  PRINT  82: 'CHOICE  FOR  CORRECT  OVERLAP  IS  'jB 

3388  PRINT  82:'AVS.  EXPERIHENTAL  DELTA  T(N=0)  IS  'iT2;'  (SEC.)' 

3390  PRINT  82: 'ST.  DEV.  IS  ";H1;'  (SEC.)' 

3400  PRINT  82:*  ’ 

3410  PRINT  82:'AV6.  TRANSIT  TIHE  FROH  EXP.  DATA  IS  'iTS;'  (SEC.)' 
3428  PRINT  82: 'ST.  DEV.  IS  ■;H2;'  (SEC.)' 

3421  PRINT  82:'  ' 

3422  PRINT  82: 'BOND  ANGLE  IS  ';C6i'  DEGREES" 

3423  PRINT  82: "TIHE  CORRECTION  FOR  BONO  IS  ';C7;'  SEC.' 

3424  GO  TO  3740 
3430  REH 

3440  REH  BOND  CORRECTION 
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345fl  REH 

3460  Din  V(38),C(3@),Q(3S) 

3479  FOR  1=1  TO  3Q 

3480  Cl={I-l)/2-J(2'5FII/36S 

3490  Z5=Z2-HZ2/2!-JTAIl(Cl)^-TAW(8.95PI)l/(Z2/Zl-TAn(S.95PI)?TflN(Cr) 

3510  C2=-2-iZ52Z3 

3513  C3=Z3sZ3-Z55Z5 

3523  C4=-23Z25Tfin(Cl)3Z3 

3538  C5=Z3523-Z2ifTAfI(Cl)5(Z25Tftn(Ci)) 

3540  V'I)=ATfl(C4/C51 

3541  V(I)=V(!  15360/ (2-iPl! 

3553  C(I)=ATn(C2/C3) 

3551  IF  C(I)>S  THEN  3555 

3552  C(I)=C';i!5P! 

3555  C(I!=CiI)f:-360/(25Pi> 

3560  C !  1 1  =y  1 1  /  (3685F2)  -C !  1 1  /  (36050. 95F2) 

3573  NEXT  I 
3580  PAGE 

358!  VIEHPORT  13,122,10,90 

3590  HINDOH  3,15,0(115550(11 /li, 3 

3591  nO'.’E  2,8 

3592  FOR  1=2  TO  16  STEP  2 

3593  ROVE  1*2,3 

3594  PRINT  “l";I-2i 

3595  NEXT  I 

3680  AXIS  2,O';ll/ie,0,e 
3618  flOVE  e,Qtl! 

3620  FOR  1=1  TO  38 
3633  DRAH  (i-l)/2,G(I; 

3640  NEXT  ! 

3650  ROVE  e,T2 

3660  DRAN  30, T2 

3661  HOVE  g,T25Hi 

3662  DRAH  15,T2+H1 

3663  MOVE  0.T2-H! 

3664  DRAi?  15,T:-H1 

3665  ROVE  §,B5O'1^'10 

3673  PRINT  "INPuT  PHASE  AnGlE  Q*  ;4TE-SE:':2V  " 

36S3  INPUT  Cs 
36  S:  HONE 
3682  PAGE 
36R0  I=25C6-1 

r7rU  ‘  ’  ;  '  7^  ^  V  ! 

37i^  pen.-  DESREES  ^ 

3720  ::]r:;E;TiGN  ^QR  BOND  IS  SEC/ 
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100  REH  PROBRAH  TO  MEASURE  STRESS  VS.  STRAIN 

120  REH  HRITTEN  BY  PROSSER  7/15/86 

130  REH 

HB  REH  VARIABLES 

NAME 

TYPE 

150  REH 

160  REH  SAMPLE  ID. 

SI 

CHARACTER(32) 

170  REH  DIRECTION  OF  STRAIN  1 

HI 

*  (12) 

180  REM  DIRECTION  OF  STRAIN  2 

Nl 

“  (12) 

200  REH  ■  "  LOADING 

LI 

■  (12) 

220  REH  DATE 

01 

CHARACTER (12 

270  REH  CROSS  SECTIONAL  AREA  OF  SAMPLE  (SQ.  IN.lAl 

• 

280  REH  EXCITATION  VOLTAGE  FOR  STRAIN  GAGE  1 

A2 

290  REH  ■  '  •  *  ”2 

A3 

300  REH  NUMBER  OF  DATA  POINTS 

Bl 

1 

330  REH  VOLTAGE  HEASUREHENTS  FOR  LOAD 

F(I) 

ARRAY 

340  REH  VOLTAGE  HEASUREHENTS  FOR  STRAIN  1 

V(I,11 

a 

350  REH  "  ■  ■  STRAIN  2 

V(I,2! 

■ 

400  REH  STRAIN  1 

L(I) 

ARRAY 

410  REH  STRAIN  2 

S(l) 

1 

460  REH  STRESS 

Ddl 

1 

520  REH  MIN  STRESS 

Nl 

■ 

530  REM  MAX  STRESS 

U2 

fl 

540  REH  HIN  STRAIN  1 

H3 

R 

550  REM  MAX  STRAIN  1 

N4 

1 

560  REM  HIN  STRAIN  2 

NS 

1 

570  REH  MAX  STRAIN  2 

U6 

R 

571  REH  SUH  Ddl 

X3 

B 

572  REM  SUH  iD(I>)^2 

X4 

R 

573  REM  SUH  L(l) 

X5 

f 

574  REH  SUH  (L(I))^2 

X6 

R 

575  REM  SUM  S(I) 

X7 

B 

576  REH  SUH  (S(I))'‘-2 

X8 

V 

577  REH  SUH  (D(I)»l(i) 

H7 

B 

578  REH  SUH  (D(I))»s(i)) 

m 

■ 

590  REM  DUMMY  VARIABLE 

ct 

CHARACTER (3: 

600  REH 

A$ 

i 

630  REH  DUMMY  VARIABLE 

A5 

R 

640  REH  INTEGER  COUNTERS 

I,J 

H 

650  REM  DUHHY  ARRAY  FOR  PLOTTING 

Pd) 

ARRAY 

660  REH  "  VARIABLES  FOR  PLOTTING 

A8,A9 

NUMERIC 

670  REH  PLOTTER  1 

B2 

B 

6S0  REM  REGRESSION  COEF.  FOR  STRESS-STRAIN  1 

R1 

B 

681  REH  "  "  ■  STRESS-STRAIN  2 

R2 

690  REH  SLOPE  FOR  STRESS-STRAIN  1 

Cl 

B 

691  REM  ■  ’  STRESS-STRAIN  2 

C2 

B 

692  REH  Y-INTERCEPT  FOR  STRESS-STRAIN  1 

El 

B 

693  REH  “  "  STRESS-STRAIN  2 

E2 

R 

700  REfl 

710  REH  DIMENSION  CHARACTER  VARIABLES 
720  REH 


-154- 


730  Dili  51(32), flS!12),f«(12),LI(i2),D»(!2) 

740  DIH  C$(321,A$i32) 

753  REr] 

763  REH  INPUT  INITIAL  VARIABLES 
773  REH 

783  PRINT  "INPUT  SAilPLE  ID." 

79S  INPUT  53 

830  PRINT  "INPUT  DIRECTION  OF  LQADINB" 

810  INPUT  LI 

320  PRINT  "INPUT  DIRECTION  OF  STRAIN  I’ 

333  INPUT  HI 

840  PRINT  "INPUT  OIRECTION  OF  STRAIN  2" 

858  INPUT  Nl 

860  PRINT  "INPUT  EXCITATION  VOLTAGE  FOR  STRAIN  BA6E  1 
383  INPUT  A2 

398  PRINT  "INPUT  EXCITATION  VOLTAGE  FOR  STRAIN  SAGE  2 

391  INPUT  A3 

908  PRINT  "INPUT  DATE" 

918  INPUT  Dl 

998  PRINT  "INPUT  CROSS  SECTIONAL  AREA  OF  SAMPLE  (SS.  IN.)" 

1308  INPUT  A1 
1300  RER 

1310  REH  DIHENSION  VARIABLES  HAXIflUH  230  DATA  PTS.  TAKEN  ABOUT  EVERY 
1323  REH  ONE  SECOND 

1333  DIH  F(230),V(283,2!,LI28e),S(28e),Di2Sg) 

1353  REH 

13i3  REH  SET  UP  INSTRQfl  AND  TAKE  DATA 
1370  REH 

1388  PRINT  “PRESS  RETURN  TO  BE6IH  TAKING  DATA" 

1390  PRINT  "PRESS  BREAK  TO  STOP  TAKING  DATA" 

1391  PRINT  "ENTER  RUN  1510  TO  RESUHE  PR08RAH" 

1438  INPUT  Al 

1413  FOR  1^1  TO  230 
1423  INPUT  a24:F(I) 

1430  INPUT  312:911,1) 

1443  INPUT  323:9! I, 2) 

1453  Bl=! 

1473  FOR  J=l  TO  130 
1483  B?=J 
1493  NEU  J 
1538  NSn  I 

1513  REH  CALCULATE  STRESSES  AND  STRAINS 
1528  RER 

1540  FOR  1=1  TO  Bl 

1553  D!I!  =  (F(I)-F!D)/Al-7li3a-i6S95 

1 563  L  ( I ) =-4-i (V  ( 1 , 1 )  -V !  1 , ! )  i  / ! A2i2 . 075) 

1570  S ( I) =-4i ! V ( 1 , 2) -V ! 1 , 2) ) / (A3«2. 375) 

I5SS  NEXT  ! 

1890  REH 

1903  REH  CALCULATE  HAX  AND  HIM  STRESS,  STRAINS 
1918  ai=! 
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1923  W2=-l 

1930  il3=l 

1940  H4=-l 

1950  «5=1 

I960  «6=-l 

1970  FOR  1=1  TQ  B1 

1980  IF  L(n)W3  THEN  2000 

1990  M3=L(!) 

2000  IF  L(I)<H4  THEN  2020 
2010  N4=L(I) 

2020  IF  D(I)>H1  THEN  2040 
2030  m=D(I! 

2040  IF  D!II<N2  THEN  2060 
2050  N2=D(I) 

2060  IF  S(I)>H5  THEN  2080 
2070  H5=Sa) 

2080  IF  S(I)<H6  THEN  2100 
2090  H6=S(I) 

2100  NEH  I 
2110  REN 

2120  REN  LINEAR  RE6RESSI0N 

2130  REN 

2140  SOSUB  7000 

2150  REN  PLOT  STRESS  AND  STRAIN 

2160  REN 

2170  SOSUB  2190 

2180  GO  TQ  3200 

2190  B2=32 

2200  FOR  1=1  TO  2 

2210  PAGE 

2211  PRINT  §B2: ‘DATE:  '.jDt;’  SAMPLE:  'iSI;"  DIR.  OF  STRESS:  °-,li 

2212  PRINT  082: "DIR.  OF  STRAIN  1:  ";N$|"  DIR.  OF  STRAIN  2:  ”;NJ 
2220  VIEWPORT  20,120,30,90 

2230  DIN  P(8) 

2240  P(1!=W1 

2250  P(2)=«2 

2260  P(3)=5 

2270  IF  1=1  THEN  2310 

2283  P(5)=M5 

2290  P(6)=H6 

2300  SO  TO  2330 

2310  P(5)=H3 

2320  P(6)=W4 

2330  P(7)=5 

2340  P5=3 

2350  SOSUB  5008 

2360  P5=7 

2370  SOSUB  5000 

2380  WINDOW  P(1),P(2),P(5),P(6) 

2390  AXIS  SB2:P(31,P(7),0,B 
2400  P5=4 
2410  A$=“HHH" 
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2428  60SUB  ii8§ 

2438  P5=:8 
2448  Ai>='’=’ 

2451  SOSUB  6880 

246S  !F  1=2  THEN  2570 

2478  HOVE  §B2;Dii) ,1(1) 

2480  FOR  .J=i  TO  Bl 

2498  DRAW  iB2:D(J) ,L(J) 

255i  NEIT  j 

2568  80  TO  266i 

2570  HOVE  §B2:D(l),S!li 

2588  FOR  J=1  TO  SI 

2598  CRAM  §B2:D(J) ,S!J) 

2650  NEH  J 

2660  IF  1=2  THEN  2780 

2678  A8=03 

2680  A9=94 

2698  eo  TO  2728 

2783  A8=t)5 

27  li  A9=116 

2720  HOVE  §B2:iil,A? 

2733  PRINT  8B2: ’HHHHHJJ"; 
2740  00  TO  2813 
2758  FOR  J=1  TO  L£N(C$! 
2763  A5=SEG!C0,J,1) 

2773  PRINT  iB2;A-4i ’HJ"; 
2788  NEU  J 
2798  HOVE  §B2:Sl,A8 
2683  PRINT  iB2;=JJJ 
2381  SO  IQ  2842 
2810  IF  1=2  THEN  2841 
2320  Of='STRAIN  2" 

2830  80  TO  2758 
2340  i:$=“STRAIN  T 
2841  30  TO  2750 
2342  CJ=”STRESS  (Pa)^ 

2843  PRINT  @B2;C« 

2853  IF  82032  THEN  2386 
2860  Z9=3G 
2378  60  TO  2890 
2330  Z9=2 
HQHE 


2*^09 

FOR  ze=l  TC  Z9 

2"li 

PRINT  iB2f^ 

292i 

NEIT  za 

2939 

PRINT  STRESS? 

Pa 

HAL 

STRESS? 

2960 

IF  1=2  THEN  2990 

2970 

PRINT  STRAIN  li 

HAL 

STRAIN  I?  LeJ4 

2980 

8Q  TO  3020 

299Q 

PRINT  STRAIN  2? 

HAL 

STRAIN  2?  Ltl6 

3020 

IF  1=2  THEN  3050 

•jSo-jO 

PRINT  OeOrSLQPE  IS?  ^;Ci; 

(Pa-1 

l  0 

-157- 


3040  60  TO  3080 

3050  PRINT  9B2:'SL0PE  IS:  ■;C2;'  (Pa''-1)" 
3080  IF  1=2  THEN  3110 
3090  PRINT  @B2:'Y-INTERCEPT:  ";E1 
3100  60  TO  3130 

3110  PRINT  «B2:'Y-INTERCEPT:  ";£2 
3130  IF  1=2  THEN  3160 
3140  PRINT  0B2:*R£BRESSION  CQEF.:  'iRl 
3150  60  TO  3170 

3160  PRINT  0B2; ’REGRESSION  COEF.:  'iRO 
3170  INPUT  A$ 

3180  NEXT  I 
3190  RETURN 

3200  HONE 

3201  PAGE 

3210  PRINT  ’HARD  COPY  OF  PLOTS?  YES=1  N0=2'’ 

3220  INPUT  A5 

3230  IF  A5<>1  THEN  3270 

3240  PRINT  "INPUT  PLOTTER  NUHBER* 

3241  INPUT  B2 
3250  6QSUB  2200 
3260  60  TO  3300 

3270  IF  A5=2  THEN  3380 

3280  PRINT  "BAD  CHOICE,  TRY  AGAIN" 

3290  50  TO  3210 

3300  PRINT  ’STORE  DATA?  YES  =1  N0=2’ 

3310  INPUT  A5 

3320  IF  A501  THEN  3350 

3330  60SUB  3000 

3340  GO  TO  3300 

3350  IF  A5=2  THEN  3380 

3360  PRINT  ’BAD  CHOICE,  TRY  AGAIN" 

3370  60  TO  3300  , 

3380  REN 

3460  PRINT  ’ANOTHER  RUN?  Y£S=1  NQ=2’ 

3470  INPUT  A5 

3480  IF  A5<)1  THEN  3500 

3490  60  TO  1370 

3500  IF  A5=2  THEN  3530 

3510  PRINT  "BAD  CHOICE,  TRY  AGAIN" 

3520  GO  TO  3460 
3530  END 

5000  REN  SUBROUTINE  FOR  PLOTTING 

5001  REN  P(P5)=N1NINUN  t  OF  TICS 
5010  Pl=(P(?5-l)-P(P5-2!)/P(P5! 

5020  P2=10^INT(LST(P1)) 

5030  P1=P1/P2 

5040  IF  PI >2  THEN  5080 

5050  IF  Pl=l  THEN  5120 

5060  P2=2*P2 

5070  60  TO  5120 

5080  IF  Pl>5  THEN  5110 
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5093  P2=55P2 
5103  60  TO  5120 
5113  P2^I35P2 

5120  REfl  ADJUST  DATA  imilHUf! 

5138  Pl  =  IllT(P(P5-2)/P2) 

5143  P3=P25(P!*2) 

5153  IF  P3(P!P5-2!  THEN  518Q 

5163  P3=P3-P2 

5178  BO  TO  515B 

5183  PfP5-2)=P3 

5tP3  R£H  ADJUST  BATA  HAAiHUrl 

5229  Pl=i;iT(P(P5-l)/P2) 

5210  P3=P2iiPl-2) 

522S  IF  P(P5-1)<P3  THEN  5253 
5238  P3=P3+P2 
5243  GO  TO  5223 

5259  P(P5-1!=P3 

5260  REH  P(P5!=ADJ!JSTED  TIC  INTERVAL 
5279  P{P5;=P2 

52BS  RETURN 
6060  REH  LABLE  AHS 
6818  P4=P!P5-i) 

6829  P':4:'=P(i) 

6338  P'Si^PlS) 

6040  ?3=ABS{P(P5-3)+P4)  RAI  ABS(P(P5-2>-P4) 

6859  P3=INT!LST(P3)-M.SE-8) 

6869  P2=I3'-P3 
6879  Pl=P!P5-2)-P4/2 
6889  ?!P5)=P1P5)+P4 
6893  IF  P';P5)>PI  THEN  6149 
6139  HOVE  iB2’.P!4),Pi8) 

6119  PRINT  AS; 

6129  PRINT  SB2;  USING  ''-D.2D,S“;P{P5)??2 
6138  60  TO  6089 

6149  IF  P3=0  THEN  6138 

6150  P(P5i=P2 

6160  HOVE  eB2:P{45,P(S) 

6178  PRINT  §B2:  USING  “2A,+FD,S":=  E";P3 
6  ISO  RETURN 
7000  REH 

7019  REM  SUBROUTINE  TO  DO  LINEAR  REGRESSION  ON  BATA 

7829  REH 

7838  73=0 

7340  74=0 

7050  75=0 

7863  76=0 

7670  77=9 

7888  78=0 

7890  i17=3 

7138  tis=g 

7119  FOR  1=1  TO  B1 
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7120  n=](3+D{I) 

7130  X4=K4+DiI)-'2 
7140  X5=X5+L(I) 

7150  Xi=X6+LiI)-'2 
71i0  X7=X7+S(I) 

7170  X8=Xa+S(I)'‘'2 
7180  W7=«7+D(I)*L(I) 

7190  iil8=«a+D([)tS(Il 
7200  NEXT  I 

7210  C1=(B1*«7-X3*X5)/{B1»X4-X3'2! 

7220  C2=  (BHt«3-X3*X7) / (BUX4-X3‘'21 
7230  El=!X4tX5-X3»«7!/(BltX4-X3-'2) 

7240  E2=(X4*X7-X3tH8)/(Bl*X4-X3'‘'2) 

7250  Rl=iBl»«7-X3tX5)/((BUX4-X3"2)*!BltX6-X5'‘-2!)'’'0.5  • 

7260  R2=(Bl*W9-X3*X7),'(':BltX4-X3‘2)*(BUX8-X7‘2))*0.5 
7270  RETURN 
8000  REM 

8010  REM  SUBROUTINE  TO  STORE  DATA 

8020  PRINT  “INPUT  TITLE  FOR  STESS-STRAIN  1  FILE  (MAX  64  LETTERS)’ 
8830  DIM  8$(64),K$(64) 

8040  INPUT  B$ 

8050  PRINT  'INPUT  TITLE  FOR  STRESS-STRAIN  2  FILE  (MAX  64  LETTERS)’ 
3060  INPUT  KS 

8070  A5=0 

8071  A5=A5+LEN(B4)+LEN(K$)+10«(B1»21+10 
3030  TLIST 

8090  PRINT  ■” 

3100  PRINT  ’INPUT  #  OF  (LAST)  FILE  TO  STORE  DATA  FOR  STESS  STRAIN  1’ 

9110  INPUT  M3 

3120  FIND  M3 

3130  MARK  1,A5 

8140  FIND  M3 

3150  NRITE  B$ 

8160  WRITE  B1 
3170  FOR  1=1  TO  B1 
8180  WRITE  D(I),L(I) 

3190  NEXT  I 
8200  CLOSE 
8210  FIND  M3+1 
8220  MARK  1,A5 
6230  FIND  M3H 
3240  WRITE  K$ 

3250  WRITE  B1 
9260  FOR  1=1  TO  B1 
3270  WRITE  D(I),S(I) 

3280  NEXT  I 
3290  CLOSE 

3300  PRINT  “DATA  FOR  STRESS-STRAIN  I  STORED  ON  FILE  ":M3 
3310  PRINT  "DATA  FOR  STRESS-STRAIN  2  STORED  ON  FILE  ";M3+1 
8320  RETURN 


In  this  appendix,  the  derivations  of  the  independent  third  order  elastic  moduli  for 
the  symmetries  of  transverse  isotropy  and  orthotropy  are  considered.  After  considera¬ 
tions  of  the  symmetry  of  the  stress  and  strain  tensors  and  the  existence  of  a  strain 
energy  function,  the  number  of  independent  third  order  moduli  for  the  most  general 
anisotropic  solid  is  reduced  to  56.  The  number  of  independent  third  order  elastic 
coefBcients  can  be  reduced  further  if  the  effect  of  the  elastic  symmetry  of  the  material  is 
considered.  As  with  the  case  of  the  linear  elastic  moduli,  the  effect  of  symmetry  rota¬ 
tions  on  the  strain  energy  function  v/hich  is  extended  to  include  terms  up  to  third  order 
can  be  used  co  determine  the  independent  moduli  for  the  different  symmetries.  Since 
only  third  order  terms  are  being  considered,  only  the  portion  of  the  strain  energy  with 
third  order  terms  (cPs)  need  be  considered.  It  is  given  by 

1  1 

*^3  T Yj'^AAAVa  "T  2  CaaB^A%  +  S  ^ABDVa^'%-  (^-1) 

A  ^  A<B<D 

The  process  is  exactly  the  same  as  was  the  case  of  reduction  of  the  second  order 
moduli.  The  strain  energj^  function  is  written  out.  Then  the  effect  of  a  specified  rota¬ 
tion  on  the  strain  tensor  components  is  determined  and  then  substituted  into  the  strain 
energy  function.  Then  the  rotated  and  original  strain  energy  functions  are  equated 
term  by  term  to  determine  the  relationships  between  the  moduli.  However,  it  would  be 
both  too  time  consuming  and  space  consuming  to  v/rite  out  the  necessary  equations  for 
each  symmetry  rotation.  Thus,  in  this  appendix  only  the  effects  on  the  moduli  for  each 
rotation  will  be  given.  These  are  determined  by  simple  algebraic  manipulations. 

Tne  orthotropic  moduli  are  determined  first  by  considering  the  effects  of  two  fold 
rotations  about  each  of  the  three  axes.  For  a  180  degree  rotation  about  the  X2  axis 
which  has  a  transformation  matrix  given  by 
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-10  o' 

[aj  =010,  (F.2) 

[  0  0-1 

24  of  the  moduli  drop  out  leaving  32  independent  moduli.  The  moduli  which  drop  out 
are  Cn4,  Cue,  Ci24,  CjQe,  Ci34)  Ci38>  <^146,  Ci66t  C224>  ^2281  <^234,  CQse,  C24S,  C266,  C334,  0335,  0345,  C356, 
C444,  C448,  C466,  C496,  C566,  and  Cgea  .  These  can  all  be  shown  to  equal  zero. 

The  effect  of  a  two  fold  rotation  about  the  Xj  axis  is  to  reduce  the  number  of 
independent  moduli  to  twenty.  The  moduli  C666,  C115,  C22S)  <^336)  *^448;  <^125;  ^issi 

C236,  C34e,  and  Csge  all  become  equal  to  zero.  This  leaves  the  remaining  independent 
moduli  as  Cm,  C1121  Cus,  C122,  C123,  C133,  C144,  Cjss,  Cijj,  C222)  C223>  ^233)  <^244»  ^265<  ^268)  C333, 
C344,  C3S6,  C3ee,  and  C45e  .  As  with  the  case  of  the  second  order  moduli,  the  third  sym¬ 
metry  rotation  about  the  X3  axis  does  not  reduce  the  number  of  independent  moduli 
further.  Thus  the  20  independent  moduli  for  an  orthotropic  material  are  those  listed 
above  with  all  other  c^bd  =  0  . 

Again,  the  transversely  isotropic  moduli  can  be  determined  by  further  reducing 
the  orthotropic  moduli  by  applying  the  additional  symmetry  condition  that  the  material 
be  isotropic  in  the  Xj,  X2  plane.  This  means  that  the  material  can  be  rotated  by  any 
angle  around  the  X3  axis  and  be  in  an  elastically  equivalent  position.  This  reduces  the 
number  of  independent  moduli  to  nine.  They  are  Cm,  C112,  C1J3,  C123,  CJ33,  C144, 
C166,  C333,  and  C344  .  The  following  relations  hold  for  the  remaining  nonzero  moduli 

C222  =  ^111’  ^355  =  C344 

<^223  ~  ^113»  <^166  ^266  ~  ~  ^112) 

^255  =  ^144r  ^223  —  ^133 
^244  =  C155,  C356  =  ~  ^123) 
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APPENDIX  G 


100  REM  PR.QSRAM  TO  MEASURE  P2L2  FRE8UENCY 

110  REM  CHANGES  AS  A  FUNCTION  OF  APPLIED  LOAD  AND  PLOT  AND 

120  REM  STORE  DATA.  WRITTEN  BY  PROSSER  2/25/8o 

130  REM 


140  REM  VARIABLES 

150  REM 

NAME 

TYPE 

160  REM  SAMPLE  ID. 

St 

CHARACTER(32) 

170  REM  MODE  OF  WAVE 

M$ 

"  (12) 

IBB  REM  DIRECTION  OF  PROPAGATION 

Nt 

"  (12) 

190  REM  •  ”  ’  POLARIZATION 

P$ 

■  (12! 

200  REM  ■  ■  LOADING 

LI 

■  (12) 

220  REM  DATE 

Dl 

CHARACTER (12! 

270  REM  CROSS  SECTIONAL  AREA  OF  SAMPLE  (SO.  IN.IAl 

fl 

300  REM  NUMBER  OF  DATA  POINTS 

B1 

B 

330  REM  FREQUENCY  MEASUREMENTS  (H:) 

F(I) 

ARRAY 

340  REM  VOLTAGE  MEASUREMENTS  FOR.  LOAD  (VQLTSl 

Vdi 

B 

400  REM  STRESS  (Pa) 

Ld! 

ARRAY 

460  REM  DELTA  F/F  AT  EACH  POINT 

Dd) 

1 

520  REM  MINIMUM  DELTA  F/F 

W1 

it 

530  REM  MAXIMUM  DELTA  F/F 

W2 

B 

540  REM  MINIMUM  STRESS  (Pal 

W3 

1! 

550  REM  MAXIMUM  STRESS 

W4 

a 

571  REM  SUM  D(I) 

X3 

1 

572  REM  SUM  (D(I))-‘-2 

X4 

n 

573  REM  SUM  L(I! 

X5 

a 

574  REM  SUM  (L(i)i"2 

X6 

" 

577  REM  SUM  (DdUfKi; 

«7 

n 

590  REM  DUMMY  VARIABLE 

Ct 

CHARACTER (32: 

600  REM 

At 

n 

630  REM  DUMMY  VARIABLE 

A5 

u 

640  REM  INTEGER  COUNTERS 

I,J 

a 

650  REM  DUMMY  ARRAY  FOR  PLOTTING 

Pd) 

ARRAY 

660  REM  "  VARIABLES  FOR  PLOTTING 

AS,A9 

NUMERIC 

670  REM  PLOTTER  » 

B2 

u 

660  REM  REGRESSION  COEF.  FOR  STRESS-DELTA  F/F 

R1 

n 

690  REM  SLOPE  FOR  STRESS-DELTA  F/F 

r\ 

U  i 

H 

692  REM  Y-INTERCEPT  FOR  STRESS-DELTA  F/F 

El 

B 

700  REM 

710  REM  DIMENSION  CHARACTER  VARIABLES 
720  REM 

730  DIM  S$(32),M$(12),NJ(12i,PI{12),LI(121,D$(12) 
740  DIM  C$(32),A$i32! 

750  REM 

7o0  REM  INPUT  INITIAL  VARIABLES 
770  REM 

730  PRINT  'INPUT  SAMPLE  ID.' 

790  INPUT  St 

800  PRINT  'INPUT  MODE  OF  WAVE" 

610  INPUT  Ml 

B20  PRINT  'INPUT  DIRECTION  OF  POLARIZATION’ 
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832  INPUT  P-l 

842  PRINT  “INPUT  DIRECTION  OF  PROPAGATION" 

85S  INPUT  N3 

862  PRINT  “INPUT  DIRECTION  OF  LOADING” 

872  INPUT  LI 

920  PRINT  “INPUT  DATE" 

?18  INPUT  Ds 

990  PRINT  “INPUT  CROSS  SECTIONAL  AREA  OF  SAHPLE  (SO.  IH.}“ 

I3SS  INPUT  Ai 
1301  REN 

1310  PER  DINENSIQN  VARIABLES  HAnHUN  2i8  DATA  PTS.  TAKEN  ABOUT  EVERY 

1320  REN  ONE  SECOND 

1330  DIN  F(20@),V(2SB!,L!28i).D(28ei 

1358  REN 

1368  REN  SET  UP  INSTRQN  AND  TAKE  DATA 
1370  REN 

1380  PRINT  “PRESS  RETURN  TO  START  TAKING  DATA" 

1398  PRINT  "PRESS  BREAK  TO  QUIT  TAKING  DATA" 

1391  PRINT  “THEN  ENTER  RUN  1510" 

1481  INPUT  AI 
1410  FOR  l=!  TO  230 
1428  INPUT  t3:F(I! 

143i  INPUT  @24:VH) 

1458  B1=I 

1478  FOR  J=!  TO  85 
1488  B9=J 
1490  NEIT  J 
1511  NEH  I 

ISIS  REH  CALCULATE  STRESSES  AMD  DELTA  F/F 

1520  REH 

1543  FOR  1=1  TO  B1 

1553  L(!)  =  (V(I)-V(li)/AH(6895'3l»i-l 
1571  D(I)  =  {F(I)-F(I1)/T(1) 

1581  NEXT  I 
1890  RER 

1903  RER  CALCULATE  RAX  AND  R!N  STRESS  AND  DELTA  F/F 

1913  Ni=l 

1923  02=-l 

1933  03=1 

1941  tl4=-I 

1970  FOR  1=1  TO  B1 

19BS  IF  L(Il>tl3  THEN  2010 

1993  !i3=L(I) 

2800  IF  LdXtlA  THEN  2020 
2010  a4=L!I) 

2021  IF  Ddi/U!  THEN  2841 
2838  Hl=Dd) 

2040  IF  DdXUZ  THEM  2iS0 
2058  H2=Ddl 
2180  NEXT  I 
2118  REH 

2128  REH  LINEAR  REGRESSION 
2138  RER 


-165- 


2140  BOSUB  7000 

2150  REH  PLOT  STRESS  VS.  DELTA  F/F 

2160  REM 

2170  SQSOB  2190 

2180  SO  TO  3200 

2190  82=32 

2210  PASE 

2211  PRINT  0B2: 'DATE;  SAMPLE:  3$;”  «AVE  MODE:  "jMI 

2212  PRINT  0B2:”DIR.  OF  PR0PA8ATIQN:  DIR.  OF  STRESS:  " 

2213  PRINT  @32: ’DIR.  OF  POLARIZATION:  ’;P$ 

2220  VIEWPORT  20,120,30,90 

2230  DIM  P(8) 

2240  P(5)=»l 
2250  Pi6)=W2 
22.60  Pf31=5 
2310  Pil)=W3 
2320  P!2)=«4 
2330  P!7)=5 
2340  P5=3 
2350  SQ3UB  5000 
2360  P5=7 
2370  BOSUB  5000 

2380  WINDOW  P':i),P(2),P!5),P(6) 

2390  AUS  SB2:Pi3i,P(7),0,0 

2400  P5=4 

2410 

2420  SOSUB'6000 
2430  P5=8 
2440  AF="" 

2450  BOSUB  6000 
2470  MOVE  i82:L(l!,Dil! 

2480  FOR  .]=!  TO  B1 
2490  DRAW  §B2;L(J),D(J) 

2550  NEXT  J 

2670  AB=«3 

2680  A.9=H4 

2720  MOVE  §B2:A8,«2 

2730  PRINT  §B2:”HHHHHJJ'; 

2740  Cf='DELTA  F/F’ 

2750  FOR  J=1  TO  LEN(C$) 

2760  A4=SES(C$,J,1) 

2770  PRINT  8B2:A4;'HJ’; 

2730  NEXT  J 
2790  MOVE  @B2:A3,WI 
2800  PRINT  ®B2:'JJJ 
2320  C4=’STRESS  IPai” 

2841  PRINT  ?B2:CT 
2850  IF  B2<)32  THEN  2880 
2860  Z9=30 
2670  SO  TO  2890 
2880  Z9=2 
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2098 

KOHE 

2908 

FOR  Z3 

;=1  TO  Z9 

2910 

PRINT 

mi"  " 

292S 

NEH  la 

2930 

PRINT 

i82rHIN.  DEL 

TA  F 

/F: 

2978 

PRINT 

iB2:‘'mN.  STRESS: 

“;li3r 

3030 

PRINT 

§B2: “SLOPE  IS 

a  3  , 

»  ? 

Cir  (Pa"-: 

30^0 

PRINT 

gB2;“Y-INTERCEPT; 

s.fri.  a  /p, 

,  _  A  p 

3140 

PRINT 

iB2rREGRE3Sr 

ON  C 

OEF.5  ';R1 

■j  i  /  V 

INPUT 

3190 

RETURN 

1 

3203 

HONE 

3201 

PAGE 

3210 

PRINT 

'HARD  COPY  OF 

PLOTS^  YES=1  1 

3220 

INPUT 

A5 

3238 

IF  A5( 

)1  THEN  3278 

3240 

PRINT 

"INPUT  PLQTTE 

R  NUMBER" 

3241 

INPUT 

B2 

3258 

60SUB 

2218 

3263 

GO  TO 

3398 

3273 

IF  A5= 

:2  THEN  3300 

3230 

PRINT 

'BAD  CHOICE, 

TRY 

AGAIN*^ 

3293 

so  TO 

3219 

3300 

PRINT 

'STORE  DATA” 

YES 

=1  N9=2' 

TT  1  Q 
j  1  li 

INPUT 

A5 

3329 

IF  A5< 

:>1  THEN  3359 

3333 

SQSUS 

S880 

3343 

GO  TO 

3338 

3353 

IF  A5= 

:2  THEN  3330 

3363 

PRINT 

^9hD  choice. 

TRY 

AGAIN" 

3373 

SO  TO 

3380 

3380 

REN 

3460 

PRINT 

“ANOTHER  RUN'!’ 

YES 

'_i  Mn-7n 

3470 

INPUT 

A5 

3438 

IF  A5< 

::'l  THEN  3500 

3490 

SO  TQ 

1378 

3590 

IF  A5= 

=2  THEM  3538 

3510 

PRINT 

“BAD  CHOICE, 

TRY 

AGAIN' 

3520 

eO  TO 

3460 

3533 

END 

5300 

REN  SUBROUTINE  FGR 

PLOT 

’TING 

5901 

RER  P 

(P5)=HINIHUR  H 

■  OF 

TICS 

5010 

P1=(P' 

(P5_i  )_p  (P5-.2) 

)/P\ 

P5) 

5120 

P2=lg’ 

'INT(LGT(P1) ) 

5030 

P1=P1: 

;  p  1 

5041 

IF  Ri; 

:>2  THEM  5880 

5059 

IF  Pl^ 

=1  THEN  5120 

5960 

P2=2'>l 

j'/ 

5079 

BO  TO 

5128 

5380 

IF  Pi: 

>5  THEN  5118 

mi. 


al  ■’ 


DELIA  F/F;  ' 
m,  STRESS; 


;H2 


:H4;"  (Pa!’ 
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5100  SC  TO  5120 
5110  P2=10tP2 

5120  REM  ADJUST  DATA  MINIMUM 
5130  P1=INT(P(P5-21/P2) 

5140  P3=P2»!Pl+2) 

5150  IF  P3<P(P5-21  THEN  51S0 

5100  P3=P3-P2 

5170  SO  TO  5150 

5180  PiP5-2)=P3 

5190  REM  ADJUST  DATA  MAXIMUM 

5200  P1=INT(P(P5-I)/P2) 

5210  P3=P2t(Pl-2) 

5220  IF  P(P5-1)<P3  THEN  5250 
5230  P3=P3+P2 
5240  SO  TO  5220 
5250  P(P5-11=P3 

5260  REM  P(P51=ADJUSTED  TIC  INTERVAL 

5270  P!P5!=P2 

5280  RETURN 

6000  REM  CABLE  AXIS 

6010  P4=Pi?5-l) 

6020  P(4>=P(1) 

6030  P(9)=P!5! 

6040  P3=ABSiP(P5-31-^P4)  MAX  ABS(P(P5-2)-P4) 

6050  P3=INTiLST!P3}+1.0E-8) 

6060  P2=10''-P3 
6070  PI=PiP5-2)-P4/2 
6080  P(P5)=P!P5)+P4 
6090  IF  PiPSDPl  THEN  6140 
6100  MOVE  gB2:P(4),P(8) 

6110  PRINT  A$; 

6120  PRINT  0B2:  USINS  •-D.2D,S'!P(P5)*P2 

6130  80  TO  6080 

6140  IF  P3=0  THEN  6180 

6150  P(P5)=P1 

6160  MOVE  gB2:P(4),P(8) 

6170  PRINT  082:  USINB  ”2A,+FD,S”:"  £“;P3 
6180  RETURN 
7008  REM 

7810  REM  SUBROUTINE  TO  DO  LINEAR  REGRESSION  ON  DATA 

7020  REM 

7038  X3=8 

7040  X4=0 

7058  X5=0 

7068  X6=0 

7090  N7=0 

7110  FOR  1=1  TO  El 
7128  X3=X3+B(I) 

7138  X4=X4+D(I)"2 
7148  X5=X5+L(I) 

7158  X6=X6+L(I)''2 
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APPENDIX  H 

130  REH  NONLINEAR  DATA  ANALYSIS  PROGRAM  (TRANSVERSE  ISOTROPY! 
110  REM  DETERMINE  7  OF  9  NONLINEAR  COEF.  BY  LEAST  S6UARES  FIT 


120  REM  0SIN6  24  SAC  MEASUREMENTS 
130  REM  WRITTEN  BY  PROSSER  11-10-86 
140  REM 

150  REM  VARIABLES 

160  REM  LINEAR  ELASTIC  STIFFNESSES  Kdl 

170  REM  ”  ”  COMPLIANCES  Sd) 

130  REM  SAC  MEASUREMENTS  Hd,l) 

190  REM  UNCERTAINTIES  IN  SAC’S  Hd,2! 

200  REM  NONLINEAR  COEF'S  Cd,l) 

210  REM  MAL  UNCERTAINTY  IN  NONLINEAR  COEF’S  Cd,2! 

220  REM  PRQB.  "  .  ”  ”  Cd,3) 

230  REM  RECALCULATED  SAC’S  Hd,3) 

240  REM  MAX.  UNCERTAINTY  IN  RECALCULATED  SAC’S  Hd,4l 

250  REM  PRGB.  •  "  •  «  Hd,5! 

260  REM  CQEFICIENTS  OF  EQ’S  (24  EQ.,  7  UNKNOWNS)  A(I,J) 

270  REM  Ad,Ji*Ad,Mi  E(I,J) 

280  REM  INVERSE  OF  Ed,J!  Fd.J! 

290  REM  PRODUCT  OF  E  AND  F  Fld,J) 

300  REM  PRODUCT  OF  F  AND  E  F2(I,J) 

310  REM  CONDITION  NUMBER  OF  E  23 

320  REM  Hd,ll+CONSTANT  TERMS  Yd) 

330  REM  INTEGER  COUNTERS  I,J,l-,M 

340  REN  PERCENT  DIFFERENCE  BETWEEN  MEASURED  AND  CALC  SAC'S  Vd! 
350  REN  "  UNCERTAINTY  IN  CALCULATED  SAC'S  T(I,2) 

1000  REM  INPUT  DATA 

1001  DIM  K(3),S(4),H(24,5),A(24,7!,C!7,3) 

1010  PAGE 


1020  PRINT  "INPUT  LINEAR  ELASTIC  STIFFNESS  MODULI  (6Pa)" 

1030  PRINT  ’11-1,12-2.44-3“ 

1040  FOR  1=1  TO  3 
1050  PRINT  "INPUT  C";I 
1060  INPUT  Kd) 

1070  NEXT  I 
1080  PAGE 

1090  PRINT  "INPUT  LINEAR  ELASTIC  COMPLIANCE  MODULI  (GPA)'’-!’ 
1100  PRINT  “11-1,12-2,13-3,33-4' 

1110  FOR  1=1  TO  4 
1120  PRINT  “INPUT  S";I 
1130  INPUT  Sd) 

1140  NEXT  I 
1150  PAGE 

1160  PRINT  "INPUT  SAC’S  AND  UNCERTAINTIES  (GPa!‘-l" 

1170  PRINT  “122-1,211-2,322-3,311-4,121-5,212-6,312-7,321-8" 

1180  FDR  1=1  TO  8 

1190  PRINT  “INPUT  SAC  “;I 

1200  INPUT  Hd,l) 

1210  PRINT  “INPUT  UNCERTAINTY' 

1220  INPUT  H(I,2) 

1230  NEXT  I 
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1241  RASE 

125i  PRINT  "SAC'S  CQNT'D." 

126S  PRINT  "123-9,213-18,323-11,313-12,131-13.232-14,132-15. 

127S  FOR  1=9  TO  li 

1238  PRINT  "INPUT  SAC  ";I 

I29S  INPUT  H(I,1I 

1383  PRINT  "INPUT  UNCERTAINTY" 

1318  INPUT  H(I,2) 

1328  NEH  ! 

1333  PAGE 

1340  PRINT  "SAC'S  CDN'T." 

1358  PRINT  "811-17,022-18,812-19,121-20,813-21,823-22,831-23 

13o8  FOR  1=17  TO  24 

1378  PRINT  "INPUT  SAC  ”;I 

1388  INPUT  Hi  I, II 

1398  PRINT  "INPUT  UNCERTAINTY" 

1488  INPUT  Hi  1,2! 

1418  NEKT  ! 

1428  REH  CALCULATE  COEF'S  OF  E3.S 
1438  A=g 

1440  Hil,l)=Si2)/(2^!(:!ll! 

1458  A!l,2!=S!ll/(2-JN(l)) 

1468  A(l,3)=S(3!/(2-«:(l!! 

1473  Ai3,l)=A(l,3l 
1488  Ai3,2)=A(l,3! 

1498  A(3,3)=Si4!/(2-JK!l!) 

1503  Ai5,l!  =  iSil)+S(2n'J8.25,''IK!l)-K{2)) 

1518  AI5,2)=-A(5,1) 

1528  Ai5,3)=3.5-JSi31/il;;!l!-K(2)) 

15-j8  A  ^ 5 , 4/ =-A  1 0 ,  j ; 

1531  Ai7,l)=Ai5,3) 

1532  A!7,2!=Ai5!41 

1548  A i7, 3) =@. 5*8 (4i / iK i 1 ! -K i2) ! 

1553  A(7,4!=-Ai7,3) 

1568  Ai9,5)=S(l!/(2!^r.!3)i 
1578  Ai9,6)=SI2i/(2^!K(3i! 

1588  AI9,7)=Si31/i2«K(3!! 

1598  Aill,5)=A(9,71 
168S  A!11,6)=A!9,7) 

1618  A(ll,7!=Si4)/i25K!3)! 

1628  Atl3,5!=A(9,6! 

1633  A(13,6)=A!9,5) 

1648  Ail3,7!=A!9,7/ 

1650  A(15,5)=Ai9,5) 

1668  Ail5,6)=A(9,6) 

1678  A(15,7)=A!9,7) 

1683  B=S(r)+SI2)+S!3) 

1698  D=2iS!3)+Si4) 

1708  A(17,l)=B/i2-flf:(l!; 

1713  Ail7,2)=A(17,li 
1728  Ail7,3l=D/(2il<(ll) 

1733  Ail9,li=@.5^B,MK(ll-K(2)) 

174B  Ail9,2)=-A(19,l) 
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1758  A(19,3)=B/(K(1)-K(2))*0.5 
1768  A(19,41=-A(19,31 
1778  A(21,51=B/(2»K(3)) 

1788  A(2l,6)=A(21,5) 

1790  A(21,7)=D/(2«K(3)) 

1800  A(23,5)=A(21,51 
1810  A(23,6)=A(21,5) 

1820  A(23,7)=A(21,71 
1830  FOR  J=l  TO  12 
1840  FOR  1=1  TO  7 
1850  A!2tJ,I)=A(2»J-l,I) 

1868  NEXT  I 
1878  NEH  J 

1880  REM  ADD  CONSTANT  TERNS  TO  SAC’S 
1890  DIN  Y(24) 

1900  Y(1)=-(H(1,11+S(2)) 

1910  Y(2)=-(H(2,11+S(2)) 

1920  Y(3)=-fH(3, 11+5(3)) 

1930  Y(4)=-(H(4,1)+S(3)1 
1940  Y(5!=-(H(5,1)+S(1)) 

1950  y(61=-(H(6,l)+S(l)) 

I960  V(71=-(H(7,11+S(3)) 

1970  Y!8)=-(H(8,1)+S(3)) 

1988  Y(9)=-(H(9,1)+S(3)) 

1990  Y(10!='(H(10,11+S(3)) 

2080  Y(11)=-(H(11,1)+S(4)I 
2810  Y(12)=-(H(12,1)+S(4)) 

2020  Y(13)=-(H(13,1)+S(l)) 

2038  Y(14)=-(H(14,1)+S(1)1 
2040  Y(15I=-(H(15,1)+S(2)1 
2050  Y(16)=-(H(16,1)+S(2)) 

2060  Y(17)=-(H(17,1)+B+1/(2»K(1))) 
2070  Y(18)=-!H(18,l)+B+l/(2*K(l))i 
2080  Y!19)=-(H(19,1)+B+1/(K(1)-K(2))) 
2090  Y(20)=-(H(20,1)+B+1/(K(1)-K(2))) 
2100  Y(21)=-(H(21,1)+D+1/(2»K(3))) 
2110  Y(22)=-'(H(22,1)+D+1/(2»K(31)1 
2120  Y(23)=-(H(23,l)+B+l/(2tKC3))) 
2130  Y(24)=-(H(24,1)+B+1/(2»K(3))) 
2140  REN  CONFUTE  A(I,N)*A(I,J) 

2150  DIN  E(7,7) 

2160  FOR  N=1  TO  7 

2170  FOR  J=1  TO  7 

2180  E(N,J)=0 

2190  FOR  1=1  TO  24 

2200  E(H,J)=E(H,J)+A(I,N)*A(I,J) 

2210  NEXT  1 
2220  NEXT  J 
2230  NEXT  H 

2240  REN  CONFUTE  INVERSE  OF  E(N,J) 
2250  DIN  F(7,7),F1(7,7),F2(7,7) 

2260  F=INV(E) 

2270  REN  CONFUTE  CONDITION  NUNBER 
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2233 

2293 

2333 

n 

i. 'j  i 

0770 

234S 

2358 

07  ip 

2378 

23S8 

oTog 

2438 

2441 

2453 

2468 

2478 

2488 

24?i 

2538 

2511 

2528 

2533 

os;  ,10 

j.  j'/iii 

2553 

2568 

0^70 
L.  u  /  a 

2500 

2593 

26i8 

2613 

2623 

OiTS 

2643 

2658 

2663 

2678 

2688 

2693 

2788 

2713 

0700 

t.  / 

0770 
L.  I  -Jil 

2743 

2753 

2768 

2773 

2788 

2793 

oDon 

li-UiUai 

28  !i 
2828 
2S3B 


Zl=-lJE+igS 
22=-lJE->ii8 
FOR  !=1  TO  7 
FOR  TO  ; 

IF  Ed.JXZl  THEfi  2349 
Z1=E(I,J1 

IF  F(I,J1<Z2  THEf!  2389 

Z2=F(LJ! 

fJElIT  .] 

ME>T  I 

Z3=ft8S(Z25Zn 
F1=E  HPY  F 
F2=F  HPY  E 

REM  CQriPUTE  HOHLIMEhR  CQEF'S 

FOR  1=1  TO  7 

C!I,l!=S 

C!l,2)=9 

C':i,3)=Q 

FOR  J=1  TO  24 

FOR  L=1  TO  7 

CiI,!)=C(I,ll^F'L,noft(J,L)sY(J) 
C(I,21=C(!,2ivABS':F(L,IiiA!J,Ll)'JH(J,2i 
NEXT  L 
NEXT  J 

C(I,3)=C(I,2!/23"2.5 
NEXT  ! 

REM  RECALCULATE  SAC  VALUES  USIH6  DERIVED  VALUES 

HdjSX-'E-ndlXSiEl+SiUiC.iZ.D+SiCXCd.D+SiTl-JCiS,!)! 

Hd,3)=H(l,3)/(2-5r;d)) 

H  !3 ,3)  =-  im.  d !  -IS  (3)  +S(3)  ■>  (C  i  1 , 1  i  +C  (2 ,  D  i  *3  (45  ->2  !3, 1  d 
Hi3,3)=H(3,3!/(2-iKd)5 

H{5,3)  =  (?:d)-K(21)3S{I)+a.2SdSd)+3i2)i-}(Cd,l!-C(2,!!! 
Hi5,3)=H(5,31+9,5-::-S(3)->iC(3,l)-C(4,l)) 

H(5,3)=-H(5,3)/(;<d)-l<(25) 

H(7,3)  =  (Kd]-K(2))5S(3)+8.5-iS(3)-dCd,l!-C(2,l)! 

H{7 ,3)  =H (7,3HS=  5-iS (4) i (C (3 , 1 5 -C (4, 1 5  5  ‘ 

Hi7,3)=-H(7,3)/(Kd)-K(2!5 

H(?,35=-(25(((3)iS!3)-fSa!5C(5ii!+S{2)4C(e,li+S(3)JC(7,l)!/(2iK(31) 

H  ( 1 1 ,35  =-  (2Hd3)-3S(4i  +S  !3}  -J  iC  (5, 1 5  +C  io ,  1 5  5  (3)  iC  (7 , 1)  5  /  (2-JK  (3)5 

Hd3,3)=-i2d(!35vS(l!-:-Sd:'5C(6,l)iS{2)'JC':5,l!+3(3X-C{7,l!)/(2JKi3)) 
H  (15,35  =-  !2-X(3:'  iS  (25  +S  (15  ->2  (5, 1 1 +S  (2)  iC  (a ,  i )  +3  (3)  -sC  (7 , 1 5 )  /  i2->K.  (3) ) 
Hll7,3!=-(H2®f:(I)5B+35(Cd,ii+C(2,l))+D^C(3,l)i/(2^K(l)i 
H(19,35  =  l-MX;(Ii-:<(25j5B 

Hd9,3)=Hd?,35+8=55(B-:H:Cd,i:'-C(2,i!)+D-i(C(3,l)-C(4,i5)! 

Hd9,35=-Hd9,35/(Kd)-X;;(25) 

H  (21 .3)  =-  { 1+2'X  (3XD-!-B-i  (C  (5, 1 )  +Ci6 , 1 5 )  +D-iC  (7 , 1 5 )  /  (2->5(  (35 ) 

H(23,35  =-!U25K(3)-iB+B5(C(5,l!+C(4,l)5-i-D-3C(7,!15/(2-iKi3)) 

FOR  !=1  TO  12 
H(2il,3)=H(2iI-!,3) 

NEXT  I 
DIM  9(7,25 
FOR  !=1  TO  2 
FOR  J=1  TO  7 
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2646  W(J,I)=ABS(C{J,1+1!/C(J,1)) 

2858  NEIT  J 
2868  NEXT  I 
2878  FOR  I=l  TO  2 
2880  FOR  J=l  TO  4 

2898  H(J,I+3)  =  CH(i,I!+W(2,I)+H(3,I)HfH(J,3) 

2988  NEXT  J 
2918  FOR  J=5  TO  8 

2928  H(J,I+3)  =  (W(l,I)+i<l(2,I)+H(3,I)+W(4,Il)*H(J,3! 

2938  NEXT  J 

2948  FOR  J=9  TO  16 

2950  K(J,I+3)  =  (H(5,r)+N(6,l)+Hi7,I))tH(J,3) 

2960  NEXT  J 
2970  FOR  J=17  TO  18 

2980  H(J,I+3)=!H(l,I)+W(2,I)+H!3,I))tH(J,3) 

2998  NEXT  J 
3008  FOR  1=19  TO  28 

3010  H-:j,I+3!  =  ’:H(l,n+lii(2,!)+W':3,I)+«(4,IlHtH(J,3) 

3028  NEXT  J 

3030  FOR  J=21  TO  24 

3048  H(J,I+3)  =  (W!5,Ii+i!i(6,Il+«(7,I))fH(J,3) 

3058  NEXT  J 
3060  NEXT  I 
3878  FOR  1=4  TO  5 
3080  FOR  3=1  TO  24 
3090  H(J,n=ABS(H(J,I)) 

3180  NEXT  J 
3110  NEXT  I 
3120  DIN  V{24) 

3130  FOR  I=l  TO  24 

3148  V(ll=ABS((H(I,l)-H(I,3))/H(I,l)) 

3150  V(!.i=V(I)*10B 
3168  NEXT  I 
3170  DIN  T(24,2) 

3180  FOR  I=l  TO  2 
3190  FOR  J=l  TO  24 

3200  T(J,I)=H(J,l+3)/H(J,3)H00 

3201  T(J,I)=ABS(T(J,I)) 

3210  NEXT  J 

3220  NEXT  I 

3230  REN  PRINT  OUT  DATA 

3240  REN 

3250  N=32 

3260  PA6E 

3270  GOSUB  3380 

3288  PRINT  ’PRINTED  COPY?  YES  =1  NO  =2  “ 

3290  INPUT  25 
3300  IF  25=1  THEN  3340 
3310  IF  25=2  THEN  4010 
3320  PRINT  "BAD  CHOICE” 

3330  GO  TO  3288 

3340  PRINT  ‘INPUT  PRINTER  S" 

3358  INPUT  N 


3363  GOSUB  3338 
3378  SO  TO  4g!3 

33S0  PRIMT  813;  "NOfiLIHEAR  DATA  ANALYSIS  FOR  TRANSVERSE  ISQTRG! 
3393  PRINT  "  = 

3409  PRINT  INPUT  DATA" 


34  li  PRIiilT 

’dS'ia 

342e  PRINT 

§N?"LIN 

EAR  ELASTIC  HQDOLr 

3430  PRIilT 

ON? "  ' 

3440  PRINT 

§N?^C11 

=  "jKfi);"  (SPa) 

Si 

1  =  ";S(l)r’  (GPa 

3450  PRINT 

§NrX12 

=  ";K{2i;"  (SPa) 

SI 

2  =  ”;S-:2);'  (SPa 

3460  PRINT 

§N?'’C44 

=  ";K(3);”  (GPa) 

SI 

3  =  ";3(3);“  (SPa 

347i  PRINT 
3483  PRINT 

•J  ” 

”;S(4);"  (SFaS  -l 

3493  PRINT 

§NrVEL 

QCiTY  DERIVATIVE  DA 

TA  ( 

6Pa)'-l" 

3500  PRINT 

3510  PRINT 

IN:"  L 

Hi22  = 

3520  PRINT 

H2il  = 

";H(2,2) 

3531  PRINT 

H322  =  =:,H(3,l)r 

‘‘  /  “ 

“■un  q; 

3540  PRINT 

ON? ^4. 

H3!I  =  '■;H(4,l)j‘= 

■r  /  - 

'';H(4,2) 

3550  PRIN- 

HI2i  =  ";Hi5,llr 

>/- 

"^H(5,2! 

3560  PRINT 

mi^h. 

H212  = 

";H(6,2) 

3570  PRINT 

H312  =  ";H(7il!:= 

>/- 

";H(7,2) 

3580  PRINT 

H321  =  ";H(8,!)r 

’(*  /  - 

";H!a,2) 

359S  PRINT 

0N?  '"P, 

H123  =  ";H(95l!j"  ■ 

3 ,  u  Q  n\ 

3600  PRINT 

§N?^IO. 

K2I3  = 

"F  /“ 

"5H(18,2) 

3610  PRINT 

§N?ML 

H323  = 

";H(11,2) 

3623  PRINT 

H313  =  "iH>:i2,l);’ 

+  /  - 

3633  PRINT 

iNs^i3, 

H131  =  “;H{13,l)r 

-!•/- 

3 , p / 1 7  n\ 

,  n  i  j  ^  i  / 

3648  PRINT 

2  *  .-I 

■dt-lj  I'f, 

H232  = 

/  - 

";H(1452) 

3653  PRINT 

i3r-‘=  2 1 

^  1  i  3  1  W  S 

Hi32  =  ';K(15,15;" 

.j./- 

";H(15,2) 

3668  PRINT 

mint. 

H231  =  "jHdo,!);" 

i/- 

";Ha6,2) 

3670  PRINT 

H011  = 

'>/- 

";H(17,2) 

3608  PRINT 

H022  =  ”;HT18,1);= 

3698  PRINT 

H312  = 

/  - 

3710  PRINT 

§r.!?=283 

H@2I  =  =;H(29il):= 

->/- 

";H(2i,2) 

3710  PRINT 

ON? "21, 

HSi3  =  ";H(21,l)r 

h'.  /  _ 

‘  / 

";H(21,2) 

3723  PRINT 

§N?"22. 

HS23  =  “;H(22,i:';  = 

/  - 

";H(22,2) 

3730  PRINT 

iN?"23, 

H831  =  “;Hi23,l);= 

t  /  - 

^;H(23,2) 

3740  PRINT  §N!“24.  HS32  =  =;H!24,!);"  “;H{24,2) 

375S  PRINT  mi"  " 

3768  PRINT  ^Ih-THIRS  ORDER  ELASTIC  CONSTANTS  DERIVED  FROf1  LEA 
3778  PRINT  8N;"SeUARES  FIT  AND  HAX.  ERROR  AND  PROD.  ERROR  (SP 
3733  PRINT  gfh'  = 


i790 

PRIflT 

ON; 

"Clll 

=  hC(L 

"lC!l 

t/-  ^ 

:8SS 

PRIflT 

“C112 

-  ”;C(2: 

";C(2 

j  Z) 

,n  . 

>819 

PRINT 

0N; 

“C113 

=  “:,C(3, 

’;C(3 

,2) 

,3  . 

:82e 

PRIfJT 

§N; 

"C123 

=  “;C(4, 

,i)r  ^z- 

=  ;C(4 

,2) 

,3  .L  /  _  n 

J 

830 

PR  HIT 

§N; 

"C144 

-  a  =  r 

,1);"  '^Z'  ^ 

HC(5 

■3 

,3  4-/- 

} 

:34e 

PR  HIT 

■§N; 

=0155 

j)r  •^Z" 

'•;C(6 

,2) 

,3  /_  El 

E  ■  f 

850 

PR  HIT 

§N; 

“C344 

=  ";C(7, 

1;;"  ^Z«  ^ 

■lC(7 

,2) 

;»  f/-  " 

86S 

PR  HIT 

IN; 

“C133, 

,C333  UfJDETERHHIED 

DUE  ' 

TO 

LACK  OF 

;C(6,3) 

:C(7.'3) 


3373  PRINT  " 

!390  PRINT  8N: ’’RECALCUIATFn  UFi  ncriY  nPRTUATiypq  iiqrMR  FTTTFFl 
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3398  PRINT  §N:'AND  HAL  AND  PROS.  ERRORS  IN  SAHE  ORDER  AS  BEFORE' 
3900  PRINT  §N:"  * 

3910  FOR  1=1  TO  24 

3928  PRINT  @N:I;’.  ”;H(I,3);'  +/-  ”;H(I,41;"  +/-  ';H(I,5) 

3930  NEXT  I 
3940  PRINT  8N:=  ’ 

3950  PRINT  @N: 'PERCENT  DIFFERENCES,  HAX.  AND  PROB.  UNCERTAINTIES” 
3960  PRINT  8N:"  ' 

3970  FOR  1=1  TO  24 

3980  PRINT  §N:Ij'.  ";V(I1;T(I,1);T(I,2) 

3990  NEXT  I 

3991  PRINT  §N:'  " 

3992  PRINT  §N: "CONDITION  #  FOR  INVERTINB  HATRIX  IS:  'iZO 
4000  RETURN 

4010  END 


lia  RErt  DETERHIME  i5  QF  28  flOfiLIMEAR  COEF.  BY  LEAST 
120  RE:t  OSIfJB  24  SAC  MEASUREHE:4TS 
13B  REH  HRITTEiJ  BY  PROSSER  11-13-86 
HO  REfi 

150  REii  VARIABLES 


iiUL,t:L/ 

SQUARES  FIT 


16S 

REF!  LINEAR  ELASTIC  STIFF 

NESSES 

K(I) 

11% 

REF!  =  '  CQHPl 

lAflCES 

Sdl 

ISS 

REF!  SAC  flEASUREFlEFITS 

Hdji} 

193 

REFi  UNCERTAINTIES  IN  SAC 

'S 

H{I,2) 

i-'-U-J 

REF!  NONLINEAR  CQEF’S 

C(I,i) 

213 

REH  HAL  UNCERTAINTY  IN 

NGflL^lEAR  CQEF'3 

CdjZ) 

22S 

REH  PR08.  “ 

C  3 

C(d3) 

77.0 

i-jij 

REH  RECALCULATED  SAC'S 

Hd,3) 

24S 

REH  HAL  UNCERTAINTY  IN 

RECALCULATED  SAC'S 

H(I*4) 

REH  PF!QB  = 

n  j 

H(d5} 

263 

REF!  CQEFICIENTS  QF  EO'S 

■■■ii-ir  t u 3 , 1 D  y y li 0 y y s ) 

Ad  j  J) 

27S 

REF!  A(I,J!JA(I,ri! 

£d,J) 

233 

REF!  INVERSE  OF  Eli,. I) 

Fd,J) 

“aQO 

REH  PRODUCT  QF  E  AND  F 

Fldd) 

303 

REri  PRODUCT  QF  F  AND  £ 

F2d,J) 

T|q 

REF!  CONDITION  MUHBER  OF 

l- 

L  J 

TOO 
i.  'il 

REFI  H(!,!L-CQKSTANT  TERN 

c 

vi 

Yd) 

330 

REFI  INTEGER  COUNTERS 

dJ,L,li 

348 

REF!  PERCENT  OiFFEREFiCE  B 

ETtlEEN  HEASURED  CP 

ic  SAC ^3  y 

■j  jI; 

REH  =  UNCERTAINTY 

IN  CALCULATED  SAC'S 

T{d21 

1000  RER  iliPOT  BATA 


•1  ont 

i  'iitj  i 

DIH  KiSijSia 

),H!24,5!,A(24,15),C!!5,3 

1013 

PAGE 

1020 

PRINT  dNPUT 

LINEAR  ELASTIC  STIFFNESS 

1030 

PRINT  dl-i,: 

22-2,33-3,44-4,55-5,66-6= 

1040 

FOR  1=1  TO  6 

1050 

PRINT  dNPUT 

C”;I 

1060 

INPUT  Kd) 

1078 

NEH  I 

less 

PAGE 

1098 

PRINT  dHPUT 

LINEAR  ELASTIC  COHPLIANCI 

1103 

PRINT  di-i,: 

22-2,33-3,23-4,13-5,12-6= 

i  M  0 
i  i  'ij 

FOR  1=1  TO  6 

1  1  IQ 

PRINT  dNPUT 

SL! 

i  i  Tp 

1  -J  •  j 

INPUT  3d) 

1140 

NEd  I 

1  i  70 

1  i  -f  ii 

PAGE 

ll6i 

PRINT  dNPUT 

SAC'S  AND  UNCERTAINTIES 

1171 

PRINT  d22-l 

5  i  1  i  "*  2  5  •->  i  2 "  -J  j  -j  1  i  **  V ,  1 2 1  ■■  J  J . 

1183 

FOR  1=1  TO  8 

1193 

PRINT  dNPUT 

SAC  d  I 

1  non 

i  i-LJU 

INPUT  Hdd) 

1210 

PRINT  dNPUT 

UNCERTAINTY^ 

1220 

INPUT  Hd ,2) 
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1240  PAGE 

1250  PRINT  “SAC'S  CONT'O." 

1260  PRINT  •123-9,213-10,323-11,313-12,131-13,232-14,132-15,231-16“ 
1270  FOR  1=9  TO  16 
1280  PRINT  "INPUT  SAC 
1290  INPUT  H(I,1) 

1300  PRINT  "INPUT  UNCERTAINTY" 

1310  INPUT  H(I,2) 

1320  NEXT  I 
1330  PAGE 

1340  PRINT  "SAC’S  CON'T." 

1350  PRINT  •011-17,022-18,012-19,021-20,013-21,023-22,031-23,032-24" 
1360  FOR  1=17  TO  24 
1370  PRINT  "INPUT  SAC  ■;! 

1380  INPUT  H(I,1) 

1390  PRINT  "INPUT  UNCERTAINTY” 

1400  INPUT  H(I,2) 

1410  NEXT  I 

1420  REM  CALCULATE  CQEF'S  OF  E0.S 
1430  A=0 

1440  A(1,4)=S(11/(2*K(2)) 

1441  A(l,5)=S(6)/(2*!(;(2)) 

1442  AIl,61=S(5!/(2iK(2)) 

1443  A(2,1)=S(6)/(2»K(1I) 

1444  A(2,2)=S(21/(2«(11) 

1445  A!2,3)=S(4)/(2tK(l)) 

1446  A(3,4)=S(5)/(2iK(2!l 

1447  A(3,5)=S(4)/(2iK(21) 

1448  A(3,6!=S(31/(2iK(21) 

1449  A(4,1)=S(5)/(2»K(1)) 

1450  A(4,2)=S(4)/(2tK(ll) 

1451  A(4,3)=S(3)/(2iK(l)) 

1452  A(5,13)=S(1)/(2*K(6)) 

1453  A(5,14)=S(6)/(2»K(6)I 

1454  A(5,15)=S(5)/(2*K(6)) 

1455  A(6,13)=S(6)/(2»K(6)1 

1456  A(6,14)=S(2)/(2»K(6)1 

1457  A(6,15)=S(4)/(2»K(6)1 

1458  A(7,131=S(5)/(2*K(6)) 

1459  A(7,14)=S(4)/(2tK!6)l 

1460  A(7,15)=S(3)/(2*K(6)) 

1461  A(8,13)=A(7,13) 

1462  A(8,14)=A(7,14) 

1463  A(8,15)=A(7,151 

1464  A(9,71=S(l)/(2tK(4)) 

1465  A(9,8)=S(61/(2*K(4)) 

1466  A(9,9)=S(5),'(2tK(4)) 

1467  A(10,10)=S(61/(2»K(51) 

1468  A(10,ll)=S(2)/(2tK(5)) 

1469  A(10,12)=S(4i/(2»K(51) 

1470  A!l!,7)=S(5)/(2tK(4)) 

1471  A(11,8)=S(4)/(2*K(4)) 

1472  A(ll,9)=S(3i/(2»K(4>! 
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1473  A(12,li)=S(5)/(25K(5)) 

1474  A(!2,11)=S(4)/(2«K(§)1 

1475  fl(12,12)=S(3)/(23K(5)l 

1476  Ail3,ig)=S(l)/(23i(f5)) 

1477  A! 13, 11) =3(6) /(25K (5)) 

1478  A(13,12)=S(5)/(23l((3)) 

1479  A(14,7)=S(6)/(25i((4>) 

1489  A(14,8)=S(2)/(25i((4)) 

1481  A(14,9)=S(4]/(2^K(4)1 

1490  A{15,71=A(9,7) 

1509  A(15,ai=A(9,8) 

1518  A(15,9)=A(9,9) 

1528  A(16,li)=A(19iiai 
1538  A(16,ll)=A(li,ll) 

1540  A(16,12)=A(18,12) 

155S  B1=S(1H-S(6)^-S(51 
1561  B2=S(61+S(2)+S(4) 

1570  B3=S(51+S(4)-i-S(3) 

1580  A(17,1)=BI/(231<(1)) 

1590  A(17,2)=82/!23l((l)) 

1680  A(17,3)=B3/(23K{1)) 

161S  A(18,4)=B1/(23K(2)) 

1620  A(18,5)=B2/(23K!2)) 

163i  A(18,6)=B3/(23K(2)) 

1640  A(19,13)=B1/(25K(6)) 

1658  A(19,1,4)=B2/(25!<(6)) 

1661  A!19,15)=B3/(25K(6)) 

1670  A(2S,13)=A(19,13) 

1680  A(2i,14)=A(19,l4) 

1690  A(20,15)=A(19,15) 

1780  A(21,l@)=01/(23!((5)) 

1718  A(21,ll)»B2/(2iK(5)) 

1728  A(2I,12)=B3/(25i((5!) 

1730  A(22,7)=81/(25K(4)) 

1740  A(22,a)=B2/(25K(4)l 
1750  A(22,9)=B3/(25K(4)) 

1760  A(23,li)=A(21,lQ) 

1770  A(23,11)=A(21,11) 

178Q  A(23,12)=A(21,12) 

1799  A(24,7)=A{22,7) 

1809  A(24,8)=A(22,8) 

1810  A(24,91=A(22,9) 

1880  REH  ADD  CONSTANT  TERRS  TO  SAC'S 
1890  DIR  7(24) 

1980  Y(1)=-(H(1,1)+S(6)) 

1910  Y(2)=-(H(2,11+S(.6)) 

1920  ■f(3)=-(H(3,l)3S(4)) 

1930  Y(4)=-{H(4, 1)38(5)) 

1940  Y(5)=-(H(5, 1)38(1)) 

1950  Y(6)=-(H(6, 1)38(2)) 

1968  Y(7)=-(H(7, 1)38(4)1 
1978  Y(8)=-(H(8, 1)38(5)) 
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1980  Y(9)=-(H(9,1)+S(5)) 

1990  Y(10)=-(H(10,1)+S(4)) 

2000  Y(11)=-(H(11,1)+S(3)) 

2010  Y(12)=-(H(12,1)+S{3)) 

2020  Y(13)=-(H(13,1)+S(l)l 
2030  Y(14)=-(H(14,1)+S(2)) 

2040  Y(15)=-(H(15, 11+8(6)1 
2050  Y(16)=-(H(16, 11+8(61) 

2060  Y(17)=-(H(17,1)+B1+1/(2»K(1))) 

2070  Y(18)=-(H(18,l)+B2+l/(2tK(2)l) 

2080  Y(19)=-(H(19,l)+B2+l/(2+K(6))) 

2090  Y(20)=-(H(20,1)+B1+1/(2»K(6))) 

2100  Y(21)=-(H(21,1)+B3+1/(2»K(5)!) 

2110  Y(22)=-(H(22,1)+B3+1/(2»K(4))) 

2120  Y(23)=-(H(23,1)+B1+1/(2*K(5)1) 

2130  Y(24)=-(H(24,11+B2+1/(2»K(4))) 

2140  REH  COMPUTE  A(I,H)*A(I,J) 

2150  DIM  E(15,15) 

2160  FOR  H=1  TO  15 

2170  FOR  J=1  TO  15 

2180  E(M,J1=0 

2190  FOR  1=1  TO  24 

2200  E(M,J)=E(M,J)+A(I,M)»A(1,J) 

2210  NEXT  I 
2220  NEXT  J 
2230  NEXT  M 

2240  REN  COMPUTE  INVERSE  OF  E(M,J) 

2250  DIM  F(15,15) 

2260  F=INV(E) 

2270  REM  COMPUTE  CONDITION  NUMBER 

2280  Z1=-1.0E+100 

2290  Z2=-1.0E+100 

2300  FOR  1=1  TO  15 

2310  FOR  J=1  TO  15 

2320  IF  E(I,J)<Z1  THEN  2340 

2330  Z1=E(1,J) 

2340  IF  F(I,J)<Z2  THEN  2360 
2350  Z2=F(I,J) 

2360  NEXT  J 
2370  NEXT  1 
2380  Z3=ABS(Z2+Z1) 

2440  REM  COMPUTE  NONLINEAR  CQEF'S 

2450  FOR  1=1  TO  15 

2460  C(I,11=0 

2470  C(I,2)=0 

2480  C(I,3)=8 

2490  FOR  J=1  TO  24 

2500  FOR  L=1  TO  15 

2510  C(I,l)=C(I,ll+F(L,I)tA(J,L)»Y(J) 

2520  C(I,2)=C(I,21+ABS(F(L,I)»A(J,L))«H(J,2) 
2530  NEXT  L 
2540  NEXT  J 

2550  C(I,3)=C(I, 21/23^0.5 
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^•Li'yU 

2373 


ECsLClILIiTE  SAC  vAlJES  USiilC  VALUES 


2353  li  i  i  5  3 )  ( 2'3 ( 2 )  os  ( 6 )  -^S  i  1 1  -JC  i  4 ,  i }  -i-S  ( 6  5  5C  i 5  i  i )  ■> a  ( 5 )  -iC  ( a ,  1 )  1 

25G1  Hil'3}=H(!531/(2o:((2)) 

2582  H (2;3) -- (2o:( ( 1  i oS ie) oS (65 oC i i  J ) og !2) oC (2,  i 5 oS (4)  oC !3,  i )  1 
2333  li(2,35=!'l(2j3)/(2oK(!5) 

2384  HI3,3)=“(2o:((2)oS!45vS(55oC(4il)-:-S(4)oC(3,i)-!-S(3)oCi6,151 
H(3, 3)41(3, 35/(202(2)) 

H(453i=-(2o!;(l)oS(3)oS(5)oCiiil5+S(4)oC[25l)->S(3)oC(3,l)! 
H(4'53)=5(4,35/(2o:(il5! 

!-H5 ,3)  ="  (2'0;(  (6)  OS  ( 1  )o.3  ( i )  oC  i  13 , 1 )  vS  (6)  oC  (14,1)  oS  (5)  oC  1 15, 1 ) ) 
;i(5i3)4!(5,31/(2o:((6;) 

5(6,3)=-(2o:((6)oS(2)oS(65oC(13,l}oS(2)oC(14ii)oS(4)oC(15,l)) 
(1(6,3)41(6, 3)/(2oK(6)) 

H(7,31="(2o:;(6)oS(4)oS(5)oC(13,l)-:-S(4)oC(14,l)-:-B(3)oC(i5,l)) 
H(7i3541(7,35  /(2o;((6)) 

1-1  (8 ,3)  =-  (20(1  (6)  os  (55  oS  (5)  oC  ( 13,  i )  oB  (4)  oC  (14,1)  -i-S  (3)00(15,1)) 
8(8,3)41(8, 3)/(2o:((6)) 

H(9,3!=-(2o:((4!oB(3)oS(l)oC(7,i)o8(6)oC(3,i)oS(3)oC(9,l)) 
H(fi3)44('?,31/(2o:((4)) 

H(i-353!^“(2o:<(5)5S(4)oS(6)oC(l’j,i)vS(2)oC(il,l)oS(4)oC(12,l)) 
HdO, 3)41(13, 3)/(2oX(51) 

H(iL3)=-(2ol((4)oS(3)oS(5)oC  (7, 1)05(4)00(8,1)  08(3)00(9, 11! 
!f'ii,3)^H(il,3)/(2o:ii4)) 

H(i2,3)^-(2oi((3)oS(3)oS(5)oC(19,l)oS(4)oC(ll,i)oS(3)oC(12,I!; 
11(12,3)41(12-3) /(2oK(51) 

1-1(  13,3)=- (2o;a3)oS(l)oS(l)oC(lS.l)vS(6)iC(ii,:i-'id)oC  (12,1)) 
(■1(13, 3)41(13, 3)/(2o;((5i) 

H  ( 1 4 , 3 1  =-  ( 25:1  (4)08(2)05(6  5  o  C  i  7 , 1 )  +  S  ( 2 )  o  C  ( 8 , 1 )  o  S '  4 )  o  C  ( 9 , 1 ) ) 
l-i(l4,3!41(14i3i'(2or(4)) 

11 1 15, 3)=-  (2iK  (4)  oS  (-6)  oS  ( l)oC  (7, 1 !  oS  (6)  oC  (8, 1 )  +3  (5)  oC  (9 , 1 ) ) 
H(I5,3!4l(15,3)/(2oK(4)) 

H(16,3)=-{2oi((5)oS(6)+S(d)oC(lS,l)oS(2)oC(lI,lloS(4)oC(12,i)i 
H(i6, 3)41(16, 3)/(2o:((5)) 


2580 

2586 

2337 

2333 

2309 

2392 

i  i 

2592 

2593 

2594 
2395 

2596 

2597 
2593 
2599 
2sC3 
26Q1 
2632 

26  j-j' 

2604 

2635 

2636 

2637 
26S3 

"■i  QQ 


2620  H(i7,3)=-(lo2oBioK(i)-)BloC(l,I)oB2oC(2,l)oB3oC(3,U) 
2639  H(i7,3)=8(!7j3)/(2o:((l)) 

2649  li ( iB,’3)  =- ( lo25l( (2) o02o0loC (4 , ! l+e2oC(3, 1 ) 0-3300(6,1 ! ) 
2633  11(18, 3)41(18, 3)/(2-oK(2)) 

2669  H ( 19, 3)  =-  ( ioZ-oK (6)  -oBB-i-BloC  (13,1) oB25C  (14,11  oB3oC  (15,1)! 
2679  H(19, 3)41(19, 3)/(2o’((6)) 

2685  H(20,31=-(io2-5:((61-2Bio3ioC(13,l)oD2-oC(14,llog3-JC(iS,l)) 
26-93  H(2S, 3)41(28, 31/(2o;((6)! 

2733  H  { 2 1 , 3 1 =-  ( ! o2o;(  ( 5)  0B30B 1  -OC  ( 1 3 , 1)  0320C  ( 1 1 , 1 )  oB 3-jC  ( 1 2 i  1 ) ) 
2713  H(2i, 3)41(21, 3)/(2-o:((5)) 

2720  H  (22 , 3 1  =-  ( 1  o25X  (4 )  oB3oill -oC  (7,1!  oB2-JC  (8,1)08300(9,1)) 
273S  H(22, 3)41(22, 31/(2-o(((4!) 

2749  H(23,31=-(lo2o!((5)-oBl->B!-oC(13i  1)08200(11, 1)08350(12,1!) 
2753  H(23,3)41(23,3)/(25:((31) 

2763  11  (24,3) =-  ( io2o!(  (4)  oB2-:'3loC  (7 , 1 !  o32oC  (S,  i )  oB3oC  (9,1)! 
2773  H(24, 3)41(24, 3)/(2o;((4)) 

2810  Dlil  11(15,2! 

2823  FOB  i=l  TO  2 
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2838  FOR  J=1  TO  15 

2840  !l!J,I)=HBS':CiJ,IHi/C':.J,lii 

2350  NEXT  J  • 

2360  NEXT  I 
2370  FOR  1=1  TO  2 

2880  H{l,r+3)=i«(4,I)rW(5,I!+H(6,!l!»H(l,3) 
238!  H!2,I+3)  =  (lil(l,I}t«!2,!)+W!3,I)i*H(2,3) 
2332  Hi3,I+3)  =  («!4,I!+W(5,n+W(6,I!UH!3,3! 


2833 

H(4, 

1+3! 

!+«(2,II+H!3, 

D) 

♦H{4,3) 

'>001 

i-UUT 

H(5, 

1+3! 

=(H(13, 

1!+«(14, 

D+Wi 

15, 

I)  )+H(5,3.- 

'JOpex 

H(6, 

1+3) 

=(«{13, 

I)+W(14, 

D+W( 

15, 

Ili+H!6,3 

2386 

H(7, 

1+3! 

=  (Sil(13, 

I)+«(I4, 

I)+i»( 

15, 

I))+H(7,3: 

2387 

H(8, 

1+3! 

=(«(13, 

I!+W(14, 

I!+H( 

15, 

I!!+H(8,3 

2S83 

H(?, 

1+3) 

)+iij(8,n 

+W{9, 

I))#H(9,3) 

2889 

H(10 

,1+3 

)  =  (if(10 

,i)+i^(ii 

,I)+i^ 

/ 1  ^ 
\  i  L 

,I!1+H(10, 

2S?9 

Hdl 

,1+3 

!=(N(7, 

i!+i4(8,I 

)+W(9 

,n 

)»H(!1,3! 

2891 

H(12 

,1+3 

,I}+«(1I 

,I)+|i| 

(12 

,I))+H(12. 

2392 

H(13: 

,1+3 

J}+i^(ll 

,I1+W(12 

,I)1*H(13, 

2393 

H!14; 

,1+3 

!!+W(8,! 

)+y(9 

!+H(14,31 

2894 

H(15, 

,1^3 

)=(W(7, 

I)+y(8,I 

}+y(9 

T  \ 

J 

)*H{15,3) 

2395 

H(16, 

,1+3 

)=(W(18 

,I)+i»l(Il 

,i)+w 

(12 

,I!!*H!16, 

2896 

H(17, 

,  I+3i 

[)+i«((2,r)+«t3 

,1) 

!*H(!7,3! 

2397 

H(13, 

.1+3: 

I!+ii(5,l 

)+«-:6 

,1! 

)+H(18,3) 

2993 

H!!9, 

L+3] 

f=(;^(i3< 

.I)+W(14 

,I)+W 

(15 

,I!!tH(19, 

2399 

H!29, 

,1+3! 

i=^N(13 

,  I)+?i(14 

,  I)+9l 

(15 

,I!!+H(20, 

2900 

H(21, 

1+3) 

'=(^(13, 

,I!+i4!il 

,I!+W 

(12 

,I!1»H(21, 

2901 

H!22. 

,1+31 

!!+il(3,! 

)+^(9, 

T  > 

,  i  1 

!tH(22,3! 

2902 

H(23, 

1+3) 

,I)+i^(Il 

,!)+«' 

!  \  n 

1  1  . 

,I!!+H!23, 

2903 

H'24, 

1+3) 

[!+W(8,l 

)+y(9, 

,n' 

)+H(24,31 

3360 

NEXT 

T 

X 

3070  FOR  1=4  TO  5 
3380  FOR  J=1  TO  24 
3390  H!J,!)=ABS!H!.],i)! 

3130  NEXT  J 
3113  NEXT  I 
3120  DIN  '7(24! 

3133  FOR  I=!  TO  24 

3140  v(!)=ft3S((HiI,l}-H(I,3!)/H(I,ll) 

3153  v;i)=V!I!»103 

3160  NEXT  I 

3173  DIN  T!24,2i 

3180  FOR  1=1  TO  2 

3190  FOR  1=1  TO  24 

3200  T(J,I!=H(J,!+3)/H!l,3)tl08 

3201  T!.:,r!=ABS(T(J,I); 

3210  NEXT  J 

3228  NEXT  ! 

3230  REN  PRINT  OUT  BATA 

3240  REN 

3250  N=32 

3260  PAGE 

3270  30SU8  3380 

3283  PRINT  "PRINTED  COPY?  YES  =l  NO  =2  ' 
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32S8 

INPUT 

7r 

lJ 

3300 

IF  Za= 

=1  THEN 

3343 

3313 

IF-  Z5^ 

=2  THEN 

4813 

3320 

PRINT 

’BAD  Ci 

HO  ICE’ 

3333 

GQ  TO 

3280 

3343 

PRINT 

3  u^pyj 

PRINT! 

TMC-IIT 

'J  ( 

N 

3360 

GGSUe 

333S 

7*:'7Q 

so  TO 

4010 

^70)7 

v  -J  W  tl 

PRINT 

8Ni ^NONLINEAR  DATA  ANALYSIS  FOR  ORTHGTROPIC  iIGDEL  ^ 

3398 

PRINT 

cn  % 

3400 

PRINT 

8N; 

INPUT  DATA” 

3410 

PRINT 

iN;" 

3420 

PRINT 

iNL'LIN 

EAR  ELASTIC  rlODLiLI’ 

3430 

PRINT 

3440 

PRINT 

8N:'C11 

=  "iKil);’  iSPa!  311  =  LSil!;"  !GPa!'--l" 

3458 

pprs^iT 

8N:^C22 

=  “;K{2);”  (OPai  822  =  ’:3{2!i”  (SPa)-!’ 

3460 

pPTMT 

8Nj”C33 

=  ";Ki3!;’  (OPa)  S33  =  ’;S!3i;”  !SPa)"-l’ 

3478 

PRINT 

8N;^C44 

=  ";K!4);"  (SPa)  323  =  ";3!4);"  (OPa)"-!’ 

3471 

PRINT 

=  ”;K;5)r  iSPa)  313  =  ^3(5);”  (SPa)-!” 

3472 

PR]  N  T 

8N;^C66 

=  ”;K(6;;”  (SPa)  312  -  ”;3(6);”  (GPal-l” 

34S0 

PRINT 

3490 

PRINT 

8N; ^VELOCITY  DERIVATIVE  DATA  ISP3)"~r 

3500 

PRINT 

"vT  a 

-■  _■  i  V 

PRINT 

8Nj^L 

U\nn  -  t  \  .  M  f  1 

Hi.,..  p  .  \  i  J  A  1  ‘  ^  1  i  ;  J.  H 

3520 

PRINT 

”2. 

uoil  -  a.Lj/n  .  .  B  j./_  '"'1 

!  t  i.  i  .  “  J  !  1  \  i,  ,  1  t  .  =  ,  Jt  !'  ■  -A-  5  ^  ' 

3530 

PRINT 

aw,  ii  7 
•r:i,  1 

H322  =  ”;H(3,l)r  +/-  ”;H(3,2) 

3540 

PRINT 

H31 1  =  1 ) ;  ”  ^/-  ‘’;H(4,2} 

3550 

PRINT 

"5, 

q  i  7 1  -  B ,  u  f  r  1  i  .  9  ^  _  !! ,  l:  c  7 

n  i  i.  i  -  J  1  :  .  J  ,  i  /  4  7  ^ 

3560 

PRINT 

uj  y » ®  i 

H212  =  ”;H(6,1);”  -/-  ”;H(6,2) 

35T0 

PRINT 

"7, 

H312  =  ”;H(7,i);'  i-/-  ”;H(7,2) 

3588 

PRINT 

■IN:  "8. 

H321  =  ”;H(8,l)r  +/“  ’’;H(3,2] 

3598 

PRINT 

my-1. 

H123  =  ”;H(9,l}r  ”;H(9,2) 

3680 

PRINT 

iNs’ie. 

H213  =  ”]H(10jl);”  ”;H(10,2i 

3618 

PRINT 

'll. 

H323  =  'iHdUI';”  dH(il,2.) 

3620 

PRINT 

§N;”12. 

UT1  7  =  11  .  L!  ;  1  7  }  i  .  :i  ^/_  ^  /  t  ?  “  ; 

l  !  V  i  J  ^  J  i  ;  ^  ■  4  _  ^  . 

3638 

PRINT 

8N:M3, 

HI31  -  ’iH(13,l)r'  "/-  ■■;H(13,2! 

3640 

PRINT 

8N; ' 14. 

H232  =  ”;H(14d');''  dH(:4,2i 

■-=030 

PRINT 

§N:"15. 

H132  =  ";H(13,Ur  W-  7H(15,2; 

7  an 

-■  j  Li  a 

PRINT 

iN;”16, 

H231  =  ';H(16,1);"  ”;Ha6,2f 

-1"  -tci 

-=0  ;  0 

PRINT 

§N:’17. 

HOil  =  dH(17,I);”  W-  ^^;H(I7,2) 

7ipni 

-■■JUV 

PRINT 

!N:”18, 

H322  =  dHdS,!);”  ";HdS,2i 

3690 

PPTMT 

-N' '  19 

H012  =  ”;Hd9 jl) ;  ”  +/-  dHd9,2) 

3700 

CpTMT 

't  fi  i  ; 

8N:^2e. 

H021  =  ’;H(28,ii;”  ’;Hf28,2; 

37:0 

PPT>-iT 

«Nr21. 

H013  =  ’;K(2i,i!;’  t/-  “;H(21,2) 

3720 

PRINT 

■§N:  ”22. 

H323  =  ’;H(22,l);"  */-  '';H;22,2) 

T'-’Tn- 

j'ii 

PRINT 

8N: ”23. 

H031  =  ’;H(23.1!r  +/-  ’;H(23,2! 

3748 

PRINT 

8N;”24, 

H332  =  ";H(24,l!;’  r/-  ’;H{24,21 

3758 

PRINT 

3763  P!-!MT  3N: "THIRD  ORDER  ELASTIC  CONSTANTS  DESIvED  FRQH  LEAST’ 
3773  PRINT  «N:'S2uARES  FIT  AND  HAL  ERROR  AND  PRQB.  ERROR  (GPa)” 
3783  PRINT  SN;"  ” 
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3790  PRINT  ■aN:"Clll  =  'iCi!,!);’  +/-  ”;C(I,2);’  +/-  ";C(1,3I 

3800  PRINT  «N:'C112  =  ";C(2,1);'  +/-  ";C(2,2);'  +/-  ";C(2,3! 

3810  PRINT  §Nt'C113  =  ";C(3,1);'  +/-  '';C(3,2))*  +/-  ";C(3,3 

3820  PRINT  aN:"C122  =  ’;C(4,1)|*  +/-  ";C(4,2);’  +/- 

3830  PRINT  3N:’C222  =  ";C(5,1I;'  +/-  ';Ci5,2);*  +/- 

3840  PRINT  aN:'C223  =  "jCIi,!);'  +/-  "58(6, 2);’  +/-  ";C(t,31 

3850  PRINT  §N:*C144  =  ';C(7,li;’  +/-  ";Ci7,2)i"  +/- 


“;C!4,3i 

j  W  J  j  •-'  / 


3851  PRINT  3N:"C244  =  ";C(8,n; 
3352  PRINT  aN;*C344  =  ”;C(9,1);' 


■;C(7,31 
+/-  ’;C(8,2);"  +/-  ";C(8,3) 
+/-  ";C(9,2)i'  +/- 


";C!9,3) 


3353 

PRINT 

AN 

"C155  = 

■;C(10 

1) 

■  +/- 

’|C!10,2);" 

+  /' 

';C(10,3) 

3854 

PRINT 

m 

■C255  = 

";C(11, 

1) 

■;C(11,2);’ 

+  /- 

■;C(11,3) 

3855 

PRINT 

’C355  = 

1) 

“  +/- 

■;C(i2,2);’ 

+  /- 

’;C(12,3) 

3356 

PRINT 

m 

•C166  = 

";C(13, 

il) 

■  +/- 

■;C(13,2)i" 

tl- 

•;C(13,3) 

3857 

PRINT 

’C266  = 

";C!14, 

,1) 

■  +/- 

’;C(14,2);" 

+  /- 

”;C(14,3) 

3353 

PRINT 

’C366  = 

•;C(15, 

,1); 

"  +/- 

”;C(I5,2);’ 

+  /- 

’;C(15,3) 

33i0  PRINT  0N:"C133,C233,C333,C123,C456  UNDETERMINED  ' 

3370  PRINT  ^N:’  " 

3830  PRINT  aN: ’RECALCULATED  VELOCITY  DERIVATIVES  USING  FITTED  TDEC'S" 
3390  PRINT  SNi’AND  MAX.  AND  PROS.  ERRORS  IN  SAME  ORDER  AS  BEFORE’ 

3900  PRINT  AN:"  " 

3910  FOR  1=1  TO  24 

3'’23  PRINT  SN:!;".  ';H(!,3);’  +/-  "iH(I,4);'  +/-  ";H!I,5) 

3930  NEXT  I 
3940  PRINT  AN:"  " 

3950  PRINT  0N; "PERCENT  DIFFERENCES,  MAX.  AND  PRQB.  UNCERTAINTIES’ 
3960  PRINT  0N:’  ” 

3R^0  FOR  IM  TO' 24 

3930  PRINT  aN;I;".  *;V(I);T(I,1I;T(I,2) 

3990  NEXT  I 

3991  PRINT  AN:’  ' 

3992  PRINT  ?N: "CONDITION  #  FOR  INVERTING  MATRIX  IS:  ’;I3 
4030  RETURN 

4010  END 
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